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FOREWORD 


This  document  provides  a  description  of  the  near-Earth  natural  space  environment  tint  is 
recommended  for  space  vehicle  development  applications.  It  is  intended  to  provide  general  at 

the  conditions  space  systems  will  encounter  arid  must  successfully  withstand  during  their  on 

orbit  It  will  be  particularly  useful  for  preliminary  design  analyses.  While  these  data  and  models  are 
considered  to  be  accurate  and  the  best  available  for  general  engineering  applications,  their  snitabiliiy^ 
use  for  final  program  decisions  depends  upon  the  specific  design  problem  involved. 

The  natural  erivironment  is  characterized  by  many  complex  and  frequently  subtle  processes.  «na«y 
more  than  it  is  possible  to  treat  in  a  general  description  such  as  this.  In  many  cases  the  characteristics 
and  interactions  among  these  processes  are  poorly  understood  or  adequate  measurements  have  never 
been  made.  In  addition,  it  is  often  impossible  to  define  a  limiting  (e.g.,  maximum  possible) 
value  for  environmental  parameters.  Likewise,  it  may  not  be  technically  or  economically  feasible  lo 
design  a  system  to  withstand  an  extreme  value  when  it  can  be  defined  if  the  probability  is  that  ssh 
a  value  will  occur  during  the  mission  lifetime. 

For  these  reasons,  good  engineering  judgment  must  be  exercised  in  the  application  of  environmeit 
data  to  space  vehicle  design  analyses.  When  environmental  considerations  become  significau  Hpqfn  cr 
cost  drivers,  environmental  specialists  should  be  consulted  to  assure  that  the  environment  was  cMrectly 
imderstood  and  used  and  that  subtle  or  infrequent  effects,  those  not  addressed  by  the  general  purpose 
information  provided  here,  are  not  present  in  a  form  that  would  compromise  the  vehicle.  Assessments 
made  early  in  the  development  program  will  prove  advantageous  in  maintaining  an  economical  progtan 
and  obtaining  a  vehicle  with  minimal  operational  sensitivity  to  the  envirotunent  For  those  parameiers 
that  need  accurate  and  timely  monitoring  prior  to  or  during  operations,  this  early  planning  wfll  permit 
development  of  the  necessary  measurement  and  communication  systems. 

This  document  is  a  follow-on  to  NASA  TM-82478.  Space  and  Planetary  Environment  Criteria 
Guidelines  for  Use  in  Space  Vehicle  Development.  1982  Revision  (Volume  I),  January  1983,  whkii  is 
still  useful  for  many  applications.  In  case  of  conflict,  the  data  in  this  document  should  be  used.  There  is 
no  intent  to  automatically  replace  references  to  NASA  TM-82478  or  other  previous  in  ay 

contract  Scope  of  Work  by  the  issuance  of  this  document 

(Questions  or  requests  for  assistance  in  the  application  of  these  natural  environment  gnufelines  to 
aerospace  vehicle  programs  should  be  addressed  to  the  Chief.  Electromagnetics  and  Enviromients 
Branch,  EL54,  Systems  Analysis  and  Integration  Laboratory,  Science  and  Engineering  Duectorate. 
Marshall  Space  Flight  Center.  AL  35812.  s  s 
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TECHNICAL  MEMORANDUM 

NATURAL  ORBITAL  ENVIRONMENT  GUIDELINES 
FOR  USE  IN  AEROSPACE  VEHICLE  DEVELOPMENT 

I.  SCOPE  AND  PURPOSE 


The  purpose  of  this  document  is  to  provide  definitions  of  the  natural  near-Earth  space  environ¬ 
ment  suitable  for  use  in  the  initial  development/design  phase  of  any  space  vehicle.  The  natural  environ¬ 
ment  includes  the  neutral  atmosphere,  plasma,  charged  particle  radiation,  electromagnetic  radiation 
(EMR),  meteoroids,  orbital  debris,  magnetic  field,  physical  and  thermal  constants,  and  gravitational 
field.  Communications  and  other  unmanned  satellites  operate  in  geosynchronous  Earth  orbit  (GEO), 
therefore,  some  data  are  given  for  GEO,  but  emphasis  is  on  altitudes  from  200  km  to  1000  km  (low- 
Earth  orbit  (LEO)). 

This  document  covers  neither  the  environment  below  100  km  altitude  (see  NASA  TM-45 1 1 , 
Terrestrial  Environment  (Climatic)  Criteria  Guidelines  for  Use  in  Aerospace  Vehicle  Development, 

1993  Revision)  nor  the  induced  environment  and  other  effects  resulting  from  the  presence  of  the  space 
vehicle  at  orbital  altitudes.  Man-made  factors  are  included  as  part  of  the  ambient  natural  environment. 
i.e.,  orbital  debris  and  radio  frequency  (RF)  noise  generated  on  Earth,  because  they  are  not  caused  by  the 
presence  of  the  space  vehicle  but  form  part  of  the  ambient  environment  that  the  space  vehicle  experi¬ 
ences.  It  is  very  important  to  take  induced  environments  into  account  because  the  net  effect  can  be  quite 
different  from  the  unperturbed  ambient  natural  environment  This  document  does  not  provide  techniques 
or  engineering  solutions  to  permit  operation  in  the  natural  environments  described  herein. 


1.1  Format  and  Use  of  the  Document 

The  format  of  this  document  is  such  that  each  section  contains  an  explanation  and  description  of 
the  natural  environment  pnenomenon  to  which  it  is  devoted.  In  general,  the  descriptions  include  the 
mean  and  limiting  values  of  each  parameter.  It  is  intended  that  these  data  will  be  sufficient  for  most 
space  vehicle  general  design  purposes;  however,  for  systems  which  preliminary  analyses  show  are  sensi¬ 
tive  to  the  environmental  parameter  being  used,  users  should  contact  personnel  of  the  Electromagnetics 
and  Environments  Branch,  ELS4,  Marshall  Space  Flight  Center  (MSFC),  AL  3S812.  For  additional 
information,  the  user  is  referred  to  the  reference  section  of  this  document.  Background  material  appears 
in  references  1, 2,  and  3. 
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n.  GENERAL  INFORMATION 


This  section  provides  frequently  used  physical  constants  and  describes  the  general  character  of 
terrestrial  space. 


2.1  Sun-Earth  Constants 

The  values  given  in  table  2-1  define  the  Sun-Earth  parameters  for  aerospace  vehicle  design  per¬ 
formance  analyses. 


Table  2-1.  Sun-Earth  physical  constants. 

1.4959787E*  km  (equals  1  AU  by  definition) 
1.37i:,‘o  W/ra2  at  1  AU 


Distance  to  the  Sun 
Solar  Constant  (Note  1 ) 

Eccentricity  of  Orbit  (Note  2) 

Orbital  Period  (Sidereal)  (Note  3) 

Radius  of  the  Earth  (equatorial) 

Mass  of  the  Earth 
Earth  Rotation  Rate 

Gravitational  Constant  for  the  Earth  dig) 
Inclination  of  the  Equator  (Note  4) 
Period  of  Rotation  (Sidereal) 

Space  Sink  Temperature 
Solar  Radiation  Pressure  at  1  AU 


=:  0.0167295 
=  365.25636  days 
=  6378.140  km 
=  5.977E24kg 
=  0.72921E-4(rads)/s 
=  3.986012EMNm2/kg 
=  23.45“ 

=  23.934  h  (86  162.4  s) 

=  3  K  (absolute) 

=  9.02xE-^N/m2  (100-percent  reflecting) 


Notes;  1.  The  spectrum  of  the  Sun  is  a  blackbody  spectrum  with  a  characteristic  temperature  of 

5762  K.  The  -i-lO  and -10  W/m^  are  due  to  the  natural  variability  of  the  solar  output  and 
measurement  uncertainty. 


2.  The  eccentricity  gives  the  noncircular  nature  of  an  orbit.  The  maximum  distance  firom  the  Sun  is 
(1+ecc)  times  tlie  average  radius;  the  minimum  distance  is  (1-ecc)  times  the  average  radius. 

3.  The  sidereal  period  is  measured  with  respect  to  the  “fixed”  stars  rather  than  with  reqiea  u)  the  Sun. 

4.  The  inclination  of  the  Equator  is  with  respect  to  the  Earth’s  orbital  plane. 


CAUTION:  The  above  are  accepted  values  for  general  use.  However,  the  various  gravitational  and  other 
models  may  have  been  developed  using  other  values.  When  this  is  ttie  case,  the  values 
associated  with  the  model  should  be  used. 
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2J1  Terrestrial  Space 


of  space  extends  from  the  base  of  the  ionosphere  at  about  60  km  above  the  surface  of 

c  boundary  of  the  magnetosphere  beyond  which  interplanetary  space  is  unaffected  by  the 

Earth.  This  distance  is  about  95  000  km  above  the  surface  of  the  Earth  (16  radii  of  the  Earth  (Re))  in  the 
sunward  direction  and  several  limes  this  in  the  anti-sunward  direction.  This  region  is  loosely  referred  to 
as  the  magnetospltere,  although  more  strictly  speaking,  this  term  means  the  (major)  part  of  terrestrial 
Space  into  which  the  Earth *$  magnetic  field  extends. 


The  morphology  is  roughly  axisymmetric  within  4  Re  of  the  Earth’s  center,  but  at  greater  dis¬ 
tances  It  becomes  wty  unsymmetric,  with  a  long  tail  extending  in  the  anti-sunwani  direction.  The  orin- 
cipal  regions  ana  their  interacting  phenomena  are  described  below  and  illustrated  in  figure  2-1. 


SOLAR 

WIND 


Figure  2-1.  Schematic  view  of  terresuial  space.'* 

f  Pislds-  The  gravitational  field  results  from  the  mass  of  the  solid  Earth  and  reflects  the 
distobuuon  of  that  mass.  It  traps  the  neutral  atmosphere,  consuains  its  motion,  and  influences  the 
motions  of  meteoroids  and  debris.  However,  it  has  litUe  effect  on  the  rest  of  terrestrial  space  because 
electnc^  forew  arc  »  much  stronger.  The  magnetic  field  has  two  sources:  (1)  currents  Lide  the  Earth 
at  produce  about  99  percent  of  the  field  at  the  surface  and  (2)  currents  in  the  magnetosjAere.  The  latter 
becomes  relanvely  more  important  beyond  a  few  Earth  radii  because  the  internal  field  deaeases  as  the 
reverse  distance  cubed  from  the  Earth’s  center.  me 

Earth’s  field  may  be  regarded  as  a  dipole  tilted  1 1.7®  from  the  rotation 
axis  and  offset  from  the  geometnc  center  of  the  Earth  by  430  km  in  the  direction  of  southeast  Asia. 


“ZT" 
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Many  phenomena  arc  related  to  magnetic  latitude  which,  as  a  result  of  the  tilt,  is  1 1.7°  greater  than 
geographic  latitude  in  the  longitude  of  eastern  North  America  and  1 1.7°  less  on  the  opposite  side  of  the 
world.  The  offset  puts  the  surface  of  the  Earth,  or  a  circular  orbit,  at  a  higher  altitude  with  re.spect  to  the 
geomagnetic  Held  in  the  region  of  the  South  Atlantic  off  the  coast  of  Brazil  than  it  is  elsewhere.  This 
region  is  called  the  South  Atlantic  Anomaly.  Since  both  the  tilt  and  the  offset  are  changing  slowly,  the 
South  Atlantic  Anomaly  is  drifting  slowly  to  the  west.  For  some  applications,  it  may  be  necessary  to 
include  the  magnitude  and  direction  of  this  drift.  If  it  is  necessary,  contact  the  personnel  of  the  Electro¬ 
magnetics  and  Environments  Branch,  ELS4,  Marshall  Space  Flight  Center,  AL  35812. 

There  are  also  electric  fields  that  result  from  the  motion  of  tte  magnetospheric  plasma.  In  the 
ionosphere,  these  electric  fields  are  perpendicular  to  the  magnetic  field  and  have  magnitudes  up  to  0.1  to 
0.5  V/m. 

2.2.2  Plasma  and  Particles.  The  neutral  atmosphere  extends  from  the  surface  to  2500-km  altitude 
01  more,  where  its  density  has  dropped  to  approximately  lO-*’  kg  m-3.  Its  density  continues  to  dpc^ease 
at  greater  altitudes,  so  an  outer  limit  is  not  rigorously  defined. 

Sunlight  and,  to  a  lesser  extent,  fast  charged  particles  ionize  some  of  the  neutral  atmosphere, 
creating  a  plasma  consisting  of  equal  number  densities  of  ions  and  electrons.  In  the  altitude  range  of 
about  60  to  1000  km.  this  is  called  the  ionosphere.  The  number  density  of  the  plasma  is  less  than  that  of 
the  neutral  atmosphere  below  about  1000  km.  However,  collisions  between  plasma  and  nftntraig  ate 
infrequent  enough  above  about  150  km  that  the  two  ate  decoupled  and  electrical  forces  dominate  the 
behavior  of  the  ionosphere  and  other  plasma.  Between  about  60  and  150  km.  the  neutral  atmosphere  and 
ionosphere  are  coupled  in  a  complicated  way. 

Plasma  fills  the  rest  of  terrestrial  space  as  well  and  is  given  various  names  in  different  regions, 
"nte  region  out  to  the  field  lines  at  (L  =  4.5)  is  sometimes  referred  to  as  the  outer  ionosphere  and  some¬ 
times  as  the  plasmasphere.  Beyond  that,  the  plasma  density  declines  to  its  interplanetary  value  of  about 
10^  ions  m-3.  The  plasma  is  electrically  neutral  (equal  ion  and  electron  number  densities),  and  the 
energy  distribution  of  its  particles  may  be  described  by  one  or  more  temperatures  (Maxwellian  distribu¬ 
tions).  The  plasma  may  have  a  bulk  streaming  velocity  as  well. 

In  addition,  there  ate  fluxes  of  fast  moving  particles,  such  as  the  trapped  radiation  (Van  Allen 
belts),  the  aurora,  and  cosmic  rays  that  have  very  nonthermal  energy  distributions.  These  are  called 
penetrating  charged  particles.  The  shape  of  the  magnetosphere,  the  motion  of  the  plasma,  and  the 
acceleration  of  the  fast  particles  ate  caused  by  the  inteiactioos  of  the  Earth’s  magnetic  field  rotating  with 
the  Earth,  the  ionosphere,  and  the  solar  wind  plasma  flowing  past  the  magnetic  field.  The  solar  wind  is  a 
very  tenuous  plasma  that  flows  radially  outward  from  the  Sun  through  interplanetary  space. 

2.2.3  Variability.  Processes  within  terrestrial  space  are  partially  controlled  by  the  level  of  solar 
activity.  The  solar  activity  varies  more  or  less  cyclically  with  an  average  period  of  1 1  years.  The 
electromagnetic  radiation  (EMR)  emitted  by  the  Sun  varies  (although  not  much  in  the  visible  portion  of 
the  spectrum)  as  does  the  solar  wind,  the  solar  magnetic  field,  and  the  production  of  solar  cosmic  rays. 
The  exact  level  of  solar  activity  cannot  be  predicted  very  accurately,  although  the  phase  within  the  1 1- 
year  period  can  be  established.  Figure  2-2  shows  solar  activity  from  17(X)  to  1992.  In  addition,  plasma, 
radio  noise,  and  energetic  particles  tend  to  be  emitted  from  localized  regions  on  the  Sun’s  surface.  These 
localized  active  regions  and  some  coronal  features  persist  longer  than  the  solar  rotation  period  of  27 
days,  and  since  they  only  aftect  the  Earth  when  they  face  it,  enhanced  solar  activity  can  be  estimated  27 
or  mote  days  in  advance. 


Figure  2-2.  Solar  cycle  as  represented  by  yearly  mean  sun^t  number  for  the  period  1700 to  1992. 
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m.  NEUTRAL  ATMOSPHERE 


The  state  of  the  neutral  aunospherc  is  most  conveniently  described  in  tenns  of  a  mea.  wnb 
spatial  and  temporal  variations  about  that  mean.  For  space  vehicle  operations,  the  neutral  atmospteie  is 
significant  because  (1)  even  at  its  low  density,  it  produces  torques  and  drags  on  the  vehicle;  (2) 
density  height  profile  of  the  atmosphere  above  100  km  altitude  modulates  the  flux  of  tn^tpediadiasioa 
encountered,  as  explained  in  section  V  of  this  dccument,  and  the  orbital  debris,  as  explained  in  seodon 
VII  of  this  document;  and  (3)  the  atomic  oxygen  both  erodes  and  chemically  changes  those 
which  are  exposed  to  it. 


3.1  Thermosphere  Region  ' 

The  region  of  the  Earl’s  atmosphere  lying  between  about  90  and  500  km  is  known  as  the 
thermosphere,  while  that  region  lying  above  500  km  is  known  as  the  exosphere.  The  temperamre  n  tlK 
lower  thermosphere  increases  rapidly  with  increasing  altitude  from  a  minimum  at  90  km.  Eveatu^  it 
becomes  altitude  independent  at  upper  thermospheric  altitudes.  This  asymptotic  temperature,  known  as 
the  exospheric  temperature,  is  consuint  with  altitude  due  to  the  extremely  short  thermal  conducticmtiine. 

The  thermospheric  gases  are  heated  by  the  absorption  of  the  solar  extreme  ultraviolet  (EUVJ 
radiation.  At  the  lowest  thermosphe;ic  altitudes,  the  absorption  of  ultraviolet  (UV)  radiation  is  also 
important  The  EUV  and  UV  radiation  initially  heats  only  the  dayside  thermo^heic,  and  although  con¬ 
ductive  and  convective  processes  att  to  redistribute  some  of  this  energy,  a  large  temperature  gtadicDt 
always  exists  between  the  daytime  ;ind  the  nighttime  thermosphere.  An  average  daytime  exospheric 
temperature  is  1060  K,  and  an  aver.ige  nighttime  exospheric  temperature  is  840  K.  The  longiindiniri 
temperature  gradient  causes  a  wind  to  fiow  from  the  dayside  to  the  nightside  thermosphere  with  sgeeds 
typically  reaching  100  m/s.  ’ 


An  additional  heat  source  for  the  thermosphere  is  the  interaction  of  the  Earth’s  magnetic  at 
very  great  distances  (several  Earth  radii),  in  the  region  known  as  the  magnetopause,  with  the  sidar  wind. 
The  solar  wind  is  a  stream  of  high  speed  plasma  emanating  from  the  Sun.  This  interaction  caa- 
getic  particles  to  penetrate  down  into  the  lower  thermosphere  at  high  geographic  latitudes  and  direedy 
heat  the  thermospheric  gas.  Thesi;  energetic  particles  are  also  responsible  for  the  aurwa  seen  « 
high  latitudes.  In  addition,  electric  fields  mapped  down  from  the  magnetosphere  onto  the  high  latitnde 
ionosphere  cause  electric  currents  to  flow.  The  ionosphere  is  a  small  fraction  of  the  thermosphere  dot 
remains  ionized  due  to  the  solar  radiation.  It  never  totally  disappears  at  night,  and  during  dayibfat  hniirc 
the  ionization  density  never  exceeds  more  than  1  percent  of  the  neutral  density.  These  electrical  cuoenis 
lose  enerjgy  through  ohmic  or  joule  dissipation  and  heat  the  neutral  thermospheric  gas.  The  ions  ai«» 
collide  directly  with  the  neutral  gas,  setting  the  whole  gas  into  motion.  At  these  high  latitudes,  the  wand 
speeds  generated  by  this  procef;s  can  be  very  large,  at  times  as  large  as  1.5  km/s.  Eventually  viscous 
effects  dissipate  these  winds,  and  their  lost  kinetic  energy  provides  an  additional  heat  somce  for  the 
neutral  thermospheric  gas. 

The  high-latitude  heat^sources  are  effective  during  both  the  day  and  night  Although  an  inteiF- 
mittent  source  of  energy  for  t|ie  thermosphere,  they  can  at  times  exceed  the  global  EUV  energy  absmbed 
by  the  thermosphere.  In  addition,  although  the  energy  is  deposited  at  high  latitudes  (greater  6flPor 
so),  the  disturbance  effects  at  e  transmitted  to  lower  latitudes  through  the  actions  of  winds  and  waves. 
However,  the  disturbance  efi'ects  at  low  latitudes  are  significantly  smaller  than  they  are  at  higlrr 
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latitudes.  The  high-latitude  ionospheric  currents  that  flow  perturb  the  geomagnetic  fleld,  so  that  such 
disturbances,  which  can  be  detected  by  ground-based  magnetometers,  are  referred  to  as  geomagnetic 
storms. 


Whenever  the  neutral  thermospheric  gas  is  heated,  it  expands  radially  outward.  Because  the 
undisturbed  thermospheric  density  decreases  with  increasing  altitude,  an  outward  expansion  of  the  gas 
results  in  an  increase  of  density  at  high  altitudes.  Thus,  the  daytime  thermospheric  density  is  greater  tkan 
the  nighttime  density,  while  during  times  of  geomagnetic  storms,  the  high-latitude  densi^  is  greater  than 
it  is  during  undisturbed  periods.  This  anisotropic  heating  leads  to  the  so-called  diurnal  and  polar  bulges, 
which  were  first  inferred  from  the  increased  drag  experienced  by  orbiting  satellites. 

Below  the  turbopause  (located  at  about  lOS  km  altitude),  the  atmosphere  is  weU  mixed  by  turbs- 
lence,  so  that  the  composition  of  the  atmosphere  does  not  vary  with  altitude.  Above  the  turbopause, 
however,  diffusion  becomes  so  rapid  that  the  altitude  variation  of  the  various  species  becomes  dependent 
on  molecular  mass,  with  the  result  that  composition  varies  with  altitude.  Thus,  the  number  densities  of 
the  heavier  thermospheric  species  (N2  and  62)  decrease  with  increasing  altiuide  much  faster  than  those 
of  the  lighter  species  (H  and  He).  This  means  that  the  heavier  molecular  qtecies  predominaie  in  the 
lower  thermosphere,  while  the  lighter  atomic  species  predominate  in  the  upper  thermosphere.  A  typical 
altitude  profile  for  the  individual  thermospheric  constituents  is  shown  in  figure  3-1.  Lifting  of  the 
thermosphere  will  cause  the  mean  molecular  weight  at  a  given  altitude  to  increase,  while  a  sinking 
motion  will  cause  it  to  decrease. 


3.2  Variations 

In  addition  to  the  diurnal  variation  in  the  neutral  mass  density  at  orbital  altitudes,  there  is  a  serru- 
annual  variation,  a  seasonal-latitudinal  variation  in  the  lower  thermosphere,  and  a  seasonal-latitudinal 
variation  in  the  helium  number  density  as  well  as  nonperiodic  perturbations  associated  with  variations  in 
solar  activi^,  atmospheric  waves,  and  thermospheric  winds.  All  of  these  variations  are  discussed  in  the 
paragraphs  that  follow. 

3.2.1  Variation.s  With  .Solar  Activity.  The  short  wavelength  solar  electromagnetic  radiation 
(EUV  and  UV)  changes  substantially  with  the  overall  level  of  solar  activity,  with  the  result  that  the 
thermospheric  density,  especially  at  orbital  altitudes,  is  strongly  dependent  on  the  level  of  solar  activityr. 
Thus,  there  is  an  1 1-year  variation  in  the  thermospheric  mass  density,  corresponding  to  the  mean  1 1-ycar 
variation  in  solar  activity.  Similarly,  there  is  also,  on  the  average,  a  27-day  variation  in  density  that  is 
related  to  the  mean  27-day  solar  rotation  period,  although  the  variation  tends  to  be  slightly  longer  than 
27  days  early  in  the  cycle  when  active  regions  occur  more  frequently  at  higher  latitude  and  slightly 
shorter  than  27  days  later  in  the  cycle  when  active  regions  occur  more  frequently  closer  to  die  Sun’s 
equator.  The  appearance  of  coronal  holes  and  active  longitudes  also  affects  this  average  27-day  varia¬ 
tion.  Changes  in  the  thermospheric  density  related  to  changes  in  the  solar  output  during  active  periods 
associated  with  flares,  eruptions,  coronal  mass  ejections  (CMEs)  and  coronal  holes  (CHs)  can  begin 
almost  instantaneously  (minutes  to  hours),  although  mote  often  a  day  or  mote  lag  is  seen.  Hgure  3-2 
shows  typical  neutral  densities  for  periods  of  high  and  low  solar  activity. 

3.2.2  Variations  During  Periods  of  Increased  Solar  Activity  (Flares.  Coronal  Mass  Erections. 
and  Coronal  Holesl.  As  previously  described,  during  periods  of  increased  solar  activity,  the  enhanced 
interaction  of  the  solar  wind  with  the  Earth’s  magnetosphere  leads  to  a  high-ladtude  heat  and  momenm 
source  for  the  thermospheric  gases.  Some  of  this  heat  and  momentum  is  convected  to  low  latitudes.  This 
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episodic  type  of  increased  solar  activity  that  usually  causes  variations  in  the  Earth’s  magnetic  Add 
varies  over  the  solar  cycle  and  usually  has  two  or  more  major  peaks,  one  during  the  rise  of  the  cyck  and 
the  other  with  Ijtfger  peaks  during  the  decline  of  the  cycle.  Also,  more  intense  solar  cycks  seem  to  haw 
more  intense  episodic  type  activity.  Finally,  there  is  a  seasonal  variation  associated  with  this  episodic 
type  activity  with  the  density  usually  being  greatest  in  March  (±1  month)  and  September  (±1  momh)  of 
each  year.  This  variaUon  is  possibly  related  to  the  path  of  the  Earth  in  its  rotaUon  around  the  Sun. 


Using  MSFC/MET  Number  Density  (cm’®) 

Figure  3-1.  Number  density  of  atmospheric  constituents  versus  altitude. 
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Figuie  3-2.  Typical  atmospheric  mass  density  profiles  at  high  and  low  solar  activity. 

(For  information  only;  not  for  design  use.) 

3.2.3  The  Diurnal  Variation.  The  rotation  of  the  Earth  with  respect  to  the  solar  EUV  heat  source 
induces  a  diurnal  (24-hour  period)  variation  (or  diurnal  tide)  in  the  thermospheric  temperature  and  den¬ 
sity.  Due  to  a  lag  in  the  response  of  the  thermosphere  to  the  EUV  heat  source,  the  density  at  orbital 
altitudes  maximizes  around  2  p.m.  local  solar  time  at  a  latitude  approximately  equal  to  that  of  tbe  sub- 
solar  point.  The  lag,  which  is  a  function  of  altitude,  decreases  with  decreasing  altitude.  Similarly,  the 
density  minimum  occurs  between  3  and  4  a.m.  local  solar  time  at  about  the  same  latitude  in  the  oppoate 
hemisphere.  In  the  lowest  regions  of  the  thermosphere  (120  km  and  below)  where  the  characteristic 
thermal  conduction  time  is  on  the  order  of  a  day  or  more,  the  diurnal  variation  is  not  a  predominant 
effect 
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The  various  constituents  of  the  thermosphere  do  not  a3  respond  to  the  diurnal  variation  of  the 
solar  EUV  heat  source  with  the  same  amplitude  and  phase.  The  time  lag  is  longer,  by  as  much  as  2  hours 
at  orbital  altitudes,  for  the  heavier  constituents  (Na.  62,  and  Ar)  than  for  O.  By  contrast,  the  lighter 
species  number  densities  maximize  in  the  early  morning  hours  (3  a.m.  and  7  a.m.  local  solar  time,  for  H 
and  He.  respectively).  This  is  due  to  dynamic^  (buoyancy)  effects. 

Harmonics  of  the  diurnal  tide  are  also  induced  in  the  Earth’s  atmosphere.  In  particular,  a  semi¬ 
diurnal  tide  (period  of  12  h)  and  a  terdiumal  tide  (period  of  8  h)  are  important  in  the  lower  themosphere 
(below  some  160  km  for  the  semidiurnal  tide  and  much  lower  for  the  terdiumal  tide).  These  tides  are  not 
important  at  orbital  altitudes. 

3.2.4  Semiannual  Variation.  This  variation  is  believed  to  be  a  conduction  mode  of  oscillation 
driven  by  a  semiannual  variation  in  joule  heating  in  the  high-latitude  thermosphere  (because  of  a 
semiannual  variation  in  the  effects  on  the  Earth-atmosphere  system  of  episodic-type  variations  in  solar 
activity).  The  variation  is  latitudinally  independent  and  is  modified  by  composition  effects.  The 
amplitude  of  the  variation  is  height-dependent  <30  to  20  percent  and  asymmetric  about  the  mean  density 
at  350  km)  and  variable  from  year  to  year,  with  a  primary  minimum  in  July,  primary  maximum  in 
October,  and  a  secondary  minimum  in  January  followed  by  a  secondary  maximum  in  April.  It  has  been 
found  that  the  magnitude  and  altitude  dependence  of  the  semiannual  oscillation  vary  considerably  from 
one  solar  cycle  to  the  next.  This  variation  is  important  at  orbital  altitudes. 

3.2.5  Seasonal-Latitudinal  Variations  of  the  Lower  Thermosphere  Density.  This  variation  is 
driven  in  the  thermosphere  by  the  dynamics  of  the  lower  atmosphere  (mesosphere  and  below).  The 
amplitude  of  the  variation  maximizes  in  the  lower  thermosphere  somewhere  between  about  105  and  120 
km,  diminishing  to  zero  at  altitudes  around  170  km.  Although  the  temperature  osciOaiicn  amplitude  is 
quite  large,  the  corresponding  density  oscillation  amplitude  is  small.  This  variation  is  not  important  at 
orbital  altitudes. 

3.2.6  Seasonal-Latitudinal  Variations  of  Helium.  Satellite  mass  spectrometers  have  measured  a 
strong  increase  of  helium  above  the  winter  pole.  Over  a  year,  the  helium  number  density  varies  by  a 
factor  of  42  at  275  km,  12  at  400  km,  and  3  or  4  above  500  km.  The  formation  of  this  winter  helium 
bulge  is  primarily  due  to  the  effects  of  global  scale  winds  that  blow  from  the  summer  to  the  winter 
hemisphere.  The  am.piitude  of  the  bulge  decreases  with  increasing  levels  of  solar  activity,  due  to  the 
increased  effectiveness  of  exospheric  transport  above  500  km  that  carries  helium  back  to  the  summer 
hemisphere.  There  is  also  a  very  weak  dependence  of  the  helium  bulge  amplitude  on  the  magnitude  of 
the  lower  thermospheric  eddy  diffusivity. 

3.2.7  Thermospheric  Waves.  Ructuations  have  been  detected  in  temperature  and  density 
measurements  throughout  the  atmosphere  from  the  ground  up  to  at  least  510  km.  Some  of  these 
fluctuations  are  caused  by  gravity  waves,  so  named  because  they  are  primarily  oscillations  of  the  neutral 
gas  for  which  the  restoring  force  is  gravity.  A  thermospheric  gravity  wave  produces  a  corresponding 
wave  in  the  ionosphere  known  as  a  traveling  ionospheric  disturbance. 

Thermospheric  gravity  waves  oscillate  with  periods  typically  in  the  range  of  30  min  to  several 
hours,  and  have  horizontal  wavelengths  in  the  range  of  hundreds  of  kilometers  up  to  about  4000  km.  The 
density  amplitudes  of  the  larger  scale  waves  are  larger  at  higher  latitudes  and  diminidi  toward  the 
equator.  At  about  200  km  altitude,  typical  values  of  these  amplitudes  are  15  percent  of  the  mean  at 
auroral  latitudes  and  5  percent  of  the  mean  at  equatorial  latitudes.  The  smaller  scale  waves  have  ampli¬ 
tudes  that  are  essentially  independent  of  latitude.  Gravity  wave  amplitudes  genet .  (iy  decrease  at  greater 
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altitudes  in  the  thermosphere  due  to  dissipation  by  molecular  processes.  The  larger  scale  waves  survive 
to  greater  altitudes  than  do  the  smaller  scale  waves. 

3.2.8  Thermospheric  Winds.  Figures  3-3  and  3-4  show  die  general  flow  patterns  of  thermo¬ 
spheric  winds  between  100  and  700  km,  as  they  are  currently  known.  The  wind  speeds  range  from  100 
to  200  m/s  at  low  latitudes  (less  than  28.5*),  while  at  high  latitudes  (greater  than  about  65*)  they  can  be 
as  large  as  1500  m/s  or  more.  Rapid  (minutes)  changes  in  wind  direction  (of  up  to  180*),  probably  driven 
by  gravity  waves,  have  also  been  observed. 

3.2.8. 1  Austral  (South  Polar  Region)  Thermospheric  Winds.  Ground-based  observations  with 
Fabry-Perot  interferometers  have  shown  that  there  is  a  direct  relationship  between  austral  thermospheric 
winds  at  F-layer  altitudes  and  the  interplanetary  magnetic  field  (IMF),  particularly  in  the  midnight  mag¬ 
netic  local  time  (MLT)  sector.  The  strength  of  the  zonal  component  of  the  austral  neutral  wind  in  the 
MLT  midnight  sector  is  directly  correlated  with  the  strength  and  direction  of  the  By  component  of  the 
IMF.  The  zonal  wind  in  the  MLT  midnight  sector  changes  direction  when  the  sign  of  the  By  component 
of  the  IMF  changes,  and  the  speed  of  the  zonal  wind  in  the  MLT  midnight  sector  is  directly  proportional 
to  the  strength  of  the  By  componenu  Conversely,  the  meridional  wind  in  the  MLT  midnight  sector  is  cor¬ 
related  with  the  B,  component  of  the  IMF.  The  meridional  wind  decreases  as  the  B^  component  turns 
northward.  When  B^  is  positive,  sunward  winds  develop  in  the  MLT  midnight  sector  and  exceed  the 
strength  of  the  zonal  winds.^ 

3.2.9  Thermospheric  Tides.  Models  of  the  tides  in  the  lower  thermosphere  (90  to  140  km)  show 
the  midlatitude  semidiurnal  tide  in  the  wind  has  an  amplitude  of  tens  of  m/s,  maximizing  at  50  to  70  m/s 
around  1 10  km.  Temperature  tide  amplitudes  typically  reach  a  maximum  of  10  to  20  K;  however,  the 
models  do  not  always  agree  with  observations  since  a  semidiurnal  oscillation  of  about  70  K  has  been 
measured  at  low  latitudes  at  120  km  altitude.  The  best  agreement  occurs  at  midlatitudes,  but  in  general, 
the  models  under  estimate  both  the  wind  and  temperature  tidal  amplitudes.  Phase  differences  are  rela¬ 
tively  small  for  the  temperature  tide,  but  are  large  for  the  tidal  winds,  as  much  as  3  h. 

Radar  measurements  of  midlatitude  semidiurnal  wind  tides  show  maximum  amplitudes  of  60  to 
70  m/s  in  the  110  to  120  km  altitude  range.  Values  of  the  maximum  midlatitude,  annual  mean  semi¬ 
diurnal  amplitudes  and  the  altitudes  at  which  they  occur  are:  (a)  temperature,  25  K  at  1 15  km;  (b)  north¬ 
ward  wind,  50  m/s  at  125  km  altitude;  and  (c)  eastward  wind,  50  m/s  at  125  km  altitude.  The  local  times 
of  these  maxima  generally  decrease  with  increasing  altitude;  temperature  leads  the  meridional  wind  by 
about  6  h,  whUe  the  zonal  wind  leads  the  temperature  by  a  few  hours. 

At  middle-thermospheric  altitudes  (300  km),  the  mean  winds  are  poleward  in  winter  and 
equatorward  in  summer,  llie  amplitude  of  the  diurnal  component  of  the  mean  wind  is  much  greater  at 
solar  minimum  than  at  solar  maximum,  showing  that  the  influence  of  ion  drag  more  than  compensates 
for  the  greater  EUV  input  at  solar  maximum.  A  typical  seasonally  averaged  wind  amplitude  for  the 
diurnal  tide  is  about  50  m/s.  The  mean  semidiurnal  wind  tide  amplitude  averaged  over  all  seasons  and 
solar  cycle  conditions  is  about  41  m/s. 

3.2.9. 1  Solar  Storm  Induced  Effects  in  the  Thermosphere.  During  magnetic  .storms  and  .sub- 
storms,  pardcle  precipitation  and  the  joule  dissipation  of  the  high-latitude  ionospheric  currents  heat  the 
neutral  gas  at  altitudes  above  some  1 10  km.  The  thermosphere  then  expands  and  densities  above  about 
120  km  increase.  The  energy  deposited  in  the  high-latitude  thermosphere  at  such  times  is  significant  and 
can  exceed  the  global  EUV  absorbed  by  the  thermosphere.  It  remains  an  important  beat  source  locally  at 
night  in  the  auroral  zones. 


Latitude  (degrees) 


Hgure  3-3.  Schematic  diagram  of  the  zonal  mean  meridional  circulation  in  the  thermosphere 
at  equinox  for  various  levels  of  magnetic  (auroral)  activity (Maximum  velocities  are 
typically  hundreds  of  meters  per  second.) 


Latiaide  (degrees) 


Hguie  3-4.  Schematic  diagram  of  the  zonal  mean  meridional  circulation  in  the  thermosphere 
at  solstice  for  various  levels  of  magnetic  (auroral)  activity.^  (Maximum  velocities 
are  typically  hundteds  of  meters  per  second.) 
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The  energy  generated  during  a  substorm  is  transferred  from  high  to  low  latitudes  principally  by 
the  meridional  winds  assisted  by  gravity  waves.  This  high-latitude  beating  sets  up  a  circulation  cell  in 
which  gas  is  lifted  in  the  auroral  zone,  flows  equatorward,  sinks  at  middle  to  low  latitudes,  and  then 
returns  to  the  heating  region  at  low  altitudes.  The  sinking  causes  compressional  heating  (figs.  3-3  and 
3-4). 


Since  the  mean  molecular  weight  decreases  with  altitude,  lifting  of  the  gas  at  auroral  latitudes 
leads  to  increases  in  the  mean  molecular  weight  at  a  fixed  altitude.  Compositional  changes  also  accom¬ 
pany  geomagnetic  storms  and  these  changes  are  even  observed  at  middle  and  low  latitudes.  In  general 
during  a  geomagnetic  storm,  there  is  a  large  depletion  in  O  and  an  enhancement  in  O2  within  the  auroral 
ovals.  The  magnitude  of  both  the  O  depletion  and  the  O2  enhancement  is  underestimated  in  current 
models,  with  the  misrepresentation  increasing  with  the  use  of  daily  values  of  Ap  instead  of  the  3-hourly 
values.7  Analyses  which  require  instantaneous  values  of  either  atomic  or  molecular  oxygen  within  the 
auroral  ovals  should  take  special  precautions  in  using  the  available  models  of  the  thermosphere. 

In  one  numerical  simulation  of  a  geomagnetic  substorm,  it  was  found  that  at  both  200  and 
450  km  the  latitudinal  variations  in  the  mass  density  during  the  substorm  were  relatively  insignificant 
and  did  not  at  all  reflect  the  total  heat  input  into  the  auroral  zone.  In  particular,  it  was  apparent  that  the 
global  density  response  was  relatively  significant  (60  percent)  when  compared  with  the  40-percent 
increase  between  the  equator  and  the  pole  at  450  km.  This  predominance  in  the  global  component  was 
even  more  pronounced  at  200  km  and  appeared  to  be  in  substs  ^tial  agreement  with  the  satellite  drag 
data.  The  maximum  in  the  density  occurred  at  0200  u.L  at  high  latitudes  and  shifted  toward  0800  u.t.  at 
the  equator,  with  the  average  time  lag  of  6  h  deduced  from  satellite  drag  data. 

In  other  simulations  of  geomagnetic  substorms,  one  finds  that  the  thermosphere  at  high  latitudes 
is  significantly  perturbed  by  the  heat  and  momentum  deposited  there,  with  wind  speeds  approaching 
I  km/s  in  the  region  of  the  polar  cap  and  temperatures  elevated  some  500  K  in  the  auroral  zones.  At 
lower  latitudes  the  effects  are  far  less  dramatic,  but  are  nonetheless  still  significant  The  response  is 
asymmetric  with  respect  to  longitude,  reflecting  the  asymmetry  of  the  high-latitude  heat  sources.  Large- 
scale  disturbances  propagate  from  each  auroral  region  toward  equatorial  latitudes  with  a  time  delay  of 
about  2  h  at  330  km  altitude  and  about  5  h  at  125  km  altitude.  At  330  km  altitude,  the  equatorial 
thermosphere  is  influenced  mainly  by  the  interactions  of  two  converging  gravity  waves,  while  at  125  km 
it  is  influenced  mainly  by  the  mean  motion.  Equatorial  temperature  oscillations  of  magnitude  100  to 
200  K  occurred  at  330  km  altitude  driven  by  oscillations  in  the  mean  meridional  wind  of  about  200  m/s. 
Temperature  enhancements  at  125  km  were  much  smaller  (60  K)  as  were  the  meridional  wind  speed 
enhancements  (8  m/s).  The  zonally  averaged  zonal  wind  enhancements  were  found  to  be  confined 
primarily  to  the  high  latitudes. 

3.2.9. 1.1  Ionospheric  Storm  Effects.  During  and  following  a  geomagnetic  substorm,  the 
ionosphere  is  altered  as  a  consequence  of  several  competing  effects.  Some  of  these  will  lead  to  increases 
in  the  ionization  density  of  the  F2-tegion  (positive  storm  effects),  and  others  will  lead  to  decreases  in  the 
F2-peak  (negative  storm  effects).  These  are  briefly  summarized  for  the  subauroral  ionosphere: 

•  Positive  Storm  Effects: 

(1)  Traveling  aunospheric  disturbances  such  as  gravity  waves  will  cause  a  sudden  uplifting  of 
the  F-layer.  Through  collisions  the  neuual  disturbances  impart  to  the  ions  a  moU'^n  parallel  to  the 
geomagnetic  field  lines.  These  effects  propagate  equatorward  at  speeds  of  many  hundreds  of  meters 
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per  second.  The  duration  of  the  upward  drift  is  about  1  h  and  is  particularly  pronounced  in  the 
noon/aftemoon  local  time  sector. 

(2)  The  large-scale  wind  circulation  will  also  impart  a  motion  to  the  ions  that  is  parallel  to  the 
geomagnetic  field.  This  effect  may  last  a  whole  day  and  is  more  likely  to  occur  in  the  daytime. 

(3)  The  expansion  of  the  polar  ionization  enhancement  toward  lower  latitudes,  due,  for  example, 
to  particle  precipitation,  will  directly  increase  ionization  densiUes.  It  is  typically  a  nighttime  phenome¬ 
non. 


•  Negative  Storm  Effects; 

(1)  Changes  in  the  neutral  gas  composition,  and  particularly  an  increase  in  the  N2/O  ratio,  will 
lead  to  a  depletion  of  ionization.  This  is  most  clearly  observed  in  the  morning  sector,  with  an  anoma¬ 
lously  low  rate  of  ionization  after  sunrise.  The  effect  may  last  from  many  hours  to  days. 

(2)  The  displacement  and  expansion  of  the  ionospheric  trough  region  toward  lower  latitudes 
creates  a  steep  ionization  density  drop  and  is  typically  observed  in  the  aftemoon/evening  sector. 

Another  ionospheric  effect  of  geomagnetic  disturbances  is  the  propagation  of  disturbance  effects  known 
as  traveling  ionospheric  disturbances  (TlDs).  They  are  characterized  by  oscillations  in  the  ionization 
density  at  the  height  of  the  F-layer,  driven  by  gravity  wave-driven  oscillations  in  the  neutral  density  and 
neutral  winds. 

3.3  Solar  and  Geomagnetic  Indices 

Various  surrogate  indices  are  used  to  quantitatively  assess  the  levels  of  solar  activity.  One  of 
these  is  the  10.7  cm  (2800  MHz)  solar  radio  noise  flux,  designated  Fjo.?.  Although  it  is  the  EUV 
radiation  that  heats  the  thermosphere,  it  cannot  be  measured  at  the  ground.  The  F10.7  can  be  measured 
from  the  ground,  and  it  also  correlates  quite  well  with  the  EUV  radiation.  Although  there  are  instances 
when  the  correlation  is  not  good,  it  appears  unlikely  that  the  F10.7  radio  flux  will  be  replaced  by  another 
index  in  the  foreseeable  future. 

An  index  that  is  used  as  a  measure  of  episodic  type  solar  activity  is  the  planetary  geomagnetic 
activity  index  Op  (or  kp,  which  is  essentially  the  logarithm  of  Op).  It  is  based  on  magnetic  field  fluctuation 
data  reported  every  3  h  at  12  stations  between  geomagnetic  latitudes  48*  and  63*  and  selected  for  good 
longitudinal  coverage.  Although  it  is  the  high  latitude  ionospheric  current  fluctuations  that  drive  the 
magnetic  field  fluctuations  as  observed  at  these  stations,  it  is  not  the  magnetic  field  fluctuations  which 
are  driving  the  thermosphere.  Therefore,  the  correlations  between  observed  density  changes  and  the  Op 
index  are  not  always  good.  The  daily  planetary  geomagnetic  index,  Ap,  is  the  average  of  the  eight 
3-hourly  Op  values  for  that  particular  ^y. 

Table  3-1  lists  the  maximum,  mean,  and  minimum  13-month  smoothed  values  for  Fjoj  and  Ap 
throughout  a  mean  1  l-year  solar  cycle.  Figures  3-5  and  3-6  show  data  from  this  table.  The  F107  data  ait 
derived  from  sunspot  records  for  the  period  1749  to  1947  with  direct  F10.7  measurements  ther^ter. 
Table  3-1  covers  a  mean  duration  cycle  length  of  1 1  years.  The  standard  deviation  about  the  mean  length 
is  1.23  years  in  the  historical  record.  “Max”  and  “min”  are  the  historical  extremes  for  each  point  in  the 
cycle  and  have  been  determined  after  the  data  have  been  13-month  smoothed  and  constrained  to  the 
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mean  duration  cycle.  The  Ap  values  in  table  3-1  and  figure  3-6  are  derived  in  a  similar  fashion  based  on 
a  data  record  that  goes  back  to  1932. 

Table  3-1.  Maximum,  mean,  and  minimum  values  of  the  13-month  smoothed  10.7  cm  solar  radio 
noise  flux  and  geomagnetic  activity  index  over  the  mean  solar  cycle  (page  1  of  4). 
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Ap 

P, 
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1 

73.3 

69.6 

67.0 

11.5 

9.5 

7.6 

2 

nA 

69.7 

67.0 

11.7 

9.6 

7.7 

3 

14.0 

70.0 

67.0 

11.8 

9.7 

7.7 

4 

74.5 

70.4 

67.0 

11.9 

9.7 

7.6 

5 

74.9 
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67.0 

11.9 

9.7 

7.4 

6 

76.2 

71.1 

67.1 

12.2 

9.9 

7.3 

7 
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67.2 

12.5 
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7.2 

8 
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67.3 

12.9 
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7,3 

9 

81.5 
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67.4 

13.3 

10.6 

7.8 

10 

84.1 

73.6 

67.5 

14.1 

10.9 

8.1 

11 

87.7 

74.5 

67.7 

15.1 

11.2 

8,2 

12 

93.4 

75.7 

67.9 

15.7 

11.5 

8.3 

13 

97.9 

77.0 

68.0 

15.9 

11.8 

8.3 

14 

101.7 

78.4 

68.0 

16.4 

12.0 

8.3 

15 
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12.3 

8.5 

16 
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82.0 

68.0 

18.4 

12.7 

8.4 

17 
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84.0 

68.1 

18.7 

12.9 

8.5 

18 

129.1 

86.2 

68.4 

18.8 

13.1 

8.7 

19 

137.6 

88.5 

68.5 

18.6 

13.2 

9.0 

20 

143.4 

91.0 

68.6 

18.3 

13.2 

9.3 

21 

147.6 

93.7 

68.8 

18.1 

13.2 

9.7 

22 

151.7 

%.3 

68.7 

18.4 

13.4 

9.5 

23 

155.7 

98.9 

68.8 

18.4 

13J 

9.3 

24 

160.1 

101.6 

69.2 

17.6 

13.5 

9.1 

25 

164.8 

104.4 

69.7 

17.1 

13.6 

9.0 

26 

169.1 

107.2 

70.1 

17.4 

13.6 

9.1 

27 

173.0 

110.2 

70.6 

17.4 

13.6 

9.4 

28 

177.1 

113.2 

70.7 

18.5 

13.8 

9.8 

29 

186.1 

116.2 

71.3 

19.9 

14.0 

10.0 

30 

191.5 

119.3 

72.2 

19.9 

14.1 

10.0 

31 

194.3 

122.0 

72.6 

19.9 

14.1 

10.1 

32 

196.9 

124.3 

'3.3 

20.1 

14.1 

10.4 

33 

199.6 

126.5 

73.9 

20.4 

14.2 

10.2 

34 

204.2 

128.6 

74.1 

20.8 

14.2 

10.3 

35 
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131.0 
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14.1 
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36 

214.8 
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74.5 

21.0 

14.0 

10.6 

37 

211.2 

135.6 

74.6 

21.2 

14.0 

10.5 

38 

221.6 

137.6 

74.5 

21.6 

14.1 

10.4 

39 

226.9 

139.6 

74.1 

22.1 

14.1 

10.6 

40 

229.9 

141.4 

73.6 

22.2 

14.0 

10.8 

41 

231.7 

143.2 

73.5 

21.0 

13.7 

10.7 

42 

233.7 

144.6 

73.6 

20.1 

13.4 

10.4 

43 

235.6 

145.6 

74.0 

19.8 

13.3 

10.5 

44 

238.8 

146.7 

75.1 

19.3 

13.3 

10.7 

45 

242.8 

147.2 

75.8 

19.2 

13.3 

10.8 

46 

245.2 

147.7 

76.5 

19.0 

13.4 

11.0 

An 

244.5 

148.1 

78.1 

18.8 

13.3 

10.7 

48 

243.3 

148.4 

80.1 

18.6 

13.4 

10.8 

49 

244.7 

148.7 

82.5 

18.6 

13.4 

10.6 

50 

245.7 

148.2 

84.0 

18.3 

13.4 

10.2 

51 

243.3 

146.8 

85.5 

18.2 

13.5 

10.6 

.  52 

239.4 

145.7 

87.9 

18.7 

13.8 

11.3 

53 

235.0 

145.1 

89.5 

19.2 

14.1 

11.4 

54 

232.9 

144.9 

92.2 

19.6 

14.2 

11.3 

55 

233.3 

144.9 

93.8 

20.3 

14.4 

11.3 

56 

233.1 

144.7 

94.9 

21.0 

14.6 

11.5 

57 

231.2 

144.2 

95.0 

21.4 

14.8 

11.6 

58 

229.1 

143.5 

94.7 

21.2 

14.8 

11.6 

59 

228.1 

142.7 

94.9 

20.4 

14.7 

11.8 

60 

227.6 

142.3 

96.5 

20.7 

14.8 

12.1 

61 

226.7 

142.1 

97.3 

21.9 

15.1 

12.2 

62 

225.6 

141.3 

96.8 

22.7 

15.2 

12.0 

63 

223.0 

140.1 

96.0 

22.7 

15.1 

11.6 

64 

218.6 

138.4 

96.0 

22.3 

15.1 

11.2 

65 

215.3 

136.8 

96.6 

^  21.7 

15.1 

11.2 

66 

212.0 

135.5 

96.7 

21.5 

15.1 

11.2 

67 

206.9 

134.3 

95.1 

22.1 

15.1 

11.2 

68 

204.0 

133.0 

95.0 

23.1 

15.5 

11.3 

69 

203.6 

131.6 

96.3 

23.5 

15.6 

11.3 

70 

200.4 

129.8 

96.5 

23.4 

15.6 

11.2 

71 

196.8 

128.3 

94.7 

23.3 

15.7 

11.1 

72 

195.7 

127.3 

93.6 

23.1 

15.5 

10.8 

Table  3-1.  Maximum,  mean,  and  minimum  values  of  the  13-month  smoothed  10.7  cm  solar  laifio 
noise  flux  and  geomagnetic  activity  index  over  the  mean  solar  cycle  (page  4  of  4). 


Cycle 

^10.7 

Ap 

Pt 

Max 

Mean 

Min 

Max 

Mean 

Min 

110 

92.7 

81.1 

69.3 

16.5 

13.6 

8.8 

111 

92.0 

80.3 

69.0 

16.7 

13.5 

8.9 

112 

91.8 

79.6 

68.8 

16.9 

13.4 

9.0 

113 

91.4 

78.9 

68.5 

17.1 

13.3 

9.0 

114 

90.8 

78.2 

68.2 

17.4 

13.3 

9.0 

115 

90.1 

77.5 

68.2 

17.7 

13.1 

9.0 

116 

89.1 

76.9 

68.2 

17.6 

12.9 

9.2 

117 

88.2 

76.4 

68.2 

17.4 

12.7 

9.3 

118 

87.0 

75.9 

68.3 

16.9 

12.5 

9.2 

119 

85.4 

75.3 

68.3 

16.1 

12.2 

9.1 

120 

83.2 

74.8 

68.3 

14.7 

11.8 

9.1 

121 

80.5 

74.2 

68.3 

13.6 

11.5 

9.1 

122 

78.5 

73.5 

67.9 

13.7 

11.2 

8.9 

123 

77.6 

72.9 

67.6 

13.4 

10.9 

8.5 

124 

77.1 

72.3 

67.4 

13.0 

10.6 

8.1 

125 

76.9 

72.0 

67.4 

12.7 

10.5 

8.0 

126 

76.7 

71.6 

67.2 

12.4 

10.3 

8.0 

127 

76.5 

71.3 

67.1 

11.7 

10.1 

8.0 

128 

76.2 

70.9 

67.0 

11.2 

9.9 

8.0 

129 

75.2 

70.6 

67.0 

11.0 

9.8 

7.9 

130 

74.2 

70.3 

67.0 

10.9 

9.1 

7.2 

131 

74.0 

70.1 

67.0 

r  11.1 

9.2 

7.4 

132 

73.5 

69.9 

67.0 

11.4 

9.4 

7.6 

This  table  may  be  extended  by  repetition  of  the  1 1-year  cycle. 


3.4  Orbital  and  Suborbital  Neutral  Atmosphere  Model 

The  MSFC  Global  Reference  Atmosphere  Model  (GRAM)  provides  a  single  continuous 
representation  of  the  neut^  atmosphere  from  the  Earth’s  surface  to  2500  km  altitude.  GRAM  consists 
of  a  four-dimensional  (latitude,  longitude,  altitude,  and  time)  model  of  the  atmosfdiere  frmn  the  sur&ce 
to  25  km  that  is  faired  to  a  Barnett  model  of  the  atmosphere  flora  30  to  1 15  km,  which  is  then  faired  to  a 
modified  MSFOJTO  model  that  extends  from  90  to  2500  km.  There  is  a  description  of  this  model  in 
section  3.8.1  of  references  2,  G,  and  H.  The  computer  program  for  the  GRAM  and  the  appropriate  input 
parameters  for  the  model,  which  depend  upon  the  date  assumed  for  application,  are  available  upon 
request 


^  Note:  The  standard  deviation  about  the  mean  period  is  1.23  years  in  the  historical  record. 

^iire  3-5.  13-month  smoothed  values  of  F\o.7  over  the  mean  solar  cycle  (use  taMe  3-1 

for  numerical  values). 


Figure  3-6.  13-month  smoothed  values  of  geomagnetic  activity  index  (Ap)  over  the  nyan 
solar  cycle  (use  table  3- 1  for  numerical  values). 
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3.5  Marshall  Engineering  Thoinosphere  Model 

The  Marshall  En^neering  Thermosphere  (MET)  model  has  been  developed  to  represent,  in  so 
far  as  is  practical  for  engineering  applications,  the  variability  of  the  ambient  mass  density.  It  is  the  stan¬ 
dard  neutral  atmospheric  density  model  used  for  control  and  lifetime  studies  for  most  NASA  spacecoft 
projects.  A  description  of  the  model  can  be  found  in  references  8, 9,  and  10.  The  computer  prognm  for 
the  MET  orbital  atmosphere  model  is  available  upon  request. 

The  MET  model  is  an  empirical  model  whose  coefficients  were  obtained  from  Mtrfiitf*  drag 
analyses.  It  is  a  static  diffusion  model  and  is  essentially  the  Smithsonian’s  Jacchia  1970  model  with  two 
additions  from  the  Jacchia  197 1  model.  Inputs  to  the  model  are  time  (year,  month,  day,  hour,  and 
minute),  position  (altitude  and  geographic  latitude  and  longitude),  the  previous  day’s  solar  radio  flux 
(Fio.t),  the  centered  solar  radio  flux  averaged  over  six  solar  rotations  (Fio  jB),  and  the  Op  index  6  to  7 
hours  before  the  time  in  question  (for  some  studies  the  daily  planetary  geomagnetic  index,  A»,  may  be 
used  instead  of  the  3-hourly  Op  value). 

With  the  exospheric  temperature  specified,  the  temperature  can  be  calculated  for  any 
between  the  lower  boundary  (90  km)  and  the  upper  levei  (2500  km)  of  the  model  from  an  empirically 
determined  temperature  profile.  The  density  for  all  poin  ts  on  the  globe  at  90  km  altitude  is  assumed 
constant,  and  mixing  prevails  to  105  km.  Between  these ;  two  altitudes,  the  mean  molecular  mass  varies 
as  a  result  of  the  dissociation  of  molecular  to  atomic  oxygen.  At  120  km  altitude,  the  ratio  of  atomic  to 
molecular  oxygen  is  assumed  to  be  1.5.  Density  betw^n  90  and  105  km  is  calculated  by  iniegrarion  of 
the  barometric  equation.  For  altitudes  above  105  km,  fhe  diffusion  equation  for  each  of  the  individual 
species  (O2, 0,  N2,  He,  and  Ar)  is  integrated  upward  ^rom  the  105  km  level.  For  hydrogen,  the  imegn- 
tion  of  the  diffusion  equation  proceeds  upward  from  ;300  km  altitude.  The  total  mass  density  is  calcu¬ 
lated  by  summing  the  individual  specie  mass  densiti|*;s. 

The  total  density  is  then  further  modified  by  the  effects  of  the  seasonal-latitudinal  density  varia¬ 
tion  of  the  lower  thermosphere  below  170  km  altiU;>de  and  seasonal-latitudinal  variations  of  helimn 
above  500  km.  These  two  effects  have  been  incorporated  in  the  MET  model  using  equations  developed 
by  Jacchia  for  his  1971  thermospheric  model. 

The  final  output  of  the  MET  model  is  toL*  J  mass  density,  temperature,  pressure,  individual  specie 
number  densities,  mean  molecular  weight,  scale,  height,  specific  heats,  and  the  local  gravitational  accel¬ 
eration.  > 

i 

_  A 

The  total  mass  density,  die  temperature,  and  the  individual  species  all  have  the  same  phase 
variation  in  the  MET  model  (i.e.,  they  all  maximize  at  the  same  local  time).  For  some  studies  inwriving 
the  effects  of  various  s|:^ies  on  an  orbiting  soacecraft  where  accurate  phases  of  the  various  species 
within  an  orbit  are  required,  it  may  be  neces,^  to  use  the  MSIS  model.^^ 

The  values  given  in  table  3-2  define  the  Earth  pressure  parameters  for  space  vehicle  per¬ 
formance  analyses. 

3.5.1  Statistical  Analysis  Mode.  7he  statistical  analysis  mode  (SAM)  of  the  MET  model  was 
developed  for  use  during  those  time  perir/ds  when  actual  values  of  the  daily  and  162-day  values 

of  the  10.7  cm  solar  radio  flux  and  the  global  geomagnetic  index,  Op,  are  unknown  or  unavailable.  e.g.. 
when  the  time  of  application  to  a  ^ci^c  problem  is  sometime  in  the  future.  This  is  particularly  mie 
during  the  design,  development,  and  toting  phases  of  a  space  vehicle  when  there  are  only  rough 

/ 
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^timates  of  the  13*inonth  smoothed  mean  values  of  these  three  parameters  and  no  guesses  at  aD  about 

combined.  The  procedure  currently  used  in  spacecraft  development  wwk  is  to  assume 
that  the  dMy  and  the  162-day  mean  values  of  the  10.7  cm  solar  radio  flux  are  both  equal  to  the  IS-mooth 
^oothed  value  predicted  using  a  technique  developed  through  a  statistical  analysis  of  sonspot  data  fium 
1855  to  the  present  {see  ref.  E  for  sunspot  to  F10.7  conversion).  The  13-month  smoothed  values  of  a-  are 
predicted  in  the  same  manner.  This  procedure  provides  useful  estimates  of  the  long-term  density 
object  pnmanly  to  the  accuracy  of  the  predicted  13-month  smoothed  input  data.  For  periods  shorter  Ln 
90  days,  the  themospheric  properties  defined  by  the  use  of  the  13-month  smoothed  inputs  xe  no  lomsr 
repr^ntative  of  the  original  model  outputs  using  the  proper  input  parameters.  For  systems  which  are 
sensihve  to  thermospheric  effects  and  variations  over  time  periods  of  a  few  days  or  less  (e.t,  control  and 
poinung  system),  representative  samples  (ref.  D)  of  the  proper  model  inputs  taken  from  tfaeiiisaoiical 
dataset  should  be  used.  Alternatively,  for  some  applications  the  variations  may  be  treated  sradstically. 


Table  3-2.  Pressure  parameters  in  Earth  orbit  (ambient  pressure  in  Pascal;  Toir  in  parenthesis).* 


Minimum 

Nominal 

Maximum 

Altitude  (km) 

4.0E-5 

(3.0E-7) 

8.5E-5 

(6.4E-7) 

2.3E-4 

(1.7E-6) 

200 

1.8E-7 

(1.4E-9) 

1.5E-6 

(l.lE-8) 

1.7E-5 

(1.3E-7) 

400 

3.0E-8 

(2.3E-10) 

3.1E-7 

(2.3E-9) 

6.5E-6 

(4.9E-8) 

500 

1.7E-8 

(1.3E-10) 

8.3E-8 

(6.2E-10) 

2.7E-6 

(2.0E-8) 

600 

4.3E-9 

(3.2E-11) 

7.5E-9 

(5.6E-11) 

1.5E-7 

(l.lE-9) 

1000 

l.OE-11 

(7JE-14) 

Geosynchronous 

*Low  and  high  pressure  values  were  computed  from  the  MET  model  using  an  oifeit 
average  value  and  assuming  a  28.5*  orbit  inclinaUon.  The  following  inputs  were  med: 
Low:  date:  August  8, 0400  u.t.,  F,o.7  =  70.  A-  =  0  High:  date:  October  27 
1400  u.t.F,o.7  =  230.  =  400. 


Nominal  pressure  values  were  taken  from  USSA  1976  (ref.  A).  Geo-synchronous  pressure  valnes  were 
taken  from  reference  1.  To  convert  from  Pascal  to  Torr  multiply  by  7.5E-K  ^ 

K  ^  is  the  imtial  step  in  developing  the  statistics  of  these  shorter  period  paturintioas 

in  orbital  altitude  density  in  addition  to  providing  limits  on  the  magnitudes  of  the  variations  that  ae 
mSd  through  the  use  of  13-month  smoothed  solar  activity  input  parameters  to  the  original  MET 

f  .•  MET-SAM  also  provides  answers  to  another  frcquenUy  asked  question-what  is  the  percent 

recommended  density  value  wiU  be  exceeded  during  the  operational  phase  of  the 
vehicle.  Or,  how  confi^nt  are  you  that  the  recommended  density  value  will  not  be  exceeded  more  than 
5  percent  of  the  time .  The  answers  are  crucial  in  the  design  of  the  guidance  and  control  the 
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selection  of  the  altitude  at  which  the  vehicle  will  orbit  and  the  reboost  strategy  for  payloads  put  into 
LEO  by  the  shuttle.  The  MET-SAM  is  also  based  on  the  premise  that  most  applications  during  the 
development  phase  primarily  require  detailed  knowledge  about  the  maximum  and  minimum  densities 
that  will  be  encountered  with  limits  on  the  magnitude  of  the  variations  that  occur  during  monthly  time 
periods.  Research  is  in  progress  to  develop  a  technique  for  including  the  perturbations  within  an  orbit 

All  available  daily  values  of  the  10.7  cm  solar  flux  and  the  three  hourly  values  of  from  1947  i 
1991  were  used  to  calculate  values  of  the  global  minimum,  mean,  and  maximum  exospheric  tempera- 
Oires.  The  statistics  of  this  new  data  set  were  calculated,  and  the  results  were  indexed  to  five  different 
levels  of  solar  active  as  defined  by  the  values  of  the  13-month  smoothed  10.7  cm  solar  flux.  This  makes 
it  possible  to  statistically  include  the  3-hou  to  90-day  variations  in  these  three  temperatures  using  only 
known  or  estimated  values  of  the  13-month  smoothed  10.7  cm  solar  flux.  MET-SAM  is  not  yet  availabl 
for  general  release,  so  an  output  summary  table  of  expected  density  variations  is  included  in  table  C- 1 , 
appendix  C.  The  data  in  this  table  present  the  statistics  of  the  glob^  (spatial)  maximum  densities  for 
averaging  periods  as  short  as  3  hours,  the  temporal  resolution  of  the  Op  data  set 

Use  of  the  data  in  this  table  also  makes  it  possible  for  the  design  engineers  to  easily  accomplish 
any  required  trade  studies  while  managers  will  be  able  to  determine  just  what  risks  they  will  be  tal^g. 
For  ^plications  which  r'.quire  density  values  other  than  the  global  maximum  or  for  questions  concern¬ 
ing  the  application  of  the  data  in  engineering  analyses  contact  personnel  of  the  Electromagnetics  and 
Environments  Branch,  EL54,  MSFC,  AL  3S812.  Table  C-1  also  shows  the  median  value  of  the  global 
maximum  density  for  the  altitude  and  the  10.7  cm  solar  flux  bin  indicated,  computed  over  the  period  of 
record.  It  also  gives  the  differences  between  the  median  value  of  the  global  maximum  density  and 
several  other  percentile  levels  of  the  global  maximum  density.  “Percentile”  refers  to  the  fraction  of  time 
the  density  was  equal  to  or  less  than  the  indicated  value.  These  data  may  be  used  directly  to  indicate  the 
magnitude  of  the  variations  in  the  global  maximum  density  which  occur  over  time  periods  of  a  few  days 
If  appropriate  for  the  application,  an  additional  increment  may  be  added  to  account  for  model 
inaccuracies  and  variations  related  to  other  sources.  To  obtain  a  conservative  (upper  limit)  estimate,  we 
recommend  assuming  these  additional  variations  are  gaussian  with  a  standard  deviation  equal  to  0.2 
times  the  total  density. 

To  understand  the  frequency  of  occurrence  of  these  variations,  table  C-2  provides  the  probabili¬ 
ties  of  encountering  intervals  of  10  and  30  days  without  exceeding  the  indicated  percentile  level  (As 
with  table  C-1,  only  solar  and  geomagnetic  variations  are  considered.)  Probabilities  for  other  time  inter 
vals  can  be  estima'.ed  from  figure  C-1.  If  one  of  these  higher  percentile  levels  is  exceeded,  it  typically 
drops  back  down  within  a  short  time  period  (hours,  a  few  days  at  most).  The  distributions  are  highly 
deewed,  and  once  a  level  has  been  exceeded,  it  is  usually  exceeded  again  within  a  few  days.  On  the  othe; 
hand,  if  the  Sun  is  quiet,  it  may  remain  so  for  several  months. 
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IV.  PLASMA  ENVIRONMENT 


Slightly  below  the  temperature  transition  region  discussed  in  section  III  that  marics  the  lower 
boundary  of  the  thermosphere — and  associated  changes  in  composition  of  the  neutral  species — there 
is  an  important  transition  related  to  the  electromagnetic  properties  of  the  gas.  At  roughly  80  /cm  alti¬ 
tude,  there  is  a  division  between  the  lower  turbulent  neutral  gas  mixture  region  where  all  the  meteo 
rological  processes  occur  and  the  upper  region  where  solar  irradiation  produces  a  partially  ionized 
plasma  composed  of  O.  N2, 02,  He,  H,  0+,  He+,  NO,  O2*,  N2*,  and  electrons.  This  upper  region 
is  electrically  neutral,  with  the  most  abundant  neutral  being  O  and  the  most  abundant  ioo  being  O 
up  to  about  1000  km  altitude  where  and  He+  become  dominant.  A  plasma  is  a  qua^neutral  gas 
of  charged  and  neutral  particles  that  exhibits  collective  behavior.  The  particles’  movements  are  con¬ 
trolled  to  a  great  extent  by  the  Earth’s  magnetic  field  and  the  solar  wind,  but  their  collective  behavior 
and  movement  generate  electric  and  magnetic  fields  that,  in  turn,  affect  ^e  particle’s  motion  and  the 
motion  of  other  charged  particles  far  away. 

Plasmas  are  usually  described  by  their  density  (expressed  as  electron  number  density, 
ne  m~3),  the  chemical  composition  of  the  ions  (often  expressed  as  percentage  of  total  ions),  and  the 
electron  and  ion  temperatures  (expressed  in  Kelvin,  T,  or  as  energy  in  electron  volts  (eV)).  The  high 
energy  particles  which  also  may  be  present  (megaelectron  volts  (MeV)  range)  cannot  be  so 
described  and  they  interact  differently  with  a  vehicle  than  the  plasma  does.  They  are  discussed  in 
section  V  of  this  document 

The  plasma  environment  may  be  conceptually  divided  into  three  regions:  the  ionosphere, 
which  is  contained  within  the  magnetosphere;  the  magnetosphere;  and  the  solar  wind.  The  iono¬ 
sphere  is  characterized  by  its  low  temperature  and  high  density  relative  to  the  other  re^ns,  as  well 
as  its  predominantly  0+  composition.  Frequently,  this  region  is  considered  to  extend  to  about  KHX) 
km.  the  altitude  where  the  ion  density  begins  to  exceed  the  neuual  density.  Alternately,  an  arbitrary 
density  criterion  of  10®  m-^  places  the  “ionopause”  at  a  few  thousand  kilometers  altitude  in  the  polar 
regions  and  at  a  few  tens  of  thousands  of  kilometers  in  the  equatorial  regions.  The  equatorial  exten¬ 
sion  of  the  ionosphere  is  predominantly  composed  of  lighter  ions,  H'*’  and  He*,  and  is  usually  referred 
to  as  the  plasmasphere.  Above  the  ionopaus^plasmapause,  but  within  the  region  of  closed  magnetic 
field  lines,  the  plasma  is  still  dominated  by  photoionized  terrestrial  particles,  but  the  density 
decreases  and  the  temperature  increases  owing  to  energy  inputs  from  the  solar  wind. 

Even  as  low  as  about  60  km  altitude,  there  are  sufficient  long-lived  fiee-charged  particles  to 
affect  the  propagation  of  electromagnetic  waves,  to  support  various  modes  of  electrostatic  waves, 
and  to  support  detectable  (or  sometimes  significant)  currents  in  the  presence  of  electric  fields.  As 
altitude  increases,  the  charged  densities  increase  through  a  series  of  layers  to  something  on  the 
order  of  10*2  m-^  (on  the  day  side)  at  the  peak  density  altitude  (250  to  300  km).  With  additional 
increases  in  altitude,  they  decrease  to  the  order  of  10^  m-^  in  the  solar  wind  and  magneiospheiic 
plasma  sheet  or  to  even  lower  densities  which  are  found  in  the  magnetospherk  tail  lobes.  -lO®  ra-3 
(fig.  2-1).  At  low  altitudes,  the  temperatures  of  these  plasmas  are  typically  300  K  to  3500  K  (0.05  to 
0.3  eV)  except  in  the  polar  auroral  regions.  In  the  auroral  regions,  an  intense,  energetic  electron  flux 
often  precipitates  from  the  plasma  sheet  The  high  altitude  plasma  is  much  more  energetic;  typically 
105  K  (10  eV)  for  ions  and  1  to  5x105  K  (10  to  50  eV)  for  electrons  in  the  solar  wind,  -3x102  K  (3 
kiloelectron  volts  (keV))  for  electrons  and  -1x10*  K  (10  keV)  for  ions  in  the  magnetosphcric  plasma 
sheet  However,  the  dynamics  of  the  intervening  region  are  such  that  temperatures  can  sometimes 


4-2 


»'  “■*»  <^»«S«io  plasma  wiU,  ft=  p,„, 

plasma  is  repUced  by  d^wlm  wiSd'pj2iiw*SnTinS^  '**  "jeopause.”  where  the  terie; 
pause”  where  the  geomagnetic  field  is  reolaced  hv  thp  im  i  and  the  “magneii 

the  subsolar  point,  the  altitude  of  the  geooause  iJilllft  or  solar  wind  magneUc  fie! 

however,  in  the  extended  tail,  the  altitude^  ihp  oLn-.  ^  the  alutude  of  the  magnetopau> 
pause  owing  to  the  penetration  of  the  solar  wind  less  than  that  of  the  magneit 

The  magnetopause  ranges  from  6  to  10  Parth  raH-*  flanks  of  the  magnetosph 

radii  in  tile  antisunward  direction  hundreds  of  E^t 


4.1  Spacecraft-Plasma  Interactions 

space  vchklSVilSS^^'w^thA^^^^^  significant  differences  betwee 

plasma  at  very  high  orbits,  anilreraSSn  **'"  •°"«*Ph®t®.  the  hot  ten 

characteristic  of  higher  altitudes  penetrates  to  LEO 

charge  from  the  plasma  in  order  tc^establish  eleciSai  po  ®  accumulates  elcL 

spacecraft  charging.  Equmbrium  5^ulS  teat 

Pla.m.  prep^rtes  abd  U.e  spacecref.  design  and  opereU„rc^S«  fnteKfSe' 

ihe  iheimal  velocity  is  inve«l”rehw  to^^Srraas5*XT  '™t 
Spacecraft  surfaces  and  structure  is  ereatpr  ‘J®  electron  flux  to  passive 

surfaces  tend  to  accumulate  negative  charee  At  ^  *  result, 

and  of  low  energy,  so  equilibrium  is  established  w  thin  a^r  *"i  **'®  “  relatively  de 

potential.  Thus,  charging  of  piTvScefifus^atv  i  plasm 

active  surfaces,  e.g.,  ^lar  aiSys  L^SSre  tieS  i w?°  HoweveV,  for 

effects  can  occur,  depending  upijn  tl»  mundta^chemp  ‘^®™’  ^  '®***®‘* 

imposed  voltages.  ^  ^  heme  and  the  magnitude  of  the  spacecraft- 

energetic  electrons^precipirating^fr(S^  *e  maeSsohei^W^^^U^  ^^®®  "®8ative  potentials  by 
are  required  to  retard  this  flux  Md  aUow  equiUbrium  ?nTnpT^“  potentials 

region,  conditions  occur  in  the  wake  of  laree  structures  nr  ^  achieved.  Also  in  this 

whicle  is  involved,  where  the  low  energ7?la^St^^c“t^  naturally  so  that  the  enti 

balMcing  the  current  from  the  high  ene^VelS^n^S  ^annl?  ®  if^’  ineffective  in 

IS  similar  to  the  situation  in  GEO  wherc^the  nlSi^iS  ve^^pjf  ®*;?*^8ing  process  is  enhanced.  This 
current  densities  can  be  three  orders  of  magndeiSl^^if  n  ^  ^EO,  thermal 

from  surfaces  can  play  a  significant  role  in  haianpSn  ^  photoelectron  emissior 

become  imponamm  te S'^'*** 

-i.  end  .boded  Ws  *™'»P 
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4.2  Ionospheric  IMasma 

SHSSSSsi.7£r»sa- 

longitude.  The  plasma  density  drops  at  nilht  S?to  the  loVofw^L  pomagnetic 

free  ions  and  electrons  The  plasma  den^hv  fi»nHc  tn  s-  maintain  the  production  n\ 

^Tm^sKo  wind^altd  can 

SfsSHSSlSrT'T^ 

Wically  k.ctea^^r^jrtalLe. 

SpS^SS!  IT  "peak, 

SSToI^i^Ja^r"  ‘^""''' "  S  ySTaa. 
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affect  the  ionospheric  plasma.  The  boundary  region,  the  auroral  oval,  is  marked  by  the  beautiful  di;. 
play  of  auroras.  Surrounding  the  magnetic  poles,  the  oval  extends  to  near  75*  geomagnetic  latitude 
noon  and  65*  at  local  midnight  On  the  night  side,  the  ovai  is  well  marked  by  a  depletion  in  electron 
density,  the  so-called  trough.  On  the  dayside,  one  finds  a  region  of  enhanced  densities  just  inside 
the  oval,  the  so-caUed  raagnetospheric  cleft  The  electron  density  at  the  tip  of  the  cleft  is  almost  an 
order  of  magnitude  greater  than  it  is  at  the  bottom  of  the  trough.  During  magnetic  storms,  the  troug 
moves  equatorward  2*  per  unit  increase  in  Kp.  The  region  inside  the  oval  is  called  the  polar  cap. 

A  wide  variety  of  ionospheric  irregularities  have  been  observed,  predominantly  at  high  lati¬ 
tudes  and  during  the  equatorial  nighttime.  The  plasma  fluctuations  range  in  scale  from  hundreds  of 
kilometers  down  to  centimeters.  Plasma  instabilities  play  an  important  role  in  the  generation  of 
med'.um-scale  (kiloroeters)  and  small-scale  (meters)  uregularities.  Examples  of  irregularities  are 
patches  of  enhanced  ionization  in  the  E  region  (sporadic  E)  and  of  depleted  ionization  in  the  F  regii 
(spread  F).  Spread  F  is  most  frequently  observed  in  the  equatorial  nighttime  ionosphere.  The  irree 
larities  cause  signal  fluctuations  in  traversing  radio  waves,  known  as  scintillations. 

Influx  of  solm  plasma  into  the  tail  of  the  magnetosphere,  sometimes  preceded  by  solar  flare, 
can  cause  complex  ionospheric  disturbances  (storms).  The  most  consistent  pattern  is  an  enhance¬ 
ment  in  D  region  ionization.  These  effects  are  most  dramatic  at  auroral  latitudes,  but  significant 
modifications  in  the  ionosphere  occur  at  all  latitudes.  A  particularly  severe  event,  called  a  polar  cap 
absorption  (PCA)  event,  is  due  to  protons  arriving  from  the  Sun  and  causing  a  radio  communication 
blackout  over  a  considerable  time  period. 

4.2.2  Temperature-  The  main  source  of  energy  for  the  terrestrial  ionosphere  is  EUV  radiatio; 
from  the  Sun.  The  electrons  arc  heated  most  efflciently»  and  their  temperature  exceeds  the  tempera 
tures  of  the  ions  and  neutrals.  Electron  temperatures  increase  from  about  300  K  at  100  km  altitude  i 
about  3500  K  at  800  km  altitude.  Ion  temperatures  are  close  to  the  neutral  temperature  below  about 
400  km  altitude  and  increase  reward  the  electron  tempera  ture  above  that  altitude.  Below  150  km 
altitude,  the  high  neutral  densities  and  the  high  collision  frequencies  re.vj:lt  in  the  same  temperature 
for  electrons,  ions,  and  neutrals.  During  nighttime,  the  temperatures  of  all  species  are  close  togethe 

In  general,  plasma  temperature  are  lowest  at  the  geomagnetic  equator  and  increase  toward 
higher  latitudes,  due  to  the  uicredsed  influence  of  henting  by  precipitating  pmticles  nt  nuroral  Inti* 
tudes.  At  low  altitudes,  however,  the  electron  temperature  peaks  at  the  magnetic  equator,  reaches 
minimal  values  at  about  ±20*,  and  then  increases  toward  higher  latitudes.  This  behavior  is  the  mirro 
image  of  the  equatorial  anomaly  of  the  electron  density  and  illustrates  the  strong  anticorrelation 
between  electron  density  and  temperature. 

Roughly  spe^ng,  the  temperatures  increase  from  an  almost  constant  nighttime  value  to  an 
^ost  constant  daytime  value.  The  most  significant  deparU're  from  this  behavior  is  the  early  morn¬ 
ing  peak  in  electron  temperature.  It  is  most  pronounced  at  jie  magnetic  equator  at  about  3(X)  km 
altitude  (the  peak  temperature  exceeds  the  daytime  value  by  a  factor  of  2  to  3);  its  magnitude 
decre^  rapidly  toward  higher  and  lower  altitudes  and  toward  higher  latitudes.  The  temperature 
peak  is  a  result  of  the  sharp  increase  in  solar  heating  coupled  with  the  still  low  electron  densities 
from  the  preceding  night 

Tae  electron  temperature  is  almost  unaffected  by  the  solar  cycle,  in  contrast  to  the  increase  o 
almost  all  other  neutral  and  ionized  parameters.  This  again  is  a  result  of  the  close  coupling  with  the 
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electron  density  which  determines  both  energy  gain  and  loss  of  the  electron  gas.  The  simultaneous 
increase  of  both  terms  leaves  the  electron  temperature  nearly  unchanged. 

4.2.3  Dynamics.  The  Sun-induced  thermospheric  winds  provide  the  energy  source  needed  to 
drive  the  so-called  ionospheric  dynamo  which  maintains  the  system  of  ionospheric  currents  and 
fields.  On  the  sunlit  side  of  the  Earth,  two  large  vortices  of  electric  current  exist  in  the  quiet  equinox 
ionosphere,  the  current  flows  counterclockwise  in  the  northern  hemisphere  and  clockwise  in  the 
southern  hemisphere  (Sq  currents).  The  concentrated  current  at  the  magnetic  equator  represents  the 
equatorial  electrojeu  Magnetic  storms  severely  affect  thermt^heric  winds  and  ionospheric  currents. 
The  thermospheric  winds  and  ionospheric  drifts  are  of  the  order  of  100  m/s  and  can  reach  1000  m/s 
and  more  during  magnetic  storms.  Ionospheric  current  densities  are  of  the  order  of  10  jiA  m-2  and 
electric  fields  are  of  the  order  of  10  mV  m-‘.  The  Earth’s  magnetic  field  strength  is  typically  3x10*5 
Tesla  at  ionospheric  altitudes. 

The  solar  wind  blowing  past  the  Earth’s  magnetic  field  creates  a  magnetospheric  dynamo 
which  drives  plasma  across  the  polar  cap.  The  empirically  found  dependence  of  the  auroral  plasma 
corivection  on  solar  wind  parameters  illustrates  the  strong  coupling  between  the  solar  wind  and  the 
high  latitude  ionosphere. 

4.2.4  International  Reference  Ionosphere  90  flRWOt.  A  computer  code,  IRI90,  that  describes 
the  ionosphere  in  the  85  to  1000  km  altitude  range  for  geomagnetic  latitudes  up  to  about  60*  (the 
auroral  zones)  is  available  (ref.  C).  It  calculates  monthly  mean  electron  and  ion  species  (H+.  He+, 
0+,  NO+,  and  O2+)  densities  and  temperatures  as  a  function  of  latitude,  longitude  (geomagnetic  or 
geocentric),  time  of  day.  day  of  year,  altitude,  and  solar  F10.7  radio  flux.  (However,  the  model  has  not 
been  extended  to  cover  the  ionosphere  during  the  intense  solar  maximum  of  the  current  solar  cycle.) 
The  error  bars  on  the  model  during  quiet  times  are  a  factor  of  2  to  4  of  the  indicated  values  at  alti¬ 
tudes  below  the  F2  peak,  and  these  are  primarily  the  result  of  small  scale  variations  in  the  iono¬ 
sphere.  Above  the  F2  peak,  problems  with  the  scale  height  could  lead  to  error  bars  of  up  to  a  factor 
of  10. 


The  IRI90  was  used  to  generates  figures  4-1  and  4-2  which  show  global  density  and  tem¬ 
perature  contours  at  400  km  altitude  on  June  21,  approximately  the  summer  solstice.  Figures  4-3  and 
4-4  show  the  variation  of  density  and  temperature  with  altitude.  IRI90  is  an  updated  version  of 
earlier  models  (Bilitza,  1986)  which  adds  several  new  options  but  retains  the  general  character  and 
range  of  parameters  from  the  previous  editions.  Since  the  different  options  are  suitable  to  different 
applications,  care  should  be  exercised  in  their  selection. 


4.3  Auroral  Oval  Plasma 

As  mentioned  previously,  the  aurora  is  primarily  produced  by  high-energy  charged  particles 
precipitating  into  the  atmosphere  along  magnetic  field  lines.  One  result  of  these  fluxes  is  the 
incre^  of  local  plasma  density  by  factors  of  up  to  100  over  regions  of  tens  of  kilometers  in  latitudi¬ 
nal  dimension  and  hundreds  or  thousands  of  kilometers  in  longitudinal  dimension  in  the  auroral 
regions  (60*  to  70*  magnetic  latitude).  These  enhancements  occur  between  about  100  and  250  km 
altitude  (fig.  4-5).  Above  250  km,  the  thermal  plasma  may  be  depleted  above  intense  auroral  in  the 
midnight  sector,  faUing  far  below  1.0E>o  m-3.  In  this  region,  there  is  no  clear  distinction  between 
magnetospheric  and  ionospheric  phenomena.  There  is  a  very  wide  range  of  scales,  both  spatially  and 
temporally,  in  the  dynamics  present.  Small-scale  (meters  to  decameters)  irregularities  exist  and 
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move  with  the  prevailing  ambient  plasma  drifts.  Ion-ion  and  ion-neutral  collisions  in  the  lower 
thermosphere  tend  to  make  the  temperature  distribution  isotropic  while  anisotropies  still  remain  at 
higher  altitudes  (near  600  km).  Plasma  property  anisotropies  are  introduced  by  the  geomagnetic 
field. 


4-3.1  Global  Morphology.  The  conditions  in  the  magnetosphere  depend  on  whether  Bz.  the 
component  of  the  interplanetary  magnetic  field  (IMF)  normal  to  the  ecliptic  plane  with  positive 
values  northward,  is  directed  northward  or  southward.  There  are  three  distinct  magnetospheiic 
states:  (1)  the  quiet  state,  (2)  the  active  polar  cap  stale,  and  (3)  the  active  auroral  oval  state 

4.3. 1.1  The  Quiet  Magnetosphere.  In  this  state,  the  energy  input  from  the  solar  wind  is 

minimized  by  a  configuration  that  minimizes  magnetic  recombination  on  the  dayside.  A  weak  but 
northward  IMF  Bz  is  typical.  The  auroral  oval  has  the  following  characteristics: 

(a)  Visually  small  in  diameter  and  circular 

(b)  Somewhat  expanded  near  the  dawn  and  dusk  regions 

(c)  Relatively  thin  near  midnight. 

4.3. 1.2  The  Active  Polar  Cap  State.  At  this  time,  characteristics  of  this  state  are  still  contio- 
versitd.  and  more  detailed  observational  data  are  required  before  a  consensus  can  be  reached.  How¬ 
ever,  it  is  believed  that  it  occurs  during  strongly  northward  IMF  Bz  conditions  and  is  characterized  by 
the  presence  of: 

(a)  Sun-aligned  arcs 

(b)  Additional  magnetic-field-aligned  current  systems 

(c)  Distortions  or  disruptions  in  the  traditional  two-cell  polar  convection  patterns. 

4.3.1.3  The  Active  Amoral  Oval  State.  The  southward  IMF  Bz  plays  a  prominent  role  in  this 
state.  Auroral  boundaries  are  active,  with  rapid  poleward  movement  of  the  poleward  boundary  and 
less  rapid  equatorward  movement  of  the  equatorward  boundary  during  the  expansion  phase  of 
magnetic  substorms.  Although  there  is  some  controversy,  it  iqipears  as  if  the  optical  signature  of 
this  poleward  movement  ends  before  the  ground  signature  of  the  electrojet  activity  ceases.  The 
equatorward  particle  precipitation  boundaries  move  equatorward  with  increasing  geomagnetic 
activity.  The  ion  boundaries  are  equatorward  of  the  electron  boundaries  in  the  disk/evening  sector, 
while  the  electron  boundaries  are  ahead  in  the  morning  sector.  The  electron  boundaries  are 
equatorward  of  the  ion  boundaries  at  all  local  times  except  the  disk  when  the  geomagnetic  activity  is 
very  low — very  quiet  times. 

Auroral  features  are  also  very  active  in  this  state.  During  the  expansion  phase  of  an  auroral 
substorra,  the  poleward  and  equatorward  motions  of  the  nighttime  auroral  can  begin  at  localized 
sites  in  local  time,  primarily  between  2200  and  2400  h  local  time.  Sometimes  an  “eye-shaped" 
structure  with  substantial  north-south  structures  embedded  can  form  as  the  expansitm  proceeds  in 
the  ept-west  as  weU  as  north-south  direcuons.  The  configuration  of  the  east-west  expansion  is  a 
function  of  the  condition  of  the  IMF.  At  the  most  extreme  portion  of  the  poleward  expansion,  the  arcs 
can  intensity  rapidly  and  form  vortices  along  the  arc’s  length.  Westward  traveling  surges  (WTSs) 
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have  been  observed  to  propagate  nearly  7000  km.  Sometimes  WTS  activity  is  die  predominant 
auroral  response  to  substorm  activity. 

4.3. 1.4  Electrodynamics  and  Convection.  There  are  many  statistical  or  empirical,  as  well  as 
several  theoretical,  models  of  the  convection  and/or  electric  fields  in  the  auroral  zones.  Specific  use 
^11  dictate  which  model  should  be  used  in  engineering  analyses.  Most  of  these  models  stress  the 
importance  of  the  IMF  By,  the  dawn-dark  component  of  the  IMF,  in  controlling  the  dayside  convec¬ 
tive  patterns.  Results  of  analyses  have  shown  that  substantial  structure  can  be  present  within  the 
distribution  of  field-aligned  currents.  Results  of  studies  have  also  shown  that 

(a)  In  the  morning  sector,  the  westward  electrojet  is  dominated  on  the  poleward  side  by 
strong  electric  fields  (with  low  auroral  luminosity  and  conductivity)  and  on  the  equatorward  side  by 
large  conductivity  values.  The  equatorward  side  is  also  the  region  of  precipitatmg  electrons  (upward 
electrical  currents). 

(b)  On  the  evening  side  in  the  eastward  elecu-ojet,  the  reverse  trend  is  observed. 

(c)  The  westward  electrojet  is  also  centered  in  the  region  of  highest  auroral  luminosity,  while 
the  eastward  electrojet  is  in  the  region  of  the  lowest  luminosity. 

(d)  The  brightest  auroral  regions  are  associated  with  upward  electrical  currents,  both  in  the 
evening  and  morning  sectors. 

(e)  The  polar  cap  potential  is  modified  rapidly  by  changing  IMF  conditiwusi  however,  the 
nightside  auroral  zone  response  is  delayed  substantially. 

(f)  Conductivity  gradients  play  crucial  roles  in  the  formation  of  auroral  surges. 

Once  again,  understanding  of  the  characteristics  of  the  global  aurora  is  progressing  rapidly. 
However,  there  are  many  differences  between  models  and  observational  data,  and  the  use  of  any 
empirical/statistical/theoretical  model  will  be  highly  dependent  upon  the  application.  Care  should  be 
utilized  in  the  selection  of  the  “proper”  model  and  experts  in  the  field  should  be  consulted. 

4.3.2  Aurora  Morphology.  Aurora  occur  in  all  three  altitude  regimes,  D,  E.  and  F.  with  the  E 
layer  dominant  in  terms  of  total  auroral  precipitation  energy.  Electron  density  profiles  in  the  auroral 
region  are  dependent  upon  the  energy  distribution  of  the  incident  precipitating  pvticles.  while  the 
auroral  emissions  are  dependent  upon  the  photochemistry  of  the  auroral  ionosidiere  as  well  as  the 
energy  distribution  of  the  precipitating  particles.  The  spectrum  of  the  precipitating  particles  deter¬ 
mines  the  altitude  at  which  the  particle  energy  is  deposited  and,  therefore,  the  i^plicable  photo¬ 
chemistry.  since  neutral  composition  and  density  depend  strongly  on  altitude.  EiAanced  electron 
densities  are  produced  by  impact  ionization  due  to  precipitating  electron  fluxes.  The  aurccal  E  region 
extends  about  2*  farther  equatorward  than  the  F  region,  and  this  extension  is  formed  by  proton 
precipitation,  in  contrast  to  the  mainly  electron  precipitation  that  forms  the  F  region.  There  is  a 
strong  seasonal  and  universal  time  (UT)  control  of  the  auroral  ionosphere. 

The  aurora  occurs  at  all  local  times  in  the  auroral  ovals  that  are  rings  at  idmut  65*  magnetic 
latitude  around  each  magnetic  pole.  The  most  intense  and  energetic  fluxes  occur  near  local  midntghL 
The  location  of  the  aurora  and  its  intensity  va^  with  solar  (and  resulting  magnetaspberk)  activity. 
Figures  4-6  and  4-7  illustrate  the  spatial  distribution.  Global  plots  of  the  average  integral  energy 
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flux  and  the  average  energy  of  precipitating  electrons  are  presented  in  polar  spectrogram  format  in  a 
magnetic  local  time-corrected  geomagnetic  latitude  coordinate  system  for  each  of  four  levels  of  Kp. 
Plots  apply  to  both  poles.  A  statistical  analysis  of  the  probability  of  encountering  various  levels  of 
aurora  particle  flux  may  be  found  in  McNeil  et  al.,‘3  and  a  detailed  discussion  of  auroral  morphology 
may  be  found  in  Jursa  (ref.  F). 

In  the  dayside  auroral  zone,  magnetic  field  lines  extend  to  the  dayside  magnetopause  and 
provide  fairly  direct  access  for  the  solar  plasma  to  the  topside  ionosphere.  The  solar  plasma  strongly 
heats  and  drives  winds  in  the  day  side  auroral  ionosphere,  generally  in  the  direction  away  from  the 
Sun  and  across  the  polar  cap  toward  the  nightside.  The  strong  ionospheric  heating  associated  with 
this  process  generates  the  largest  outflows  of  ionospheric  plasma  anywhere  at  the  Earth  and  forms  a 
global  scale  “ion  fountain”  aaoss  the  polar  cap.  Global  electric  systems  flow  along  magnetic  field 
lines  into  and  out  of  the  auroral  ionosphere.  Where  these  are  most  intense  and  in  the  upward  (tirec- 
tion,  they  exceed  the  capacity  of  the  plasma  to  carry  them  with  thermal  motions,  and  discharges 
occur  in  which  electrons  are  accelerated  downward  and  ions  are  accelerated  upward.  These  dis¬ 
charges  greatly  enhance  the  brightness  of  the  aurora  where  the  electrons  are  incident  and  produce 
outward  flowing  ion  beams  which  are  far  out  of  local  thermodynamic  etprilibrium  and  unstable,  gen¬ 
erating  natural  plasma  wave  noise  at  a  variety  of  frequencies. 


4.4  Polar  Cap  Plasma 

In  the  polar  cap,  the  differences  between  magnetosphere  configurations  are  more  pronounced 
than  in  the  auroral  oval  region.  The  configurations  depend  on  whether  the  IMF  component  nwmal 
to  the  ecliptic  plane  with  positive  values  northward,  is  directed  northward  or  southward.  Current 
understanding  of  the  polar  cap  convection  under  southward  IMF  conditions  is  well  developed  eveti 
though  the  convection  is  complex  and  multiple  flow  entry  regions  exist.  However,  observed  traveling 
convection  vortices  in  the  polar  cleft  region  are  inconsistent  with  present-day  current-closed  models. 

4.4, 1  TnnnCTrfiftre  and  Polar  Wind.  Elensity  Structures  in  both  the  E  and  F  regions  show  that 
in  darimcM  the  ratio  of  enhanced  density  to  the  background  density  remains  constant  for  many  tens 
of  hours  and  that  the  vortices  only  disappeared  when  they  were  convected  into  regions  of  sunlight  or 
auroral  precipitation.  There  are  systematic  differences  between  the  winter  electron  density  signa¬ 
tures  in  the  topside  ionosphere  in  the  southern  and  northern  polar  regioas  with  the  southern  hemi¬ 
sphere  having  the  lower  densities.  The  region  above  the  polar  ionosphere  contains  mainly  iono¬ 
spheric  plasma  flowing  generally  away  from  the  Sun  across  the  polar  cap  and  upward  into  the  mag¬ 
netosphere.  Only  weak  fluxes  of  energetic  components  of  the  solar  wind  plasma  enter  the  polar  cap 
region,  sometimes  referred  to  as  “polar  rain.”  The  light  ion  component  of  the  ionospheric  plasma, 
which  flows  upward  into  the  magnetosphere  even  without  solar  wind  energy  inputs,  is  referred  to  as 
the  polar  wind.  The  heavy  ion  component  re^onds  to  energy  inputs  from  the  solar  wind,  to  heat- 
driven  thermospheric  dynamics,  the  solar  cycle  in  UV  inputs  to  the  theimoqrhere,  and  inagnetic 
activity.  This  combined  light  and  heavy  ion  source  of  plasma  to  the  magnetosfdiere  provides  much  of 
the  material  from  ndiich  the  energetic  particle  populations  are  generated.  Theoretical  research  on  this 
area  has  shown  that  results  are  highly  model  dependent  and  that  two-*nensional  models  produce 
different  results  than  three-dimensional  models.  Caution  in  use  of  modds  is  strongly  advised. 

The  polar  cap  region  becomes  active  with  Sun-aligned  auroral  ares  when  the  IMF  is  north¬ 
ward.  Field-aligned  clectrodynamic  systems  which  close  in  the  ionosphere  are  associated  with 
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these  arcs.  Plasma  drift  inside  the  arcs  is  primarily  antisunward.  while  outside  the  arcs  both  san- 
ward  and  antisunward  drifts  occur. 

Models  of  the  ionosphere  have  not  progressed  to  the  poiitt  where  they  can  be  used  for  reli¬ 
able  forecasts  of  future  conditions  especially  those  related  to  magnetic  storm  and  substorm  rfyam- 
ics.  Major  questions  of  convection  patterns  are  still  unanswered.  Even  the  static  conditions  yt- 
dicted  by  the  IRI90  model  frequently  do  not  match  observational  data.  Great  care  must  be  exeicised 
in  the  use  of  all  ionospheric  models. 


4.5  Geosynchronous  Altitude  Plasma 

The  geosynchronous  altitude  plasma  environment  is  very  complex  and  dynamic.  The  flmes  in 
GEO  can  be  quite  energetic  and  are  highly  variable  with  magnetic  activity  e^cially  during  geom^- 
netic  substorms.  The  values  given  in  table  4-1  are  an  estimate  of  the  90th  percentile  worst  changing 
case  environment  assuming  a  single  Maxwellian  representation  of  the  environment  Figure  4-S  pro¬ 
vides  a  recommended  time  variation  sequence  suitable  for  modeling  the  worst  effects  of  a  geomag¬ 
netic  storm.** 

Table  4-1.  Worst-case  plasma  environment  in  geosynchronous  Earth  orbit 


Electron  number  density,  rtg  m~^ 

1.12  E* 

Electron  temperature.  T^,  eV 

1.2  E* 

Ion  number  density,  n,-  m-^ 

2.36  E5 

Ion  temperature,  7,-,  eV 

2.95  E-* 

4.6  Transient  Fluxes  in  Low  Earth  Polar  or  Geosynchronous  Orbit 

In  polar  LEO.  the  important  transient  and  energetic  fluxes  occur  in  the  auroral  zone.  Figne 
4-9  shows  typical  aurora  fluxes  from  rocket  measurements  at  about  320  km,  and  figure  4-10  shows 
extreme  values  for  space  vehicle  design  derived  from  satellite  measurements  made  at  840  km.  Dis¬ 
tribution  functions  of  electrons  (left)  and  ions  (right)  for  very  intense  auroral  fluxes  are  shown  in 
figure  4-10.  These  particles  are  not  very  penetrating  but  may  be  significant  fot  charging  at 
above  250  km. 
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Figure  4-1.  Plasma  density  (m-3)  at  400  km  for  June  21  and  solar  minimum  conditions  (F la? 
0:0:0  G.m.L,  results  from  IRI90  with  default  options  selected). 
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Figure  4-2.  Electron  energies  (eV)  at  400  km  for  June  21  and  solar  minimum  coodiffims 
(Fio.7  =  70, 0:0:0  G.m.L,  IRI90  results  with  default  options  selected),  fi 
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longitude  plane,  IRI90  results  with  default  options  selected). 
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longitude  plane,  IRI90  results  with  default  options  selected). 
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Figure  4-5.  Ionospheric  electron  density  in  an  aurora  (xlO^ 


Figure  4-6.  Global  plots  of  the  average  integral  energy  flux  (keV/cm^  s  sr)  of  precipitating  electrons  are  presented  in 
polar  spectrogram  format,  in  a  magnetic  local  time  corrected  geomagnetic  latitude  coordinate  system, 
for  four  levels  of  Kp.  Plots  apply  globally  to  both  poles  (page  1  of  2).^^ 


Figure  4-6.  Global  plots  of  the  average  integral  energy  flux  (keV/cm^  s  sr)  of  precipitating  electrons  are  presented  in 
polar  specU’ogram  format,  in  a  magnetic  local  time  corrected  geomagnetic  latitude  coordinate  system, 
for  four  levels  of  Kp.  Plots  apply  globally  to  both  poles  (page  2  of  2). 


re  4-7.  Global  plots  of  the  average  energy  (keV)  of  precipitating  electrons  are  presented  in  polar  spectrogram  format, 
in  a  magnetic  local  time  corrected  geomagnetic  latitude  coordinate  system,  for  four  levels  of  Kp. 

Plots  apply  globally  to  both  poles  (page  1  of  2).^^ 


Figure  4-7.  Global  plots  of  the  average  energy  (keV)  of  precipitating  electrons  are  presented  in  polar  spectrogram  format, 
in  a  magnetic  local  time  corrected  geomagnetic  latitude  coordinate  system,  for  four  levels  of  Kp. 

Plots  apply  globally  to  both  poles  (page  2  of  2).i^ 
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V.  PENETRATING  CHARGED  PARTICLES 

.*  **^"^*"^**"8  particles,  often  referred  to  as  charged-panicle  radiation  or  corpuscular 

adiauon,  presents  a  significant  challenge  to  the  design  and  operation  of  a  spacecraft  This  is  because 

mf  reriak  Vh spacecraft  A  high  level  of  radiation  will  significant  affect 
materials,  chemical  proce^.  and  living  organisms,  especiaUy  the  crew,  ft  will  also  affect  electronics  bv 
( 1 )  causing  soft  upsets  (referred  to  as  single  event  upsets  (SEUs)),  (2)  degrading  performance,  and  (3)  ^ 
producmg  permanent  damage.  In  addition,  ionizing  radiation  will  affect  Se  proSoZfS^toSish 
Optical  matenals  by  altering  their  optical  properties.  The  quantitative  lelaUonship  between  tte^ount^of 
lomzanon  and  ihe  effects  depends  upon  the  tSnre  of  Jnfiadon  in  a  c^aStd^ 

It  is  convenient  to  divide  the  charged  particles  into  two  groups:  magnetospheric  particles  and 
cosmic  rays.  Magnetospheric  particles  arc  accelerated  from  thermal,  low-cSrgy  plasma  by  processes 
inside  the  magnetosphere  and  occur  only  within  terrestrial  space.  Cosmic  rays  exist  in  interplanetarv 

Wnd!  ?"‘®''.‘®*^strial  space  from  outside.  Within  terrestriai  space,  the  moS  of  Si 

kinds  of  charged  particles  is  controlled  by  the  geomagnetic  field. 

5.1  Magnetospheric  Partides 

Trapped  Radiaiion-  These  particles  are  trapped  in  the  geomagnetic  field  and  perform  manv 

ticlw^eTurawrSi?iSl**fi*^^  and  around  the  Earth  during  their  lifetime.  (Other  terms  £  these  par-^ 
^  ^  trapped,  geomagneacaUy  trapped,  or  Van  Allen  radiation.)  The  range  of  energies  is 
ather  l^ge  and  is  centered  in  tiie  tens  of  ke  V  for  electrons  and  MeV  for  protons  The  fluxes  aip  oenor 

P  magnetic  storms,  and  the  isoflux  contours  have  axisymmeu^ 

mound  UieEar^.  (More  prca.sely.  the  flux  in  three-dimensional  space  may  be  mappedta^SSon  of 
nvo  vanablw:  field  strength R  and  £  or B  and  a. TTie £-shell  value  is  definfdkSn  vS oS 
ocumenL  l^is  means  that  the  flux  contours  in  a  meridian  plane  can  represent  the  entire  three-dimen 
sional  morphology.)  Both  protons  and  electrons  occur  (see  reference  1^  a  general  discussion)! 

fhP  radiation  is  composed  of  both  protons  and  electrons.  These  particles  may  originate  on 

the  Sun,  m  the  atoosphere,  or  may  be  produced  by  cosmic  ray  albedo  neutron  decay.  DecaviMneutrons 
projected  upward  from  the  atmosphere  are  the  major  source  for  trapped  protons  Once  thev^w  hrjin 
£t”th  they  are  guided  and  transported  by  geomagnetic  and  electric  fields.  In  the  absena  of  fluctua 
uons,  the  ^>^cles  retmn  their  energy  and  move  in  stable  orbits  on  a  constant  L  shell.  Particles  leave  the 

radiation  belts  when  Aey  collide  with  the  neutral  atmosphere.  This  is  more  likely  to  occur  the  deeper  an 
orbit  penetrates  into  the  neutral  atmosphere.  'c  u»ciy  lo  occur  me  deeper  an 

The  general  shape  of  the  Van  Allen  belts  foUows  the  shape  of  the  geomagnetic  fieW  twc  m»onc 
vehicle  penetrates  most  dwply  into  the  belts  in  the  i^ion  of  L  South  AU^&iwT^ 
Because  the  flux  is  incrcssuig  with  altitude  in  the  region  of  300  to  looo  Wn  ttisa 

r  ^  ~  V-  figureWe  SLlSo^f’S^e’cSo'^^^^^^ 

raendional  pl^e  are  shown  m  figures  5-2  and  5-3,  and  the  profile  of  proton  flux  with  Stimdc  S? 
minimum  IS  shown  in  figure  5^.  TTie  very  steep  inner  gradient  is  conftoUed  by  thTex^?StiS,i^^n- 

nf-  "h  atmosphenc  density  which  varies  with  solar  activity.  The  atmosphere  is  more  ^ 

extended  (higher  density  at  a  given  altitude)  when  the  Sun  is  active  Thus  ar  . 

flux  is^eater  when  the  Sun  is  quiet  The  profile  atthe  inJ^r  eir^JSS^Sy  Itited  b^  ' 
^o^here  modulated  by  solar  activity.  These  changes  do  not  show  cleariy  on  the  figures  listed  above 
However,  they  can  be  shown  by  a  contour  plot  against  B  and  £  as  shown  in  figure  5-5. 


5-5 


The  trapped  flux  environment  at  solar  maximum  and  minimum  is  defined  the  computer  codes 
AP8MAX  and  AP8MIN  for  protons  and  AE8MAX  and  AE8MIN  for  electrons.^’  »  These  compmer 
codes  are  available  upon  request  It  is  essential  that  the  magnetic  field  used  with  the  proton  modds  be 
the  model  of  1%5  (Epoch,  1964),  i.c.,  IGRF  1965  SO-term  projection  to  1964  model  using  the 
coefficients  in  table  5- 1,  for  solar  minimum  and  1970  (Epoch,  1970),  i.e.,  the  Hurwitz  US  C&GS  1970 
field  and  the  coefficients  in  table  5-2,  for  solar  maximum.^^  Use  of  current  (1991  model)  field  values,  or 
those  projected  for  the  1990’s,  will  give  flux  values  far  higher  than  real  at  projected  operating  altitudes. 
This  is  because  the  Earth’s  field  is  decaying,  which  lowers  the  mirror  point  resulting  in  loss  of 

particles  to  the  atmosphere  at  these  low  altitudes.  Since  the  magnetic  moment  value  of 0.31 165300  ws 
only  appropriate  near  1960,  the  ALLMAG  routines^*  for  computing  B  and  L  and  the  ORP  program  for 
looking  up  fluxes  in  the  AP  and  AE  models  must  be  modified  by  replacing  the  constant  magnotir 
moment  value  with  the  following: 


A/  =  (C2j2+G222+C22i)1/2/100000  (5-1) 

ivhere  are  the  field  expansion  coefficients:  The  revised  computer  code  is  available  upon  request 

In  calculating  total  dose  in  LEO,  it  will  be  found  that  txapped  protons  contribute  nearly  the  endue 
unount  with  three  exceptions:  at  the  lowest  althudes  (below  about  300  km),  the  contribudon  fiom 
rapped  particles  becomes  so  small  that  galactic  cosmic  rays  (GCRs)  make  the  largest  contribution.  For 
'ery  thin  shields  (<0.3  gfcm2),  trapped  electrons  are  more  important  than  trapped  protons.  At  hi^ 
nclination  orbits,  solar  flare  event  pardcles  make  significant  contribudons.  To  repeat,  because  of  the 
apid  variation  with  altitude,  it  is  essendal  dial  the  correct  numerical  models  be  used  with  the  conect 
nput  parameters.  It  should  be  borne  in  mind  that  these  models  are  correct  only  to  withia  a  factor  of  3  fin- 
rapped  proton  fluxes  and  within  a  factor  of  5  for  tra(^d  electrons.  Vaiiadons  can  occur  during  very 
arge  magnetic  storms  that  last  a  few  days. 


Table  5-1.  Spherical  harmonic  coefficients  of  the  International  Geomagnetic  Reference 

Reid  (KjRF)  1965  (page  1  of  2). 
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-118 
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5. 1.2  Nonisotropic  Effects.  The  proton  flux  is  nonisotropic  in  the  South  Atlantic  Anonaaly 
region,  a  factor  which  will  affect  the  results  of  calculations  involving  complex  shielding  geomeiiy.  Hie 
flux  is  nonisotiopic  because  the  protons  follow  a  helical  path  about  a  nu^tic  field  line.  As  the  Held 
intensity  increases,  both  the  diameter  and  the  pitch  of  the  helix  decrease  until  the  pitch  becomes  zera 
The  point  with  zero  pitch  angle  is  called  the  mirror  point,  and  the  center  of  the  helical  p»th  is  ratipH  the 
guiding  center.  From  here  the  helix  reverses  direction  a:id  protons  travel  up  the  field  line  toward 
decreasing  field  intensity  and  away  from  the  Earth.  In  the  South  Atlantic  Anomaly,  almost  all  the 
protons  observed  are  near  their  mirror  points.  Thus,  the  flux  is  anisotropic,  with  most  of  the  flux  arriving 
from  a  narrow  band  perpendicular  to  the  local  geomagnetic  field  direction.  Any  protons  there  which  are 
not  near  the  mirror  point  will  travel  deep  into  the  atmosphere  and  be  scauered  or  stopped  by  atmo^ihcTic 
interactions. 

Atmospheric  interactions  also  affect  the  proton  angular  distribution  in  another  fashion.  Protons 
that  are  observed  traveling  eastward  are  following  guiding  centers  above  the  observation  point,  and 
protons  traveling  westward  are  following  guiding  centers  below  the  obsovation  point  The  gyroradius 
for  energetic  protons  in  the  anomaly  is  on  the  same  order  as  the  atmospheric  density  sctiic  height  Thus, 
westward  traveling  protons  encounter  a  significantly  denser  atmosphere  and  are  mote  likely  to  suffer 
atmospheric  inieractions  and  be  lost  The  resulting  energy-dependent  anhotropy  is  called  the  east-west 
effect  For  ty{»cal  values  of  the  magnetic  field  parameters  and  atmoroheiic  scale  height  in  the  anomaty. 
the  ratio  of  eastward  traveling  to  westward  traveling  flux  is  3  for  50-MeV  protons  and  22  for  40O-MeV 
protons.  Ev  iluatioii  of  the  anisotropic  effects  is  complex,  but  a  computer  code  for  use  in  rnnjuHTjnp 
with  the  AP8MAX  and  AP8MIN  models  is  available  upon  request 

5.2  Cosmic  Rsqrs 

For  the  purpose  of  this  locument  the  term  “cosmic  rays"  applies  to  electrons,  protons,  and  the 
nuclei  of  all  elements.  Outside  the  Earth’s  magnetosphere,  ail  Giactic  Cosmic  Rays  (GCRs)  may  be 
treated  as  isotropic  (to  within  about  10  percent).  Within  the  Earth’s  magnetic  field,  the  less  energetic 
nuclei  are  deflected  by  the  field,  and  at  lower  altitudes,  the  flux  becomes  anisotropic  while  the  energy 
spectrum,  the  total  flux,  and  the  directional  characteristics  are  modified  by  geomagnetic  <hu»Ming  For 
every  space  vdiicle  mission  phase,  this  effect  must  be  evaluated  and  int^ted  over  many  oibits. 

The  source  of  cosmic  rays  is  either  galactic  or  solar.  GCRs  originate  outside  the  solar  system  and 
are  thought  to  permeate  our  galaxy.  Since  the  condition  of  the  interplanetary  medium  tluough  which 
cosmic  rays  propagate  depends  on  solar  activity,  the  number  reaching  the  Earth's  orbit  is  modulated  by 
the  1 1-year  solar  cycle.  During  the  maximum  activity  phase,  the  intensity  near  the  Earth  rfeoyaistK.  and 
during  the  minimum  activity  phase,  it  increases. 

Solar  panicle  events,  in  contrast,  originate  in  the  Sun  and  are  produced  in  solar  flares.  They  are 
lower  in  energy  than  GCRs  (1  MeV  to  1  GeV/nucleon)  and  ate  mostly  pnions  and  alpha  particles^  Th^ 
are  a  significant  hazard  during  years  of  maximum  solar  activity.  Althou^  the  occurrence  of  large  flares 
is  not  yet  predictable,  especially  with  long  leads  times,  it  is  now  known  that  large  flares  lend  to  occur  in 
large,  highly  sheared,  magnetically-complex  active  regions;  so,  prediction  of  events  with  <24-h  notice 
ma'y  someday  be  possible. 

5.2. 1  Galactic  Cosmic  Ravs.  GCRs  consist  of  the  nuclei  of  the  elements  plus  about  2  percent 
which  are  electrons.  Their  ener^es  cover  the  range  from  10  MeV  per  particle  to  above  iOis  MeV  per 
particle.  The  flux  is  highest  during  periods  of  minimum  solar  activity.  At  these  energies,  the  nuclei  arc 
completely  ionized.  A  small  flux  of  x-rays  and  F-rays  from  celestial  sources  is  sometimes  inrfinVH  in 
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GCRs,  but  for  space  vehicle  design  purposes,  only  cosmic-ray  nuclei  aie  significant,  and  this  sectioa  wi 
be  ccufined  to  them.  GCRs  consist  of  the  nuclei  of  the  elements  from  hydrogen  through  iron  in  roughly 
the  same  proportions  as  are  found  in  the  solar  system,  but  with  the  heavier  nuclei  more  abundant  in  the 
cosmic  rays.  Nuclei  heavier  than  nickel  are  extremely  rare. 

Figure  5-6  gives  the  relative  abundances  and  energy  spectra  of  GCRs  of  interest.  Ihe  figure 
covers  the  energy  range  of  importance  in  essentially  aU  space  vehicle  radiation  effects  studies.  In  spite  ( 
their  small  number,  the  heavy  elements  are  very  important  due  to  dieir  densely  ionizing  tracks.  Th^  arc 
responsible  for  many  effects  in  detectors  and  microelectronics.  From  figure  5-6,  it  can  be  seen  that  the 
flux  of  each  nuclear  specie  decreases  rapidly  with  increasing  energy.  Tte  lowest  energies  are  observed 
outside  the  magnetosphere,  where  the  flux  is  limited  by  magnetic  fields  carried  by  the  solar  wind.  The 
energies  observed  and  the  flux  at  these  energies  vary  inversely  with  the  solar  cycle. 

Although  the  contribution  from  GCRs  to  the  total  dose  in  rads  inside  a  spacecraft  is  typicalfy  les 
than  15  percent  for  most  geocentric  orbits,  these  nuclei  are  responsible  for  such  effects  as  *^^s”  and 
“latch-up”  in  microcircuits  (large-scale  integrated  (LSI)  and  very  large-scale  integrated  devices 
(VLSIDs)).  Along  with  the  trapped  radiation-belt  protons,  the  nuclei  are  also  responsible  for  the  inducei 
radioactivity  in  most  materials  in  orbit  Noise  induced  directly  by  ionization  in  sensitive  devices  such  as 
charge-coupled  devices  (C(3Ds)  and  via  Cherenkov  and  fluorescence  radiation  in  photomultiplier  tubes 
are  other  effects  of  GCRs  that  must  frequently  be  considered.  The  designer  should  also  consider  the 
possible  effects  of  GCRs  on  materials  as  well  as  the  probability  of  production  of  secondary  particles  anc 
their  effects. 

The  available  data  have  been  fit  to  an  empirical  model  which  is  summarized  in  pan  1,  appendix 
B.  Figures  5-7a,  b,  and  c  show  the  differential  spectra  of  protons,  a-particles,  and  iron  nuclei  at  various 
levels  of  sclar  activity  and  the  fit  of  the  model.  This  model  should  be  used  to  define  the  GCR  environ¬ 
ment. 


5.2.2  Solar  Particle  Events.  Solar  particle  events  are  sporadic  {tenomena  with  durations  of  a  fev 
days.  They  occur  at  irregular  intervals  from  a  few  weeks  to  a  few  years  depending  on  the  intensity 
threshold  considered.  Their  most  significant  components  consist  (tf  protons  at  energies  from  a  few  Me  V 
to  a  few  hundred  MeV.  The  proton  spectra  may  exhibit  inten^ies  ^t  vary  over  several  orders  of  mag¬ 
nitude,  both  from  event  to  event  and  within  the  time  profile  of  an  individual  event  The  solar  particles 
will,  however,  envelop  die  Earth  within  minutes  after  a  solar  flare  event  reach  peak  intensity  in  a  few 
hours,  and  decay  in  1  or  2  days.  They  are  less  likely  to  occur  during  srdar  cycle  minima.  Within  the 
Earth’s  magnetosphere,  the  protons  reach  LEO  most  freely  in  die  polar  regions  at  magnetic  latitudes 
above  about  63*  because  the  magnetic  energetic  cutoff  goes  to  near  zero  at  higher  latitudes  (see  below). 

An  empirical  representation  of  the  probability  of  solar  particle  events  has  been  summarized  in 
part  2,  appendix  B.  Figure  5-8  shows  a  history  of  solar  proton  events  over  two  solar  cycles,  and 
figure  5-9  gives  differendai  spectra  integrated  over  various  events. 

5.2.3  Geomagnetic  Shielding  and  Cutoff  Rigidity.  The  Earth’s  magnetic  field  deflects  incoming 
cosmic  rays  (solar  and  galactic)  to  a  degree  which  depends  on  the  energy  of  the  particles,  preventing 
those  with  lower  energies  from  penetrating  deep  into  the  magnetoqihere  (cut-off  rigidity),  hi  LEO,  the 
minimum  rigidity,  P/Z,  that  a  particle  can  have  and  reach  a  latitude.  A,  fiom  the  vertical  is 


raHli  anH  J  •  ^  distaiicc  fioin  the  center  of  the  Earth  in  units  of  Earth 

radu,  an-  A  is  the  latitude.  The  minimum  cutoff  drops  to  near  zero  somewhat  faster  than  the  fnrmnin 

de,  as  can  be  seen  from  the  r/J?^  dependence.  Thus,  spacecraft  orbiting  the  Earth  may  either  be 
exposed  to  the  unattenuated  interplanetary  cosmic  ray  intensities  (e  g  high-altitude  ool 
Z  •^‘"‘^"ni^enWexposed (e'g.. low-altitSf '^roS)  or^ 

Swri^H  ^  V"  Of  approximately  10  MeV  for  low-inclination,  low-altitude 

e*^,m  ^  f  10  gJgaelectron  volt  (GeV)  for  equatorial  orbits.  To  approximate  the  average 

SSLf  ^  t^o^ic  rays,  calculations  must  be  performed  over  many  orbital  ^ 

revolutions.  An  approach  to  this  is  nmmarized  in  part  3.  appendix  B.  ^ 


Figure  5-6.  Cha^mtics  of  GCR&  Differential  energy  spectra  for  OCRs  outside  the  magnetosphere 
at  maximum  and  minimum  solar  activity  (as  defined  by  sunspot  number).  ^ 
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Hgure  5-la.  Differential  energy  spectrum  of  hydrogen  (mostly  protons).  Ihe  data  are  selected  to  sho\ 
tte  solar  maximum  and  solar  minimum  fluxes.  The  dashed  curve  is  a  worst-case  qiectmm. 
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Rgurc  5-7c.  Differential  energy  spectrum  of  iron.  Hie  data  are  selected  to  show  die  solar 
and  solar  minimum  fluxes.  The  dashed  curve  is  a  worst-case  !tp«»rtnim 


Figure  5-8.  Event-integrated  proton  fluxes  above  30  MeV  for  the  major  solar  events  of  the  19th  and  20th  solar  cycles. 
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Hguie  S-9.  Event-integrated  proton  differential  energy  spectra.  For  F^,  a  typical  ordinary  event 
Fmt.  a  worst-case  ordinary  event  (90-percent  confidence  level);  and  Fa, 
an  anomalously  large  solar  event 


VL  ELECTROMAGNETIC  RADIATION 


Important  sources  of  electromagnetic  noise  e.dst  over  the  entire  fnequency  s|iectrum  firom  direct 
rrent  (dc)  to  x-ray  at  low-Eaith  orbital  altitudes.  Hiese  noise  sources  broadly  squate  into  four 
tegories:  galactic,  solar,  near-Earth  natural  plasma,  and  manmade  radio  noise.  The  highest  power 
nsity  occurs  at  the  visible  and  IR  frequencies  and  amounts  to  the  solar  constant  of  137  C,oW/m2.  This 
ction  discusses  the  nature  of  these  noise  sources  and  the  manner  in  which  the  nataal  ionospheric 
vironment  may  affect  the  propagation  of  these  signals  to  the  q)ace  vehicle.  Only  natural  and  remote 
uimade  electromagnetic  radiations  (EMR’s)  are  considered. 

The  peak  of  the  ionosphere  with  an  electron  density  of  n,  =  10*°  to  10*^  m-^  occurs  at  200  to  300 
1  altitude  (see  section  IV  of  this  document).  These  densities  correspond  to  a  peak  in  the  natural 
:ctron  plasma  frequency  (Langmuir  oscillations)  oi  about  0.8  to  8  MHz.  EMR  bdow  this  fiequency 
mot  propagate  through  the  peak  electron  density  region  without  strong  attMinatinn  Thus,  sub- 
lospheric  and  ground-based  sources  below  this  frequency  wiQ  generally  not  reach  orbiting  space 
Glides  except  for  ducted  propagation  for  some  plasma  waves.  Note,  however,  that  the  density  and 
sition  of  this  peak  electron  density  vary  significandy  with  solar  activiQr  and  time  of  day. 

High  levels  of  power  in  the  RF  spectral  region  are  primarily  a  resnh  of  Earth-based  and  space 
tide  onboard  radiation  sources  that  are  manmade.  The  major  ground-based  transmitters  resulting  in 
;h  power  levels  in  the  ionosphere  are  in  the  100  MHz  to  5  GHz  range  used  by  radars.  They  generally 
/e  very  narrow  spatial  distributions. 


Galactic  Radio  Noise 

Galactic  noise  reaching  orbital  altitudes  extends  from  frequencies  of  IS  MIb  to  100  GHz.  The 
irce  is  broadly  directional  toward  and  perpendicular  to  the  galactic  plane  and  varies  py  about  10 
W/Hz.  Figure  6-1  is  a  spectrum  of  the  galactic  noise  as  a  function  of  frequency.  Space  vehicle  radio 
jineers  will  find  that  this  galactic  noise  is  a  dominant  factor  for  a  typical  radio  receiver  operating 
ween  40  and  250  MHz  (fig.  6-2).  For  comparison  between  cosmic  noise  levels  and  natural  near-Earth 
se.  see  section  6.3. 


Solar  ESectromagnetic  Radiation 

Solar  electromagnetic  energy  flux  near  the  Earth  but  outside  the  aimoqdierefat  1 AU)  is  about 
'1-10  W/m2.  The  irradiance  spectrum  approximates  a  black  body  curve.  T= 5762  K.  which  peaks  in 
visible  frequency  range  (fig.  6-3).  The  integrated  irradiance  of  the  nearUV  portion  of  the 
:tiomagnetic  flux  (0. 18  pm  <  A  <  0.4  pm)  is  approximately  1 18  W/m2.  Le.  about  8.6  percent  of  the 
u-  constant  The  far  UV  portion  (A  <  0. 18  pm)  contributes  about  0.023  W/ni2.  Inadiaiice  levels  for  the 
'erne  high  energy  spectrum  are  summarized  in  taMe  6-1. 

The  UV  portion  (A  <  0.3  pm)  of  the  electromagnetic  spectrum  is  (^particular  importance 
etermining  the  effects  of  solar  radiation  on  material  properties.  Fw  “tenestrial  space."  the  UV 
ation  is  independent  of  orbit  geometry  above  the  atmosphere;  however,  the  qraoecraft  «irfarp 
metry.  orientation,  and  orbit  geometry  will  affea  the  exposure  level  over  a  given  orbital  lifetime. 
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Figure  6-3 
outside 


'.  Normally  incident  solar  radiatioo  at  sea  level  on  veiy  clear  days,  solar  socctral  i 
the  Earth  s  atmosphere  at  1  AU.  and  blackbody  qiectial  irradianceoave  at  T= 5 

(nonnalized  to  1  Al^. 


5762K 


grcaUy  improved  the  data  base  for  this  portion  of  saectmm  Tliesc 

reference  22.  have  been  incorporated  into  this  n^^tioa 

?•'  ^  a  8^  indicaticm  of  the  imdiance  vanririlty  over  the  27-dhv 

ar  lomtion  penod  and  somewhat  longer  •ime  scales;  however,  the  time  ^ancovoed  by  eood  oan^ 

3kr  cycle.  confidence  in  estimates  of  inadiance^aiiaioK  over^tte  11-^^ 


Table  6-1.  Solar  electromagnetic  radiation. 


Type 

Waveiengdi(nm) 

Uvel(W^ 

LTV 

100-150 

7.5E-J 

EUV 

10-100 

2EJ 

X-Rays 

1-10 

5E-5 

Flare  X  Rays 

0.1-1 

lE-» 

Cosmic  Rays* 

- 

lE-V 

Solar  Wind* 

1E-* 

*Nol  elcciromagnetic  bm  inclaled  ft>  cnmpiin<KM  24 
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Table  6-2.  Parameters  for  estimating  irradiance  variability  over  the  1 1-year  solar  cyde.^ 


Wavelength  Range  (nm) 

i?5^(inax/imn) 

175-190 

1.020  (±0.020) 

1.07 

190-210 

1.030  (±0.020) 

1.06 

210-240 

1.026  (±0/)20) 

1.03 

240-300 

1.005  (±0.020) 

1.01 

is  the  ratio  of  a  27-day  mean  inadunce  near  solar  maxffliuui  to  a  27-day  stem  at  s(te  miniiiniiiL  is 

the  ratio  of  ibe  maximum  irradiance  observed  during  a  typical  27-day  solar  foiatioa  10  the  minimum 

inadiance  during  this  rotation. 

The  Lyman-alpha  (Ly-or)  emission  line  at  0.1216  pm  is  the  primary  emission  line  bdow 
0. 1300  pm.  and  its  integrated  irradiance  varies  by  roughly  a  factor  of  2  over  the  1  l-year  sofar  cycle  ariifa 
an  uncertainty  of  ±15  percent  For  estimates  of  the  Ly-aflux  at  any  level  of  solar  activity.  iKfoBow^ 
relationship  is  the  best  available: 


q(Ly-a)  =  2.25E“+0.014E**  (F,o.7-65)  (6-1) 

where  q  is  expressed  in  photons-cm-^  s-^. 

The  variation  in  the  integrated  0.130  to  0.175  pm  wavelength  flux  over  the  solar  cyckisafador 
of  2  with  error  bars  cn  the  recent  measurements  of  ±20  percent  The  percent  rhany  of  the  dm-cm 
variation  is  less  than  that  at  Ly-o;  and  most  values  lie  between  7.6E11  and  8.7Eii  photons  cs-^.^ 
solar  activity  levels  prevailing  in  1982  to  1984. 

It  is  now  generally  accepted  that  the  spectral  region  0.175  to  0.210  pm  varies  with  tlB27-dlay 
solar  rotation  and  the  1 1-year  solar  cycle.  In  this  spectral  region,  the  toot  mean  square  deviations  fiom 
the  Reference  Solar  Irradiance  Specnum  in  table  ^3  are  less  than  10  percent  in  most  cases  tie 

quoted  uncertainly  in  any  single  measurement  is  typically  ±15  to  ±20  percent. 

Typical  uncertainty  limits  for  the  measurements  in  the  0.210  to  0.3275  pm  spectral  r^kmare 
±15  to  20  percent;  however,  once  again  root  mean  square  deviations  from  the  reference  solarantiiame 
spectrum  ^SIS)  are  less  than  10  percent  Between  0.175  and  0.210  pm.  the  maximum  irracfimceieak 
to  be  1.05  to  1.07  limes  the  minimum  value  with  a  peak  ratio  ^pro^hing  1.10.  Betwemi  0.210  tad 
0.250  pm.  the  maximom  to  minimum  irradiance  ratio  is  typically  1.03  and.  b^ond  0.260  jim  Aty 
1.01  or  less.  Table  6-2  presents  estimates  of  the  solar  cycle  variability  in  selected  wavelength  bins 
spanning  the  region  0.175  to  0.300  pm.  The  quantity,)?^  is  ameasute  of  the  ratio  of  solar 
inadiance  to  solar  minimum  irradiance.  R27  ^  estimate  of  the  27-day  irradtance  variation.  Ihe  atm  of 
the  maximum  to  minimum  irradiance  during  a  solar  cycle  is  the  product  of  R27  and  ^se- 
and  minimum  iiradiaiices  over  a  solar  cycle  are  related  to  the  mean  reference  ctycle  of  table  6-3  bf: 

.  »2) 

where /I27  and  are  Obtained  from  table  6-2. 


Table  6-3.  Reference  solar  iiradiance,  Rayleigb  scattering,  oxygen,  ozooe  cross  secticn^ 

and  energy  fhix  (page  1  of  6).* 


Wavelength  Ranget 
(nm) 

Solar 

Irradiance 

(phot/cm^-s) 

Raylei^ 
Scalier  ajts 

(an2) 

Oxygen 

Hei^rg 

oHz(02) 

(cm2) 

Ozonctt 

0(Q3) 

(cm^) 

Eoe^ 

FbK 

(Wfefi 

1  175.439-176.991 

1.74E-hll 

6.79E-25 

4.61E-24 

8.11E-19 

IS6E-3 

2  176.991-178.571 

2.10E+11 

6.49E-25 

5.03E-24 

7.99E-19 

234E-3 

3  178.571-180.180 

2.38E+11 

6.20E-25 

5.46E-24 

7.86E-19 

2.64E-3 

4  180.180-181.818 

3.04E+11 

5.93E-25 

5.88E-24 

7.63E-19 

335B-3 

5  181.818-183.486 

3.19E+11 

5.66E-25 

6.29E-24 

7.29E-19 

3.47E-3 

6  183.486-185.185 

2.93E+11 

5.41E-25 

6.68E-24 

6.88E-19 

3.16E-3 

7  185.185-186.916 

3.62E+11 

5.16E-25 

7.04E-24 

6.40E-19 

3.87E-3 

8  186.916-188.679 

4.73E+11 

4.93E-25 

7.36E-24 

5.88E-19 

5.O0E-3 

9  188.679-190.476 

5.61E+11 

4.70E-25 

7.64E-24 

5.31E-I9 

5.88E-3 

10  190.476-192.308 

6.63E+11 

4.49E-25 

7.87E-24 

4.80E-19 

6.88E-3 

11  192.308-194.175 

6.90E+11 

4.28E-25 

8.04E-24 

4.38E-19 

7.09E-3 

12  194.175-196.078 

9.56E+11 

4.08E-25 

8.14E-24 

4.11E-19 

9.74E-3 

13  196.078-198.020 

1.15E+12 

3.89E-25 

8.17E-24 

3.69E-19 

1.16E-2 

14  198.020-200.000 

1.27E+12 

3.71E-25 

8.13E-24 

3.30E-19 

1.27E-2 

15  200.000-202.020 

1.52E+12 

3.53E-25 

8.0IE-24 

3.26E-19 

ISOE^ 

16  202.020-204.082 

1.78E+12 

3.36E-25 

7.84E-24 

3.26E-19 

1.74E3 

17  204.082-206.186 

2.20E+12 

3.20E-25 

7.63E-24 

3J1E-19 

2.13E2 

18  206.186-208.333 

2.69E+12 

3.05E-25 

7.33E-24 

4.11E-19 

2.S8E2 

19  208.333-210.526 

4.54E+12 

2.90E-25 

6.99E-24 

4.84E-19 

431E^ 

20  210.526-212.766 

7.14E+12 

2.76E-25 

6.45E-24 

6.26E-19 

6.7DE2 

21  212.766-215.054 

8.35E+12 

2.62E-25 

5.81E-24 

8J7E-19 

7.75E3 

22  215.054-217.391 

8.39E+12 

2.49E-25 

5.23E-24 

1.17E-18 

7.71E3 

23  217.391-219.780 

1.08E+13 

2.36E-25 

4,71E-24 

1.52E-18 

9.S2E-2 

24  219.780-222.222 

1.18E-I-13 

2.24E-25 

4.26E-24 

1.97E-18 

1.06E-1 

25  222.222-224.719 

1.60E+13 

2.13E-25 

3.80E-24 

2J5E-18 

1.4QE-1 

26224.719-227.273 

1.34E+13 

2.02E-25 

3J5E-24 

3J4E-18 

1.SSE-: 

27  227.273-229.885 

1.41E+13 

1.92E-25 

2.90E-24 

4.00E-18 

1.23E-1 

From  reference  22. 

Wavelength  range  for  spectral  intervals  1  to  49  conc^nds  to  S(X)  wavenombers.  Wavdin^ianse 
for  remainder  of  specti^  intervals  is  5  nm. 

7’=203Kandr=273K. 


Table  6-3.  Reference  solar  iiradiance,  Rayleigh  scatterii^,  oxygen,  ozone  cross  sectiaas, 

and  eneigy  flux  (page  2  of  6).* 


Wavelength  Ranget 
(nm) 

Solar 

Irradiance 

(phot/cm^-s) 

Rayleigh 
Scatter  0)t5 
(cm2) 

Oxygen 

Heizberg 

<Jlfa(02) 

(cm2) 

Ozonett 

0(03) 

(cm2) 

Energy 

FhK 

(W/tai^ 

28229.885-232.558 

1.57E+13 

1.82E-25 

2.45E-24 

4.83E-18 

1J5E-1 

29232.558-235.294 

1.38E+13 

1.72E-25 

2.05E-24 

5.79E-18 

L17E-I 

30235.294-238.095 

1.60E+13 

1.63E-25 

1.69E-24 

6.86E-18 

L34E-1 

31 238.095-240.964 

1.45E+13 

1.54E-25 

lJOE-24 

7.97E-18 

L20B-1 

32240.964-243.902 

2.20E+13 

1.46E-25 

0.93E-24 

9.00E-18 

L80E-1 

33  243.902-246.914 

1.99E+13 

1.38E-25 

O.OOE-00 

l.OOE-17 

L61E-1 

34246.914-250.000 

1.97E+13 

1.31E-25 

- 

1.08E-17 

L58E-1 

35250.00&-253.165 

1.94E+13 

1.23E-25 

- 

1.13E-17 

1.53E-1 

36253.165-256.410 

2.91E+13 

1.17E-25 

- 

1.15E-17 

2.27E-1 

37256.410-259.740 

4.95E+13 

l.lOE-25 

- 

1.12E-17 

3.81E-1 

38259.740-263.158 

4.53E+13 

1.04E-25 

— 

1.06E-17 

3l44E-1 

•  From  reference  22. 

t  Wavelength  range  for  spectral  intervals  1  to  49  corresponds  to  SOO  wavenumbos.  Wasdenedi  range 
for  remainder  of  specti^  intervals  is  5  nm. 
ttT=203Kandr=273K. 
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Table  6-3.  Reference  solar  irradiance,  Rayleigh  scattering,  oxygen,  ozone  cross  sections. 

and  energy  flux  (page  3  of  6).* 


Wavelength  Ranget 
(nm) 

Solar 

Irradiance 

(phot/cm^-s) 

Rayleigh 
Scatter  <trs 
(cm2) 

Ozone 
0(03) 
7=203  K 
(cm2) 

Ozone 

0(03) 

7  =  273K 
(cm2) 

Energy 

Flux 

39  263.158-266.667 

1.07E+14 

9.78E-26 

9.59E-18 

9.65E-18 

8.02E-1 

40  266.667-270.270 

120E+14 

9.22E-26 

8.31E-18 

834E-18 

8.88E-1 

41  270.270-273.973 

I.IOE+H 

8.68E-26 

6.89E-18 

6.92E-18 

8.03E-1 

42  273.973-277.778 

1.04E+14 

8.17E-26 

5.35E-18 

5.42E-18 

7.49E-1 

43  277.778-281.690 

8.24E+13 

7.68E-26 

3.91E-18 

4.02E-18 

5.85E-1 

44  281.690-285.714 

1J2E+14 

7.22E-26 

2.67E-18 

2.77E-18 

1.06 

45  285.714-289.855 

1I5E+14 

6.78E-26 

1.73E-18 

1.79E-18 

1.48 

46  289.855-294.118 

3.48E+14 

6.36E-26 

1.04E-18 

1.09E-18 

2L37 

47  294.118-298.507 

3.40E+14 

5.97E-26 

5.85E-19 

6.24E-19 

228 

48  298.507-303.030 

322E+14 

5.59E-26 

3.16E-19 

3.43E-19 

2.13 

49  303.030-307.692 

4J3E+14 

5.24E-26 

1.66E-19 

1.85E-19 

2.75 

50  307.692-312.5 

4.95E+14 

4.90E-26 

8.67E-20 

9.80E-20 

3.17 

51  312.5-317.5 

5.44E+14 

4.58E-26 

4.33E-20 

5.01E-20 

3.43 

52  317.5-322.5 

5.93E+14 

4.28E-26 

2.09E-2O 

2.49E-20 

3.68 

53  322.5-327.5 

6.95E+14 

4.01E-26 

9.37E-2I 

1.20E-20 

425 

54  327.5-332.5 

8.15E+14 

3.75E-26 

4.7IE-21 

6.17E-21 

4.91 

55  332.5-337.5 

7.81E+14 

3.52E-26 

1.98E-21 

2.74E-21 

4i>3 

56  337.5-342.5 

835E+14 

3.31E-26 

7.77E-22 

1.17E-21 

428 

57  342.5-347.5 

8.14E-t-I4 

3.11E-26 

1.77E-22 

S.88E-22 

429 

58  347.5-352.5 

833E+14 

2.92E-26 

- 

2.66E-22 

424 

59  352.5-357.5 

9.i7E+14 

2.75E-26 

- 

L09E-22 

5.13 

60  357.5-362.5 

838E+14 

2.60E-26 

- 

5.49E-23 

422 

61  362.5-367.5 

1j04E-»-15 

2.45E-26 

- 

- 

526 

62  367.5-372.5 

l.lOE+15 

2.31E-26 

j 

- 

521 

63  372.5-377.5 

9.79E+14 

2.19E-26 

- 

5.19 

64  377.5-382.5 

1.I3E+15 

2.07E-25 

- 

5.91 

*  From  reference  22. 

t  Wavelength  range  for  spectral  intervals  1  to  49  corresponds  to  500  wavenumbers.  Wavdenglh  range 
for  remainder  of  spectr^  intervals  is  S  nm. 
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Table  6-3.  Reference  solar  uradiance,  Rayleigh  scattering,  oxygen,  ozone  cross  sections, 

and  energy  flux  (page  4  of  6).* 


Wavelength  Rangef 
(nm) 

Solar  Irradiance 
(phot/cm2-s) 

Rayleigh  Scatter 
0)15  (cm2) 

Ozonett  0(03) 

(cm2) 

Energy  Flux 
(W/m2) 

65  3825-387.5 

8.89E-»-14 

1.96E-26 

- 

4.59 

66  3875-392.5 

1.14E+15 

1.86E-26 

— 

5.81 

67  3925-3975 

9.17E+14 

1.76E-26 

4.61 

68  3975^2.5 

1.69E+15 

1.67E-26 

— 

8.39 

69  402.5-407.5 

1.70E+15 

1.59E-26 

8.34 

70  4075^12.5 

1.84E+15 

151E-26 

2.91E-23 

8.91 

71  4125-417.5 

1.87E+15 

1.44E-26 

3.14E-23 

8.95 

72  417.5-422.5 

1.95E+15 

1.37E-26 

3.99E-23 

9.22 

73  422.5-427.5 

1.81E+15 

1.30E-26 

6.54E-23 

8.46 

744275-432.5 

1.67E+15 

1.24E-26 

6.83E-23 

7.72 

75  432.5-437.5 

1.98E+15 

1.18E-26 

8.66E-23 

9.04 

76  437.5^2.5 

2.02E+15 

1.13E-26 

1.25E-22 

9.12 

77  442.5-447.5 

2.18E-t-15 

1.08E-26 

1.49E-22 

9.73 

78  447.5-452.5 

2.36E+15 

1.03E-26 

1.71E-22 

1.04E+1 

79  452.5-457.5 

2.31E+15 

9.85E-27 

1  2.12E-22 

l.OlE+1 

80  457.5-462.5 

2.39E+15 

9.42E-27 

3.57E-22 

1.03E+1 

81  462.5-467.5 

2.38E+15 

9.01E-27 

3.68E-22 

1.02E+1 

82467.5-472.5 

2.39E+15 

8.63E-27 

4.06E-22 

l.OlE+1 

83  4725-477.5 

2.44E+15 

8.26E-27 

4.89E-22 

1.02E+1 

84477.5-482.5 

2.51E+15 

7.92E-27 

7.11E-22 

1.04E-)-l 

85  482.5-487.5 

2.30E+15 

7.59E-27 

8.43E-22 

9.42 

86  487.5-492.5 

2.39E+15 

7.28E-27 

8.28E-22 

9.69 

87  492.5-497.5 

2.48E-I-15 

6.99E-27 

9.09E-22 

9.95 

88  497.5-502.5 

2.40E+15 

6.71E-27 

1.22E-21 

954 

89  5025-5075 

2.46E+15 

6.44E-27 

152E-21 

9.68 

90  507.5-512.5 

2.49E+15 

6.19E-27 

1.58E-21 

9.70 

91  512.5-5175 

2.32E+15 

5.95E-27 

1.60E-21 

8.95 

92  5175-522.5 

2.39E+15 

5.72E-27 

1.78E-21 

9.13 

93  522.5-5275 

2.42E+15 

5.50E-27 

2.07E-21 

9.16 

94  527.5-532.5 

2.55E+15 

5.30E-27 

2.55E-21 

956 

95  532.5-537.5 

2.51E+15 

5.10E-27 

2.74E-21 

9.32 

*  From  iirfereiice  22. 

t  Wavelength  range  for  ^wctial  intervals  1  to  49  corresponds  to  5(X)  wavenumbers.  Wavelength 
range  fm  remainder  of  spectral  intervals  is  5  run. 
ttT  =  203Kandr=273K. 


Table  6-3.  Reference  solar  irradiance,  Rayleigh  scattering,  oxygen,  ozone  cross  sections. 

and  energy  flux  (page  5  of  6).* 


6- 


Wavelength  Ranget 
(nm) 


96  537.5-542.5 


97  542J-547.5 


98  547 J-552.5 


99  552.5-557.5 


100  557 J-562.5 


101  562.5-567.5 


102  567.5-572.5 


103  572.5-577.5 


104  577.5-582.5 


105  582.5-587.5 


106  587.5-592.5 


107  592.5-597.5 


108  597.5-602.5 


109  602.5-607.5 


I10607.5-<il2.5 


111  612.5-^17.5 


112617.5-622.5 


113622.5-627.5 


114  627.5-632.5 


115  632.5-637.5 


116637.5-642.5 


117  642.5-647.5 


118  647.5-652.5 


119  652.5-657.5 


120  657.5-662.5 


121  662.5-667.5 


122  667.5-672.5 


123  672.5-677.5 


124  677.5-682.5 


125  682.5-687.5 


126  687.5-692.5 


Solar  Irradiance 
Q)hot/cm2-s) 

Rayleigh  Scatter  ( 
ORSicrn^) 

2.49E+15 

4.91E-27 

2.55E+15 

4.73E-27 

2.53E+15 

4.56E-27 

2.54E+15 

4.34E-27 

2.50E+15 

4.18E-27 

2.57E+15 

4.04E-27 

2.58E+15 

3.90E-27 

2.67E+15 

3J6E-21 

2.67E+15 

'  3.63E-27 

2.70E+15 

3.51E-27 

2.62E+15 

3.39E-27 

2.69E+15 

3.28E-27 

2.63E+15 

3.17E-27 

2.68E+15 

3.06E-27 

2.66E+15 

2.96E-27 

2.59E+15 

2.87E-27 

2.69E+15 

2.77E-27 

2.61E+15 

2.68E-27 

^62E+15 

2.60E-27 

2.62E+15 

2.52E-27 

i63E+15 

2.44E-27 

2.60E+15 

2.36E-27 

2.55E+15 

2.29E-:7 

2.48E+15 

2.22E-27 

2.57E+15 

2.15E-27 

2.61E+15 

2.09E-27 

2.61E+15 

2.03E-27 

2.62E+15 

1.97E-27 

2.62E+15 

1.91E-27 

2.57E+15 

1.85E-27 

(cn»2) 


2.88E-21 


3.07E-21 


3.17E-21 


3.36E-2! 


3.88E-21 


4.31E-21 


4.67E-21 


4.75E-21 


4.55E-21 


4.35E-21 


4.42E-21 


4.6iE-21 


4.89E-21 


4.84Er-21 


4.54E-21 


4.24E-21 


2.52E+15 


1.80E-27 


3.43E-21 


3.17E-21 


2.74E-21 


2.61E-21 


2.42E-21 


2.20E-21 


2 


1.85E-21 


1.67E-21 


1.54E-21 


1.42E-21 


1.25E-21 


1.12E-21 


Energy  Hux 
(W/ni2) 


9.16 


9.30 


9.14 


9.09 


8.87 


9.04 


8.99 


9.22 


9.14 


9.17 


8.82 


8.98 


8.71 


8.80 


8.66 


8.37 


8.62 


8.30 


8.26 


8.20 


8.17 


8.01 


7.79 


7.52 


7.73 


7.80 


7.74 


7.71 


7.66 


7.45 


I  7.25 


*  Rom  reference  22. 

t  Wavelength  range  for  ^lectial  intervals  1  to  49  corresponds  to  500  wavenumbers.  Wavelength 
range  for  remainder  of  spectral  intervals  is  5  n- 
ttr=203Kandr=273K. 
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Table  6-3.  Reference  solar  iiradiance,  Rayleigh  scattering,  oxygen,  ozone  cross  sections, 

and  energy  flux  (page  6  of  6).* 


Wavelength  Ranget 
(nm) 

Solar  Irradiance 
(phot/cm2-s) 

Rayleigh  Scatter 

0)JS  (cm2) 

Ozonett  0(03) 

(cm2) 

Energy  Flux 
(W/ni2) 

127  692.5-697J 

2.60E+15 

1.75E-27 

1.02E-21 

7.43 

128  697.5-702J 

2.58E+15 

1.70E-27 

9.20E-22 

732 

129  702.5-707J 

2.52E+15 

1.65E-27 

8.40E-22 

7.10 

130  707.5-712J 

2.51E+15 

1.60E-27 

7.70E-22 

7.02 

131  712J-717J 

2.48E+15 

1.56E-27 

6.90E-22 

6.89 

132  717.5-722J 

2.45E+15 

1.52E-27 

6.30E-22 

6.76 

133  722.5-727.5 

2.48E+15 

1.47E-27 

5.70E-22 

6.79 

134  727.5-732.5 

2.45E+15 

1.43E-27 

5.25E-22 

6.67 

135  732.5-737J 

2.44E+15 

i.39E-27 

4.75E-22 

6.60 

136  737.5-742J 

2.39E+15 

1.36E-27 

4.47E-22 

6.41 

137  7423-747J 

2.40E+15 

1.32E-27 

4.20E-22 

6.40 

138  747.5-752.5 

2.41E+15 

1.29E-27 

3.75E-22 

638 

139  752.5-757.5 

2.40E+15 

1.25E-27 

3.25E-22 

632 

140  757.5-762J 

2.38E+15 

1.22E-27 

2.92E-22 

632 

141  762.5-767J 

2.34E+15 

1.19E-27 

2.76E-22 

6.08 

142  767.5-772J 

2.32E+15 

1.16E-27 

2.70E-22 

5.99 

143  772.5-777J 

2.30E+15 

1.13E-27 

2.80E-22 

5.89 

144  777.5-782J 

2.33E+15 

l.lOE-27 

2.85E-22 

5.93 

145  782.5-787J 

2.34E+15 

1.07E-27 

2.52E-22 

5.92 

146  787 J-792J 

2.29E+15 

1.04E-27 

220E-22 

5.76 

147  7915-797J 

2.29E+15 

1.02E-27 

1.82E-22 

5.72 

148  797J-802J 

2.27E+15 

9.90E-28 

1.63E-22 

5.64 

149  802J-807J 

2.27E+15 

9.66E-28 

1.75E-22 

5.60 

150  807J-812J 

2.20E+15 

9.42E-28 

1.90E-22 

5.40 

151  812.5-817.5 

2.22E+15 

9.19E-28 

1.85E-22 

5.41 

152  817.5-822J 

2.18E-I-15 

8.96E-28 

1.70E-22 

538 

153  822.5-827J 

2.20E+15 

8.75E-28 

1.52E-22 

530 

154  827.5-832J 

2.14E+15 

8.53E-28 

1.42E-22 

5.12 

155  832J-837J 

2.14E+15 

8.33E-28 

1.40E-22 

5.09 

156  837.5-842J 

2.13E+15 

8.13E-28 

1.40E-22 

5.04 

157  842.5-847.5 

2.09E+15 

7.94E-28 

1.42E-22 

4.91 

158  847J-8525 

2.05E+15 

7.75E-28 

1.45E-22 

4.79 

*  From  reference  22. 

t  Wavelength  range  for  spectral  intervals  1  to  49  corresponds  to  5(X)  wavenumbers.  Wavdength 
range  for  remainder  of  spectral  intervals  is  S  nm. 
ttr  =  203  K  and  7=273  K. 


Natural  Environmental  Electromagnetic  Radiation 

J*‘®™;S"®*osphe';p-wnosphere  system  is  filled  with  natural  plasma  emission  sources  which 
equator.  The  frequencies  of  these  emissions  generally  extend  down  tro; 
Ae  ireak  electron  ptoma  frequencies  (1  to  10  MHz)  to  direct  current.  Die  electro^  plasma  frequenc 
found  simply  from  the  electron  density  as:  nequenc 

and  m,  is  U,. 

^  electron  plasma  frequency  is  independent  of  temperature.  A  wide  variety  of  nai 
pSh  emissions  ^related  to  electron  cyclotron  waves  from  wave-particle  phenomeiL  in 

.w?!i  f  emissions  are  generally  in  the  frequency  Lige  of  0  5  to  30  kP? 

mclude  very  low-frequency  whistler-mode  noise  such  as  hii  and  chores.  TTiey^  Wghly  Son 
dqiendent  since,  although  they  are  electromagnetic  waves,  they  have  a  plasma  compiwnt  which  aff 
the  wave  propagation  vector.  At  typical  LEO  velocities  in  the  5  to  10  km/s  range,  one  can  expect  nat 
ptema  structures  ^small  scale  size  (i.e..  O.OI  to  10  km)  to  introduce  an  effective  noise  so^  at  so 

SiniaSSrF?’  ^1?  phenomena  such  as  auroral  arcs  and  ionospheric 

inegulanues.  Hgure  6-4  presents  an  overview  of  the  natural  plasma  noise  levels  from  near-EarSi 
^  cosm«:  sounds.  Hie /.gure  shows  pom  nux  lewis  for  ,^ous  fr^w^^SS 

M  *^P'““*"™‘"lheEarth'sedvlronioen?sodta  as^Awicsl 

1 

As  mentioned  above,  noise  sources  below  about  1  MHz  are  Hkelv  m  iv»  tiAoiimkiw  «» -« 
wtele  if  ihey  haw  sources  below  200  bn.  However  “  L^hb  fSrST^'USMX; 
wh  as  li^uing.generaled  sferics  lo  penetrate  the  ioiiosphere.  In  particular,  recent  iniliniiinnt  are  ih 
*te«  hghlning-getieraled  waveforms  with  frequencies  above  aboutl  tamS.  !5!  T 

Helds  up  to  about  50  mV/m  in  the  ionosphere.  In  geS!^TSfteld  s^SeS^C£r  a  tS 
Sin '?vr?'*K^  lai^t  Held  strengths  are  for  quasi^slalic  Held  .structures  whteh  occasionallv  reach  0 

“ r  '***^  '« >» “mv/m.  such  large  quasi-staticTKlSS^L, 

poleward  of  50"  lautude  in  associaUon  with  the  aurora.  Hiese  dLci  cnn<aIrn.M.  .t,™ii!2™  “ 
the  apparent  direa  ^nt  elecuic  Held  due  to  the 

6.4  Manmade  Noise 

ri,  1  . manmade  noise  sources  below  the  peak  plasma  frequency  of  1  to  10  MHz  will  nor 

foJTI  MHz attenuation,  it  is  likely  fahat  many  nan^u^btmd  rourcc 
to  300  GHa  wtll  be  detectable  and  perhaps  important  noli  somees  tor 
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and  plasma  waves.** 
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The  electric  fleld  E  (m V/m)  at  a  distance  z  (kilometers)  from  any  RF  transmitter  can  be  estimated 
with  the  formula: 


£  =  (30ERP)ifl/z  .  (5.5) 

where  ERP  (waits)  is  the  effective  radiated  power  which  equals  the  product  of  the  transmitter  power  and 
the  antenna  gain. 

Thus,  a  radar  with  a  power  of  10  kW  and  an  antenna  gain  of  40  dB  (10  000)  will  lesuh  in  fields 
of  over  100  mV/m  at  500  km  altitude.  Presumably  a  source  such  as  a  radar  will  have  a  very  narrow  tw»a«n 
width  (perhaps  less  than  a  degree)  but  will  be  wept  over  wide  loci  of  look  tfirections.  Occasionally, 
ionosphere  research  transmitters  are  also  operated  which  produce  measurable  »■»«*««;  of  the  kmospberk 
plasma  and  have  power  levels  in  the  same  range  as  the  defense  radars  but  wider  beam  patterns. 

In  June  1983.  the  U.S.  Air  Force  Space  Division  requested  the  Department  of  Defense  Electro- 
magnetic  Compatibility  Analysis  Center  (ECAC)  to  identify  the  worldwide  RF  electromagnetic 
environment  (EME)  for  hypothetical  spacecraft  at  various  LEO’s.  The  lesahs  of  this  task  are 
documented  in  the  “Worldwide  Spacecraft  EME  Definition.”  report  number  ECAC-CR-85-065.  USAF 
Space  Division  (refs.  B  and  I). 


Vn.  METEOROIDS  AND  ORBITAL  DEBRIS 


Once  in  orbit,  a  space  vehicle  will  encounter  meteoroids  and  orbital  debris.  Either  object  can 
pose  a  serious  damage  or  decompression  threat  to  the  vehicle.  Meteoroids  are  natural  in  origin,  and 
debris  is  the  result  of  manmade  material  remaining  in  Earth  orbit  General  information  and  discussion 
appears  in  references  26  and  27.  The  debris  model  presented  here  has  been  recommended  for  most  if 
not  all.  NASA  engineering  applications  since  mid-1990.  Comments  have  been  added  where  recent 
measurements  or  studies  have  confirmed  or  raised  issues  with  model  results.  Hie  model  is  still 
recognized  as  the  best  available  and  valid  within  stated  uncertainties;  no  improvements  or  iqidaies  are 
expected  in  the  near  future. 

For  historical  reasons  related  to  measurement  method,  the  meteoroid  and  debris  environments  are 
usually  specified  as  a  time-averaged  flux,  Fr.  against  a  single-sided,  randomly  tumbling  surface.  Flux  is 
defmed  as  number  of  intercepted  objects  per  unit  time  and  area.  For  Fr,  the  relevant  area  is  the  actual 
surface  area  of  the  satellite.  One  may  also  define  a  “cross-sectional  area  flux,”  Fc  for  a  randomly 
tumbling  satellite,  where  the  relevant  area  is  the  time-averaged  cross  sectional  area.  A  useful  theorem 
which  is  obvious  for  a  tumbling  sphere  but  which  holds  for  objects  of  arbitrary  non-self-shielding  shape 
(no  concave  surfaces),  is  that  Fc  =  4  Fr. 

For  spacecraft  which  fly  with  a  fixed  orientation,  the  meteoroid  and  orbital  debris  fluxes  are 
treated  as  vector  quantities,  F,  and  the  effects  of  directionality  must  be  carefully  evaluated.  Some  effects 
of  impact  will  be  direction  dependent.  To  simply  evaluate  the  expected  number,  N,  or  probabili^  of 
impacts  from  either  meteoroids  or  debris  (or  both),  one  may  use  a  “il-factoi”  method  and  the  approfviate 
Fr,  such  that: 


f  Zk,F^,d,  . 


(7-1) 


where  the  summation  is  over  the  i  surfaces  of  the  spacecraft,  each  of  area  Aj,  andl^^Fris  the  amiai  flux 
on  surface  A,-.  The  calculation  of  kj  is  discussed  in  section  7.3. 


Once  an  A/  has  been  determined,  the  proba'  ility  of  exactly  n  impacts  occurring  on  a  surfiure  in  the 
corresponding  time  interval  is  found  from  Poisson  statistics,  thus; 


ni 


(7-2) 


7.1  Meteoroids 

The  meteoroid  environment  encompasses  only  particles  of  natural  origin.  Nearly  all  meteoroids 
originate  from  comets  or  asteroids.  Meteoroids  that  retain  their  parent  body  orbit  and  create  periods  of 
high  flux  are  called  streams.  Random  fluxes  with  no  apparent  pattern  are  called  sporadic. 

The  average  total  meteoroid  environment  presented  is  comprised  of  the  average  qioradic 
meteoroids  and  a  yeariy  average  of  stream  meteoroids.  The  mass  density  for  meteoroids  spans  a  wide 
range,  from  approximately  0.2  g/cm3  or  less  for  a  portion  of  the  population  to  values  as  large  as  8  gfcm3. 
The  values  for  average  mass  density  vary  widely,  so  that  a  value  can  only  be  estimated.  Recommended 


mean  values  are  2  g/cm^  for  meteoroids  smaller  than  10-<  g.  1  g/cnP  for  meteoroids  between  10^  and 
0.01  g.  and  0.5  g/cm3  for  masses  above  0.01  g.  The  uncertainty  in  density  is  not  too  serious  because  the 
model  presented  below  was  derived  from  crater  and  impact  data  and  will  provide  a  good  representation 
of  expected  damage  even  though  absohite  mass  calibration  could  conceivaUy  be  in  error  by  as  much  as 
a  factor  of  10  for  the  smallest  sizes. 

Because  of  the  precession  of  asatellite’s  orbit  and  the  tilt  of  the  Eardi*s  equatorial  plane  with 
respect  to  the  ecbptic  plane,  the  meteoroid  environment  can  be  assumed  to  be  mnniditectional  relative  to 
Earth  for  design  applications.  Ifowever.  it  becomes  directional  relative  to  a  qncecraft  moving  through 
the  environment  with  most  meteoroids  coming  from  the  direction  of  motion.  The  directionality  derives 
from  the  vector  summation  of  the  spacecraft  velocity  vector  with  the  meteon^  velocity  (fistribution.  An 
additional  directionality  factor  is  introduced  by  the  shielding  {m)vided  by  Earth. 

The  normalized  meteoroid  velocity  distribution  with  respea  to  Earth  is  illustrated  by  figure  7>1. 

It  is  given  by  the  expressions  (number  per  km/s): 

n(v)  =  0.112,  11.1  S  v<  16.3  km/s, 

«(v)  =  3.328x10*  v-5^.  16.3Sv<55km/^  (7-3) 

n(v)  =  1.695xl0-»,  55  ^  v  S72.2  km/s. 

This  distribution  has  an  average  velocity  of  about  17  km/s;  relative  to  an  orbiting  spacecraft  the 
average  velocity  is  about  19  km/s.  To  determine  the  velocity  and  direction  distribotions  relative  to  any 
surface  on  an  orbiting  .^cecraft,  the  vector  relationship  between  the  meteoroid  velociqr  and  the 
spacecraft  velocity  should  be  used  as  Ascussed  in  paragra{rfi  7.3. 


Figure  7-1.  Normalized  meteoroid  velocity  distribution  from  eqpiation  (7-3). 


Meteoroid  flux  is  given  in  terms  c.  the  integral  flux  Fr,  the  number  of  particles/m^/year  of  mass 
m  or  greater  against  a  randomly  tumbling  surface.  Fv,.  the  interplanetary  f.ux  at  1  au  is  described 
mathematically  as  follows  (for  #n  :s  10  g): 

F'>,{iw)  =  Co  {  {CiillOJ06+C2H-38+C3(m+C4m2+C5m4)-fl-36  +C6(m+C7m2)-0.85  ) ,  (7. 

where: 

Co  =  3.156x102 
c,  =  2.2x103, 

C2=15. 

C3  =  1.3x10-9 
C4=10n. 

C5=  1027. 

C6=  1.3x10-16 
C7=106. 

To  convert  the  meteoroid  flux.  Fv^,  stated  above  to  that  in  Earth  orbit.  Fr,  both  Earth  shielding 
and  focusing  factors  must  be  applied,  Fr  =  s/Gf  FVr.  The  formula  for  shielding  is: 

shielding  factor  =  S7^=(l+cos  17 )/2  ,  (7.57 

where 


sm  »7  =  ReHRe^H)  , 

where  Re = Earth  radius  + 100  km  atmosphere  (6478  km)  and  =  height  above  Earth’s  atmosphere 
(Earth’s  atmosphere  height  to  be  taken  as  100  km  for  this  purpose).  Consequently,  the  Earth  shielding 
factor  vanes  from  0.5  just  above  tte  atmosphere  to  1.0  in  deep  space. 

The  factor  Ge  represents  the  focusing  effect  of  Earth’s  gravitational  field  which  attracts 
meteoroids  smd  increases  their  flux.  The  factor  ranges  from  a  value  of  2.0  just  above  the  atmosphere  to  a 
value  of  1.0  m  deep  space.  The  focusing  factor  is  represented  by  the  following  equation: 

focusing  factor  =  Gf  =  l+(F£/r)  (7^ 

where  Re = Earth  radius  + 100  km  atmosphere  (6478  km)  and  r  =  orbit  radius.  The  meteoroid  environ- 
ment  at  SOO  km  is  illustrated  m  figure  7-2. 

7.1.1  UKeitainiY  m  He  Meieoroid  Environment  Except  for  small  msmir  Hnc^  grainc 
OTflected  from  the  stratosphere,  tbe  physical  properties  of  meteoroids  must  be  determined  by  relatively 
uidirect  means,  examination  of  impact  craters,  optical  scattering,  etc.  They  are  also  known  to 
from  comets  (apparenfly  composed  of  low  density  “ices”  and  dust)  and  asteroids  which  are  rock-Bke 
Therefore,  there  is  considerable  uncertainty  in  their  properties.  In  particular,  the  uncertainty  in  mass 
lends  to  dominate  the  uncertainties  in  the  flux  measuremeuL  For  meteoroids  less  dian  10^  g  the  is 
uncertam  to  within  a  factor  fium  about  0.2  to  5  times  the  estimated  value,  which  impUes  the  flux  is 
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ncertain  to  within  a  factw  of  0.33  to  3  at  a  given  mass.  For  meteoroids  above  this  size,  the  flux  ii  weU 
eflned  but  the  associated  mass  is  even  mote  uncertain.  This  implies  an  effective  oocenainty  in  Ae  flux 
U  a  set  mass)  of  a  factor  from  0.1  to  10  (because  of  the  slope  of  the  function^  rdatioiidup)L 


DUnMlw(Gai) 


Figure  7*2.  Comparison  of  meteoroid  and  orbital  debris  fluxes,  asa  fimctioa  of  size. 

2  Orbital  Debris 

7.2.1  fiackgQMUl- The  natural  meteoroid  flux  discussed  above  lepneseats.  at  any  instant,  a  total 
about  200  kg  of  mass  widiin  2000  km  of  the  Earth’s  surface,  most  of  it  conoentraied  in  the  0. 1  mm 
:teoroids.  Within  this  same  2000  km,  there  is  an  estimated  1.5  to  3  million  1^  of  niamnade  ortiimg 
jects  as  of  mid*1988.  Most  of  these  are  in  high  inclination  orbits  where  they  sweep  p««t  each  oAer  at 
average  speed  of  10  kmAs.  About  1500  ^nt  rocket  stages,  inactive  payloadi;  and  a  few  active  pay- 
ids  account  for  most  of  this  mass.  These  objects  ate  currently  tracked  by  the  USAF  Space  r-nmmitnit 
are  about  4500  others  totaling  20  000  kg.  mostly  fragments  of  satellites  or  ofeer  orbitiag  hmdnme. 
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Recent  observations  indicate  a  total  mass  of  about  1000  kg  for  orbital  debris  with  diameters  of  1  as  or 
tmaiip.r  and  about  300  kg  for  orbital  debris  sizes  smaller  than  1  mm.  This  distribution  of  mass  makES  the 
orbital  debris  environment  more  hazardous  than  the  meteoroid  envirooment  in  most  spacecraft  appfica* 
ions  below  2000  km  altitude. 

The  detuis  model  presented  in  the  following  represents  an  extenaon  of  the  model  presented  by 
Kessler  et  al.^^  Rrst,  a  curve  fit  to  the  “current”  debris  environment  was  developed  based  on  the  best 
;xperimental  data  available.  This  was  then  coupled  with  additional  terms  which  represent  a  projection 
he  expected  environment  change  into  the  future.  The  year  1988  was  selected  as  tl^  oase  year  for  the 
'current”  environment.  The  applicable  data  sets  include  those  referenced  by  Kessler  et  aL,27  ^ 

malysis  of  panels  returned  from  the  Solar  Max  satellite,  data  from  the  MTT  ETS  telescopes,  and  USAF 
Space  Cnmnnand  cataloged  and  uncaialoged  data  sets.  In  addition,  updates  were  added  from  recent 
jptical  measurements  made  with  the  GEODSS  telescope  systems  and  from  biphase  radar  mcamreroens 
Tiadft  with  the  Arecibo  and  Goldstone  radars.  Quick-look  data  from  the  Long  Duration  Exposme  FacSag 
TDEF)  appears  to  be  consistent  with  these  data  sets,  but  analysis  was  not  sufficiently  compkie  to  allaw 
ncorporation  of  this  information.  The  model  is  valid  for  debris  diameters  of  1  fim  and  larger  (althonghh 
s  known  that  some  debris  smaller  than  1  pm  does  exist).  Measurements  made  with  the  Haystack  radar 
ifter  this  model  was  developed,  through  1992,  fit  well  within  the  uncertainty  bands  of  the  model, 
■iowever,  they  do  seem  to  incficaie  that  there  is  some  additional  debris  in  low  inclination  orUt  that  was 
lot  accounted  for,  and  that  the  size  of  the  material  in  the  3-  to  10-cm  size  range  was  overesthaatedby  a 
'actor  of  about  1.4.  These  data  ate  stiU  being  studied  so  that  no  changes  in  the  current  model  are 
;xpected  in  the  near  future. 

Model  development,  especially  prediction  of  future  trends  in  the  debris  environment,  is  difficult 
md  subject  to  substantial  uncotainty.  The  problems  may  be  grouped  into  two  categories  winch  must  be 
seated  somewhat  differently.  The  categories  are; 

(a)  Uncertainties  in  the  current  environment.  These  uncertainties  shift  the  flux  by  a  factor 
vhich  is  independent  of  time  (a  change  in  intercept  on  a  flux  versos  time  plot).  They  include  miceismn- 
ies  of  measurement,  statistical  limitations  of  the  data  sets,  and  voids  in  the  data  sets.  Le..  size  nd  alti- 
ude  ranges  where  no  measurements  have  been  made,  limitations  of  debris  shape  and  density  Monaa- 
ion,  etc. 

(b)  Uncertainties  related  to  trend  projection.  These  factors  alter  the  slope  of  a  flux  vtcas  lime 

)lot  They  derive  primarily  from  die  assumptions  which  must  be  made  to  predict  future  trends  in  hanun 
kctivities  (launch  rate  for  example,  tiriiich  historically  has  deviated  significantly  from  the  traffic  model 
vojections)  and  from  assumptitHis  needed  to  overcome  the  lechnkal  uncertainties  listed  abore:  The 
itate-of-tbe-art  is  such  that  understanding  is  still  lacking  in  several  important  areas  necessary  for  com- 
)lete  engineering  analysis  and  numerical  modeling  of  the  enviroomenL  Key  examples  include  mcoai- 
)lete  knowledge  of  satellite  and  rocket  body  fragmentation  mechanics  (e.g.,  fragment  size  and  velod^ 
listributions)  and  uncertainty  in  the  cause  and  intensity  of  many  fragmentation  events.  There  me  also 
mportant  limitations  in  modeling  capability  due  to  the  necessity  keeping  computation  time  and  roodtl 

:ompIexity  within  reasonable  limits.  These  all  have  important  influence  on  the  final  results. 

Before  discussing  these  issues  in  detail,  we  present  the  numerical  formalism  of  the  model  in  the 
oUowing  three  sections.  Then,  in  section  7.2.5,  the  uncertainties  and  assumptions  are  listed  in  detaB 
dong  with  a  quantitative  estimate  of  their  importance. 
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■".Wing  surt«»  ml  ihe  «|oc»y  dlaniwlico. 

>=?*  «  SfelMncSMlurtte  *■“' 

Sivenl^  the  foUowing  equation:  ‘“^*^^’*"“«»wacimty  was5l3arpMoriotis 


where 


s  flux,  iapacts  per  square  meter  of  smfice  per  year, 

d = orbital  (id)ris  diameter  in  cm,  (I(H  iSrfs  500) 
rs  date  (year). 

A  s  altitude  m  km  (A  ^  2000  km), 

S=  13-month  smoothed  solar  radio  flux  F-7  for  M  year  ,„  in4  r  -  -  . 

i = inchnatioB  in  degrees. 


H(d)  =S 

P(A,5)  =  p,(Wy(^,(A^l) 


Fi(d)=L22xlO~^d~^ 

=  8.1xl0"(rf+700)'® 

P — the  assumed  annual  growth  rate  of  mass  in  oitit = 0.05, 

^ and  the  estimated  growth  rate  of  fiagm^t  g  -  on?,  y-  a  oi  -n,  _ _ _ 

nsedfor2011andlaierdates.^^  ^ 

Pivi  =  (l+g)(*-t««)  for  r<  2011- 


forr<2011. 


^i(l)  =  (l+g)23(l4^')<'-2oii)  forri2011. 
fc(»)*l+p(M988). 
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The  inclination-dependent  function,  y/ii),  defines  the  relationship  between  the  flax  oa  a  spoce- 
aaft  in  an  orbit  of  inclination,  i,  and  the  flux  incident  on  a  spacecraft  in  the  current  population's  uwen^ 
nclination  of  about  50*.  Values  for  p(i)  are  as  follows; 


Inclination  i  (degrees! 

28.5  0.91 

30  0.92 

40  0.96 

50  1.02 

60  1.09 

70  1.26 

80  1.71 

90  1.37 

100  1.78 

120  1.18 


An  example  orbital  debris  flux  is  compared  with  the  meteoroid  flux  from  equation  (7-4)  in  figure 
-2  for  A  =  500  km,  r  =  1995,  i  =  28.5*,  and  5(t-l  yr)  =  97.0. 

7.2.3  Average  Shape  and  Mass  Density.  The  state  of  knowtedge  of  A»Jiris  <hap<»  anrf  Omaay  k 
ery  scant  Actual  diapes  are  irregular,  including  flat  plates,  rods,  hollow  structures,  and  crunpled 
letaL  As  size  decreases,  the  objects  tend  to  be  somewhat  less  irregular.  For  the  purposes  of  lUs  rerndd, 
le  objects  are  assumed  to  be  spherical,  witli  a  aze-dependent  mass  densi^  fimction  to  iippmiimtr 
lese  irregularities  and  the  probability  that  they  may  impact  with  any  tuientation. 

The  average  mass  density  tm  debris  0.62  cm  and  larger  is; 

p  =  2.8<H)-74  (g/cm3,  d  in  cm)  .  (7-8) 

or  debris  smallm  than  0.62  cm, 

p  =  4.0g/cm3  (d<  0.62  cm)  . 

hat  is,  for  small  objects  the  mean  density  should  be  assumed  to  be  a  amstant,  4.0  g/cm^.  imfrirrndpn* 
f  size.  Actual  shape  and  density  distributions  are  very  broad.  This  issue,  in  particular  deareity 

:  opposed  to  object  density,  is  addressed  further  in  section  7.2.5.3. 

7.2.4  Velocitv  and  Direction  Distribution.  Av<»rap<»rf  ftv<»r  all  aititiiHx  Hu. 

liliskm  velocity  distribution,  i.e..  the  number  of  impacts  with  velocities  between  v  and  vHh,rda^ta 
^lacecraft  with  orbital  inclination,  i,  and  considering  the  debris  is  in  a  locally  tmrrrontal  fdiire,  ispiven 
'  the  following  equations; 

/(v)  =  {2vV(rv2}{Cexp(-((v-Avoy(flvo))V^exp(-((»^o)®vjf}-i-flC(4vvo-v2)  .  (7-9) 

iiere  V  is  the  collision  velocity  in  km/s,  A  is  constant,  and  B.  C.  D,  E.  F.  G,  and  Vg  are  fuiKiia»  of 

i  orbital  inclination  of  the  spacecraft  The  values  fmr  these  constants  and  p»fatn«»ifry  are  as 
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A  =  Z5 


fO.5 

|■<60 

0.5-0.01(1-60) 

60<i<80 

lo.3 

i>80 

^_/0.0125 

i<100 

"  \0.012540.00125(i-100) 

/>100 

D*  13-0.01(1-30) 

£=035+0.005(1-30) 

|0.3+0.0008(i-50)* 

i<50 

£=<0.3-0.01(WO) 

50<i<80 

lo.o 

i>80 

nsj 

i<60 

G=<  18.7+0.0289(r-60)* 

60<i<80 

1250.0 

{>80 

»  =  1.0-0.0000757(1-60)2 

/7.25+0.015(i-30) 

i<60 

II 

i>60  . 

Wlicn/(v)  is  less  than  zero,  the  function  is  to  be  reset  equal  to  zero.  The  user  may  find  itcon- 
vement  to  numerically  normalize/(v)  so  that; 

■  an 

I  f{y)dv 

When  normalized  in  this  manner,/(v)  over  any  1  km/s  velocity  interval  becomes  the  fraction  of  ddns 
impacts  within  a  1  km/s  incremental  velocity  band.  The  function  is  fflustrated  in  figure  7-3.  Any 
velocity  moment  may  be  defined  as  ^  ^ 


V" 


JvV'(v)rfv  . 


(7-11) 
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Figure  7-3.  Nbnnalized  collision  velociQr  distribution  as  function  of  the  debris  veloci^ 
for  a  spacecraft  with  orbital  inclinations  of  28.5%  57*.  and  98*. 

Frequency  of  impact  from  a  given  direction  can  be  estimated  by  using  this  velocity  distributioo. 
Direction  of  impact  is  assumed  to  be  specified  by  tite  intersection  of  the  spacecrafi  velocity  vector  md 
another  circular  orbit.  That  is,  the  relative  velociQr  vectors  may  be  obtained  by  vector  addition  in  a 
plane  tangent  to  the  Earth’s  surface.  Since  a  spacecraft  velocity  of  7.7  km/s  was  used  to  calculate 
relative  velocity,  the  diiection  of  the  relative  velocity  vector  is  given  by  the  relationship: 


cos(±a)=-^  =  - 


(7-12) 


where  a  is  the  angle  between  the  impact  velocity  vector  and  the  spacecraft  velocity  vector  in  a 
coordinate  system  fixed  with  respect  to  the  space  vehicle,  v  is  the  impact  velocity,  and  is  die 
maximum  possible  velocity  difference  between  the  debris  and  the  spacecraft.  The  coordinate  system  is 
illustrated  in  flgure  7-4. 
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DEHNmONS: 

Plane  A  represents  a  siuface  of  the  ^ncecraft 
N,  is  the  unit  vector  normal  to  the  plane  A. 

X  is  the  direction  of  travel  of  the  spacecraft 

is  the  debris  velocity  relative  to  the  spacecraft 
O  is  the  tangent  plane  (horizontal)  to  the  spacecraft’s  orbit 
A  right-handed  coordinate  system  (positive  x.  y.  z)  is  defined  in 
x:  direction  of  spacecraft  travel 
y:  90*  from  x  and  in  plane  O  (port  direction) 
z:  Earth  vertical  (up) 

ANGLES: 


a  is  the  angle  between  x  and 

6  is  the  zenith  of  with  respect  to  the  Z  axis  in  thii'  reference  frame. 

^  is  the  azimuth  of  with  respect  to  the  spacecraft  direction  of  travel. 

Hgure7-4.  Orbital  debris  reference  frame. 
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Growth  Curvesl.  For  orbital  debris  sizes  larger  dian  10  cm  in  diameter,  the  envirooment  is  generally 
measured  by  ^und  radars.  Ute  most  extensive  measurements  are  made  by  the  USAF  Space  PfimmaiMi, 
which  also  maintains  acatalog  of  the  debris  population.  While  these  data  provide  an  adequate  description 
of  the  distributions  of  large  debris  with  re^ct  to  altitude  and  inclinatioa.  and  of  histotical  trends,  anafy- 
sis  of  GEODSS  optical  telesct^  data  has  shown  that  the  radars  are  detecting,  and  the  Space  Command 


is  cataloging,  less  than  half  of  the  population  in  this  size  range.  This  information  has  been  incorporated 
in  the  model  presented  here,  so  that  the  model  represents  the  current  environment  in  this  range  accu 
rately  to  within  the  range  1.5  to  0.5  times  the  flux,  Lc.,  the  •‘90-p2rcent  confidence”  upper  limit  flux 
equals  1.5  times  the  flux  from  equation  (7-7),  etc.  Table  7-1  summarizes  the  uncertainties  and  accuracy 
limitations  in  the  orbital  debris  flux  model. 


Me^urement  of  the  debris  flux  at  the  other  size  extreme,  sizes  smaller  than  0.05  cm,  is  made  by 
analysis  of  impact  craters  on  pieces  of  space  hardware  returned  from  orbit.  Meteoroid  impacts  are  dis¬ 
tinguished  from  debris  impacts  by  analysis  of  the  chemical  elements  retained  in  the  crater.  For  these 
sizes,  the  flux  has  only  been  measured  on  hardware  flown  at  about  500  km;  at  this  altitude  the  environ¬ 
ment  is  known  within  the  range  2  to  0.5  times  the  flux. 

Until  lecendy,  the  only  measurements  between  the  two  debris  size  extremes  was  a  limited  set  of 
opticaj  telescope  data  from  the  MIT  ETS  telescopes.  This  provided  an  indication  of  the  cumulative  flux 
for  objects  believed  to  be  2  cm  and  larger.  For  intermediate  sizes,  the  environment  was  estimated  by  a 
simple  straight  line  interpolation  on  a  log-log,  flux  versus  size  plot,  as  in  figure  7-2.  This  practice  was 
retained  for  the  current  model,  but  in  this  case  the  interpolation  is  conflrmed  by  recent  measurements  by 
Arecibo  and  Goldstone  radars  in  the  important  midrange  between  0.2  and  2  cm.  These  show  a  detection 
rate  which  is  consistent  with  the  current  model,  but  both  systematic  and  random  errors  in  these  meas¬ 
urements  leave  the  environment  uncertain  within  3  to  0.33  times  the  flux  for  these  sizes.  Measurements 
were  only  made  between  500  and  600  km  altitude.  Between  2  and  10  cm,  no  measurements  exist,  tat 
interpolation  of  the  data  from  either  side  and  modeling  this  region  also  yielded  an  estimated  uncertainty 
between  3  and  0.33  times  the  best  estimate. 


12.5.2  Trend  Projection  (Factors  Which  Alter  the  Slone  of  Hux  Growth  rurvesV  As  is  the  case 
with  any  analysis  or  model,  the  results  hold  only  so  long  as  the  underlying  assumptions  remain  valid. 
The  following  are  the  key  assumptions  upon  which  the  model  rests: 


(1)  It  is  assumed  that  the  rate  of  accumulation  of  mass  in  LEO  Ls  constant,  with  the  annual 
increase  equal  to  5  percent  of  the  amount  accumulated  by  1988.  This  matches  the  historical  trend  over 
the  last  few  decades. 


(2)  The  rdadve  use  of  different  orbits  is  assumed  to  remain  constant  For  example,  the 
history  of  launches  by  the  USSR  has  been  such  that  80  percent  of  their  payloads  re-enter  within  2  years 
of  launch.  These  do  not  contribute  significantly  to  the  debris  environment.  If  this  practice  changed  with 
increased  use  of  higher,  longer  life  orbits,  the  population  of  objects  in  orbit  would  grow  at  a  proportion¬ 
ally  increased  rate. 

(3)  It  is  assumed  that  the  efforts  to  minimize  fragmentation  of  satellites  in  orbit  will  con¬ 
tinue  such  that  fragmentation  events  will  continue  at  the  rate  of  only  one  per  year  in  LEO  (29  percent  of 
past  practice).  In  the  last  decade,  intentional  (or  apparently  intentional)  fragmentation  of  satellites 
accounted  for  about  71  percent  of  the  known  fragmentation  events.  Apparently,  recent  publicity  and 
increased  awareness  of  the  hazards  associated  with  orbital  debris  has  generated  policy  shifts  among  the 
space-faring  nadons.  No  intendonal  fragmentation  events  above  300  km  have  been  observed  in  the  last  5 
years.  Indeed,  for  a  period  of  about  2  years,  there  were  no  significant  events  at  all.  However,  from 
October  1990  through  Fdmiary  1991  there  have  been  14  unintentional  breakup  events,  about  half  of 
which  must  be  considered  significant  This  is  roughly  equivalent  to  the  breakup  rate  observed  in  the 
early  1980’s. 
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Table  7-1.  Uncertainties  and  accuracy  limitaticns. 


Uncetlainty  in  Cui 

Tent  Environment  (Intercept  Shift) 

Treatment 

Percent  ConfidemV* 

Notes 

Flux  Meastuements 
(d>10cm) 

BestEst 

lJto0.Sx(1988FIux) 

Flux  Measurements 
(0.0S  cm  <  d  <  10  cm) 

BestEst 

3  to0.33x(198S  Flux)  in  this 
size  range 

Due  to  statistical  and  measuaemcnt  limi¬ 
tations  in  portions  of  range,  data  missing 
(interpolatioo  used)  in  rest  of  ranee. 

Flux  Measurements 
(d<  0.05  cm) 

BestEst 

2  to0.5x(1988Flux)  in  this 
size  range 

Altitude  Distribution 
(d<  10  cm) 

BestEst 

S  to  0.2x( 1988  Flux)  per  200 
km  away  from  500  km 

Due  toififiiculty  in  detenmnn^  flux  in 
highly  elliptical  orbits 

Altitude  Distribution 
(d  2 10  cm) 

Smoothed 

BestEst 

2lo0.5x(1988Flux) 

See  Figure  7-5 

Unceminty  is  somewhat  worse  in  800 
and  1000  im  regions 

Debris  Density 
(d<  1cm) 

Simplified 

BestEst 

0.10  @  1.8/03918 
0.10943/OJ9  8.9 

Estimmed  typical  lieav/*  (fistiibutioa. 
Insufficient  data  to  devdop  a  true 
unccftnty  limit  estimaie. 

Debris  Density 
(d>  1cm) 

Simplified 

BestEst 

2  to  0.5  Mean  Density 

Mean  values  are  faiily  weU  defined  but 
number  vs.  density  detribation  is  broad 

Debris  Shape 

Simplified 

Nonconserv. 

Spherical  sbape  is  assum 'xL  actual  debris 
willbeciegular. 

Velocity  Distribution, 
Fraction  <  S  kmA 

BestEst 

0.5  to  3x(Slow  Baction) 

Distiibacion  of  orbii  inclinations  could  be 
in  error  or  change  with  time. 

Uncertainty  in  Trei 

Ml  Projection  (Slope  Shift) 

Treatment 

*<90  FerccMCoondcncc^ 

Notes 

Launch  Rate 

Orbit  Use  Profile 

BestEst 

BestEst 

} 

p  a  0.04  to  “comp  p"  «  0. 1 
qs0to2p 

**cooipp’*  implies  g2  =  (lip)!*’*^^). 
Worst  case  assunies  combined  effect  01 
increased  traffic  and  iocseaKd  use  of 
LEO  above  400  km 

Fragmentation  Rate 

Bagmentation 

Mechanics 

BestEst 

BestEst 

} 

qs0to0.10 

Assuming  no  changes  in  projected 
launch  rate  and  oihit  nse  profile. 

Statistical  Variation  of 
Fragmentation 

BestEst 

03  to  1.5x(CniTentFlux) 

Solar  Activity 

BestEst 

Substitute  “Max”  and  “Min” 
s  values  fona  Table  3-1 

Modd  tends  to  ovcrestimaie  variation 
with  solar  cycle  io  these  woirid  be 
extieaie  limits 

“LocaT  Fragmentation 
Events 

Nonconserv. 

44x(1988Rnx)  fori  Year 

Difficult  K>  assess,  depends  strongly 
on  type  of  event  and  pnnimity  to 
statian  orbit  Ignored  in  canent 

modcL 

(4)  It  L  assumed  that  the  debris  size  distribution  is  independent  of  altitude:  One  would 
expect  small  debris  in  circular  orbits  to  decay  faster  than  large  debris,  implying  an  altitude-dependent 
distribution.  However,  consideration  of  the  population  of  small  fragments  in  elliptical  orbits,  assuming 
trend  similar  to  the  one  for  large  fragments,  leads  to  a  dependent  distribution  with  the  opposite  trend. 
Therefore,  pending  further  measurements  and  research,  an  intermediate  assumption  been  made. 
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(5)  It  is  assumed  that  future  solar  cycles  will  affect  debris  orbit  decay  uniformly  for  all 
sizes.  The  unceiiainties  associated  with  these  factors  differ  from  those  discussed  above  in  that  they  alter 
the  growth  term.s  in  equation  (7-7)  rather  than  entering  as  multipliers  of  the  total  flux.  Thus,  the 
variations  are  expressed  as  variations  of  the  p,  q,  and  S  terms. 

The  fust  two  of  these  assumptions  relate  to  the  predicted  future  accumulation  of  objects  in  LEO. 
This  is  a  key  determinant  of  the  expected  debris  growth  rate— the  value  of  p.  The  combination  of  a 
decreasing  launch  rate  for  the  United  States  with  an  increasing  rate  for  the  rest  of  the  world  has  led  to  the 
relatively  constant  historical  trend.  It  is  not  clear,  however,  that  this  trend  will  continue.  Expected 
launch  rates  are  subjea  to  political  and  economic  influences  which  may  change  unpredictably,  and  many 
new  countries  are  becoming  involved  in  space  activities.  Current  traffic  models  extended  to  the  year 
2010  predici  that  future  growth  will  be  ata  compounded  rate  between  5  and  10  percent  per  year.  The 
lower  limit  of  these  models  corresponds  to  a  value  of  p  =  O.OS  compounded  annually.  These  models 
represent  the  projects  that  are  planned,  and  since  some  projects  are  either  canceled  or  postponed,  the 
actual  rate  has  always  fallen  below  the  traffic  models.  Therefore  the  baselined  rate  is  a  constant  (not 
compounded)  p  =  0.05.  A  substantially  lower  rate  would  only  be  expected  in  the  case  of  a  worldwide 
economic  depression  or  similar  event.  A  higher  rate  is  possible,  especially  if  an  increa<%d  launch  rate  is 
coupled  with  an  increased  use  of  higher,  long-life  orbits.  A  compounded  10  percent  per  year  increase,  or 
p  =  0.1,  represents  the  “90-percent  confidence  level”  upper  limit. 

Tne  value  of  q,  which  represents  the  expected  growth  rate  of  small  debris,  primarily  depeiHls  on 
the  frequency  of  expected  satellite  breakup — assumption  (3).  Breakups  may  be  intentional,  or  they  may 
result  from  accidents  and  random  collisions.  Thus,  Ae  breakup  rate  is  partially  controllable,  partially  not 
The  value  selected  for  this  model,  q  =  0.02.  assumes  no  intentional  breakups,  and  an  accidental  breakup 
rate  of  1  per  year.  The  range  of  possible  q  values  is  from  0.0  (random  collisions  are  not  important  and 
improved  precautions  lower  the  accidental  rate  below  historical  values)  to  0.06  (both  accidental  and 
intentional  fragmentation  rates  match  the  1980  to  1987  rates).  Technically,  negative  values  of  q  are 
possible  if  all  fragmentation  events  are  prevented,  but  this  is  not  considered  a  credible  possibility.  At  the 
other  extreme,  if  fragmentation  rates  rise  above  the  1  per  year  and  the  rate  of  accumulation  of  mass  in 
orbit  increases,  on-orbit  collisions  become  important,  <  ad  ^  will  approach  2p. 

Unintentional  fragmentation  events  can  result  from  either  explosions  or  collisions  between 
objects.  The  first  of  these  represents  a  simple  linear  source:  the  population  growth  is  directly  related  to 
the  fragmentation  rate.  The  second  is  exponential  in  nature,  rince  &e  number  of  fragmentation  events  is 
a  function  of  the  square  of  the  population.  Thus,  the  coupling  between  p  and  q  uncertainty  limits  noted 
above.  Modeling  these  processes  indicates  that  the  first  procesa  will  be  dominant  until  about  the  year 
2010.  About  that  time,  under  the  assumptions  of  the  current  model,  the  second  process  becomes 
significant,  and  the  small  particle  population  will  grow  at  an  increasing  rate.  If  the  assumptions  of  the 
current  model  hold  and  current  practices  continue,  q  is  expected  to  increase  to  0.04  beginning  about  the 
year  2010. 

Uncertainties  in  the  mechanics  of  fragmentation  events.  i.e.,  fragment  size,  number,  and  velocity 
distribution,  impact  the  capability  to  model  and  analyze  the  debris  environment  Fragmentation  events 
have  multiple  possible  causes  and  may  vary  widely  in  intensity.  Direct  data  from  simulations  of  these 
processes  are  quite  limited.  In  the  analvris  supporting  the  current  moiiel.  these  limitations  were  over¬ 
come  by  tuning  the  fragmentation  moc :  >  so  that  the  historical  fragmentation  record  reproduced  the 
current  environmenL  This  approach  provides  the  “best  estimate”  fragmentation  models  for  the  analyris, 
but  it  is  very  limiting  in  the  sense  that  there  is  no  independent  check  of  the  analysis.  Since  there  are  only 
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recent  measurements  of  the  small  fragment  environment,  there  is  no  second  point  to  check  the  analysis 
against  The  associated  uncertainty  has  been  included  within  the  2p  factor  described  above. 

Assumption  (4)  is  necessary  because  there  were  no  measurements  of  debris  smaller  than  2  cm  at 
other  than  in  the  500  to  600  km  altitude  range  when  tins  model  was  developed,  and  because  poor 
knowledge  of  the  fragment  velocity  distributions  (and  computer  run-time  limitations)  makes  it 
impractical  to  model  small  fragments  in  highly  elliptical  orbits.  The  recent  Haystack  measurements 
indicate  this  assumption  is  good  up  to  1000  km,  but  falls  apart  for  higher  altitudes.  No  data  are  available 
at  1500  kin  or  above.  Additional  study  is  needed  in  this  area.  Analysis  indicates  that  the  actual  flux  of 
small  material  could  trend  either  above  or  below  the  large  object  distribution  with  altitude,  depending 
upon  the  assumptions  used.  The  model  presented  here  assumes  that  the  distribution  with  altitude  for  the 
small  material  matches  the  catalog  distribution.  The  actual  amount  that  these  fluxes  differ  could  be  as 
high  as  a  factor  from  5  to  0.2  for  every  200  km  away  from  the  500  km  altitude.  The  distribution  with 
respect  to  altitude  is  also  assumed  to  be  smooth.  Actually,  the  USAF  Space  Command  data  (sizes  larger 
than  10  cm)  give  fluxes  at  800  and  1000  km,  which  approach  the  level  predicted  by  the  recommended 
flux  model,  as  shown  in  figure  7-5,  not  allowing  for  the  correction  factor  from  the  GEODSS  telescope 
study. 

With  respect  to  assumption  (5),  the  possible  variance  of  solar  activity  can  be  well  defined  based 
on  the  historical  rec  .rd.  A  high  solar  cycle  will  increase  the  depletion  rate  for  debris  in  low-altitude 
orbits,  compared  *0  mean  or  low  cycle,  but  prediction  of  the  solar  activity  level  beyond  about  1  year  is 
highly  uncertai:  '.able  3-1  provides  profiles  of  maximum  and  minimum  solar  activity,  and  these  values 
may  be  used  in  eolation  (7-7)  to  estimate  the  range  of  variation. 

An  impof'-'ti  rj)  term  factor  mat  is  not  included  in  the  model  and  thus  contributes  an 
additional  uncermirny  is  i.  flux  arising  from  the  intentional  or  inadvertent  fragmentation  of  a  satellite  in 
an  orbit  at  or  near  tiic  nr-pot  ed  space  vehicle  operating  altitude.  In  the  region  of  the  breakup,  an 
enhanced  flux  may  be  apparent  for  a  considerable  period  of  time,  depending  upon  the  altitude  of  the 
breakup  and  the  size  and  velocity  distribution  of  the  debris.  Analysis  and  modeling  of  various  scenarios 
indicates  that  such  an  event  would  probably  cause  increases  in  the  flux  environment  by  factors  of  a  few 
tens  of  a  percent  for  a  year  or  more,  although  a  factor  of  4  may  be  possible  as  an  extreme  worst  case. 

7.2.5.3  Uncenainlies  in  Direction.  Velocity  Distribution,  and  rtendry  Th<.  fart  th^t  nrhs.oi 
debris  objects  are  not  in  exactly  circular  orbits  will  introduce  a  small  error  in  direction.  As  a  result  of  the 
currentiy  small  eccentricities  of  debris  orbits,  the  actual  directions  of  impacts  are  within  1*  for  most 
velocities  derived  from  section  7.2.4.  For  velocities  less  than  2  km/s,  the  uncertainty  is  much  larger,  with 
a  significant  fraction  being  more  than  20*  from  the  direction  derived  from  section  7.2.4.  These  errors  in 
direction  can  be  in  the  local  horizontal  plane,  or  they  can  appear  as  direction  errors  above  or  below  this 
plane. 


Uncertainty  in  the  distribution  of  debris  orbit  inclinations  leads  to  an  uncertainty  in  the  velocity 
distribution  which  can  affect  penetration  analysis.  Since  inclinations  are  only  known  for  large  (cata¬ 
loged)  debris,  the  sm^  fragments  may  have  a  different  distribution,  or  the  distribution  may  change  with 
time  as  a  result  of  orbit  selection  and  fragmentation  events.  These  considerations  imply  the  slow  fraction 
of  the  population.  Le.,  the  fraction  of  debris  objects  with  a  relative  speed  less  than  5  km/s  with  respect  to 
a  vehicle,  could  shift  or  be  in  error  by  a  factor  from  0.5  to  3. 
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Attitude  (km) 

Figure  7-5.  Comparison  of  model  flux.  Fc  with  catalog  flux  not  corrected  for  GEODSS  results. 

The  expression  for  debris  density  given  in  section  7.2.3  for  objects  larger  than  0.5  cm  has  been 
verified  by  direct  measurements  of  actual  objects,  studies  of  orbit  decay,  and  fragmentation  experiments. 
Thus,  it  is  believed  to  be  a  good  representation  of  the  mean  density  of  the  debris  population  (within  a 
factor  of  0.5  to  2),  especially  for  sizes  above  30  cm.  However,  it  represents  the  mean  of  a  broad  distri¬ 
bution;  Use  density  of  individual  objects  can  vary  widely. 
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For  small  fragments,  the  density  issue  is  more  difficult  because  information  is  extremely  scant 
To  lUustrate  the  problem,  consider  the  following  materials  density  profile  based  on  a  summary  review  of 
space  shuttle  materials  usage  (neglecting  tiles):  ^ 


Estimated  Vidume 
Fraction 

Specific  Gravity 

Representative 

Materials 

0.65 

2.8 

Aluminum,  Glass 

0.15 

1.8 

Epoxy-Glass,  Rubber 

0.05 

4.5 

Utanium 

0.15 

7.8  to  8.9 

Copper,  Steel 

One  might  expect  that  this  would  be  typical  of  many  spacecraft  but  it  may  underestimate  the 
taction  of  dense  maten^als  for  several  reasons.  First  many  objects  involved  in  fragmentation  events  arc 
beheved  to  have  had  a  higher  fraction  of  dense  materials  used  in  their  consuuction.  The  Delta  second 
stages,  for  example,  contain  about  70-percent  steel,  20-percent  aluminum,  and  10-pcrcent  titanium. 
Second,  low-density  objects  are  more  affected  by  drag  and,  thus,  decay  from  orbit  more  rapidly.  How¬ 
ever,  fragment  sh^  is  also  an  important  determinant  of  effective  density.  The  thirifn<»gg  of  the  Delta 
w^ls  is  between  0.2  and  0.5  cm,  so  fragments  larger  than  this  would  be  irregular  in  shape  with  an 
effecuve  density  less  than  that  of  steel.  Since  definitive  studies  have  not  been  done,  uncertainty  bounds 
c^ot  be  defined  at  this  time.  For  study  purposes  the  recommended  “heavy”  distribution  is- 10  percent 
(by  voluine)  at  1.8. 50  percent  at  2.8, 10  percent  at  4.5,  and  30  percent  at  8.9  g/cm3.  TTiis  gives  a  mean 
density  of  4.7  g/cna^  for  solid  spheres. 


73  Evaluation  of  Directionality  Effects 


The  parameter  k  which  appears  in  equation  (7-1)  is  defined  as  the  ratio  of  the  flux  against  an  ori¬ 
ented  surface  to  the  flux  against  a  randomly  tumbUng  surface.  Introduction  of  the  k  factor  allows  appli¬ 
cation  of  the  flux  (equation  (7-7))  to  evaluate  the  expected  number  and  probability  of  impacts  on 
surfaces  flying  iwth  fixed  orientation.  Evaluation  of  penetration  probabUity  should  be  by  numerical 
techniques  which  account  for  the  directional  dependence  of  both  the  penetration  equations  and  the 
^teoroid  ai^  debris  fluxes,  although  the  k  factor  may  be  useful  for  quick  approximations.  The  value  of 
k  ^  tteorettcally  range  from  0  to  4;  a  value  of  4  can  only  be  achieved  when  a  surface  normal  vector  is 
onented  in  the  direction  of  a  monodirectional  flux.  It  depends  on  the  orientation  of  the  surface  with 
i^ct  to  the  Earth  vertical  and  the  spacecraft  velocity  vector.  If  the  surface  is  randomly  oriented,  then 


must  be  taken  in  evaluating  k  factors  and  other  directional  effects  because  of  the  complex 
directional  nature  of  the  meteoroid  and  debris  fluxes.  Unlike  most  fluxes  with  which  the  engineer  and 
physicist  deals,  ^  meteoroid  and  debris  fluxes  do  not  have  a  unique  direction  associated  with  them  at 
My  given  point  in  space.  Meteoroids  are  equally  likely  to  appear  from  any  direction  (except  where  the 
Earth  provides  shielding)  m  a  reference  frame  fixed  with  respect  to  the  Earth;  they  tend  to  appear  from 
the  ram  direction  on  an  orbiting  satellite.  The  relative  velocity  with  respect  to  a  randomly  tumbling 
^^raft  IS  about  20  km/s.  The  approach  for  evaluating  k  for  meteoroids  wUl  be  similar  to  the  approach 
for  (tebns  which  is  presented  below.  As  an  illustration  of  the  expected  effect,  the  rara-to-lee  ratio  of  the 
number  of  impacts  was  found  to  be  about  7  to  1  in  preUminary  analysis  of  LDEF  data  for  constant  par¬ 
ticle  size,  or  about  18  to  1  for  constant  crater  size.2«  It  is  expected  that  meteoroids  were  dominant  inthis 

C3SC. 
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For  orbital  debris  the  directionality  in  a  reference  frame  fixed  with  respect  to  the  space  vehicle  is 
defined  by  combining  equations  (7-10)  and  (7-12).  That  is,  the  diiectionali^  can  be  written  as  a  function 
of  velocity  alone,  so  differential  pseudo- vector  fluxes  can  be  defined  in  terms  of  the  velocity 
distributions  such  that: 


-rfF+(v)  =  F{(r(v)/2)</v+l 

and  (7-13) 

-dFJv)  =F{(f(vV2)dvJ, 

where /'(v)  is  defined  by  equation  (7-10)  and  the and  -  subscripts  are  associated  with  the  +aand-a 
angles  of  equation  (7-12),  i.e.,  the  left  and  right  lobes  of  the  “butterfly-shaped”  debris  distribution  (sym¬ 
metric  about  the  direction  of  flight).  Our  sign  convention  is  such  that  dF  is  positive  in  the  minus  dv 
direction  (fig.  7-4).  By  solving  the  problem  of  the  flux  against  a  ^)here,  it  can  be  shown  that 
F  =  4  F,  =  where  F^  is  defined  by  equation  (7-7). 

To  find  the  expected  rate  of  impacts  on  a  surface.  A,  flying  with  fixed  orientation  one  must  solve 
the  following; 


F(A)  =  J  J  dF^{v)  N^da  +  J  J  dFJy)  N„da  ,  (7-14) 

V  A  V  A 

where  is  the  outward  unit  vector  normal  to  the  surface  element,  da.  IMPORTANT:  The  limits  on 
the  surface  integrals  must  be  such  that  all  of  the  surface  where  the  dot  product  is  positive  is  included, 
and  portions  where  it  is  negative  are  excluded.  (A  negative  dot  product  corresponds  to  flux  leaving  the 
surface.)  The  k  factor  for  the  surface  A  is  simply: 

K  F(A) 

AF,  • 

Figures  7-6a  and  b  and  7-7  illustrate  k  factors  for  a  flat  plate  and  a  right  circular  cylinder  at 
various  orientations  and  inclinations. 


AZIMUTH  A»IGt£.  (DEGREES) 


Hguie  7-6i>.  K  factor  for  a  single-sided  flat  plate  for  28.5*,  57*.  and  98*  inclination  orbits  (6=  90*X 
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Vra.  MAGNETIC  FIELD 


As  di%ussed  in  previous  sections,  the  Earth  has  a  magnetic  field  which  greatly  affects  the 
environment  in  the  LEO  region.  The  magnetic  field  traps  charged  particles  (discussed  in  sections  IV  and 
V)  and  deflects  low-energy  cosmic  rays.  The  field  is  basically  a  dipole  field,  but  the  magnetic  axis  is  not 
coincident  with  the  rotation  axis.  The  south  (north  seeking)  pole  of  the  magnetic  field  is  rotated  1 1  7* 
^ong  the  69*  west  longitude  tine  from  the  rotational  north  pole.  The  strength  of  the  dipole  results  in  a 
Ireld  strergth  at  the  Earth’s  surface  of  approximately  0.3  gauss  (G)  at  the  equator  and  0.6  G  at  the  poles. 
The  magnetic  field  direction  at  the  equator  is  horizontal  pointing  north.  At  the  magnetic  south  (north 
seeking)  pole  it  points  down  into  the  Earth.  Beyond  an  altitude  of  approximately  2000  km  strong 
currents  in  the  magnetosphere  cause  deviations  from  the  near-Earth  field. 

The  particle  L-shell  is  most  simply  defined  for  a  dipole  field  like  the  Earth’s  by  the  eouation  for  a 
magnetic  field  line 

r  =  Z,/?£Cos2  A  .  (8.1) 

where  A  is  the  magnetic  latitude,  r  the  radial  distance  to  the  frshell,  and  Re  the  radius  of  the  Earth.  Thus 
the  L  value  of  a  field  line  is  its  distance  in  the  magnetic  equatorial  plane  from  the  center  of  the  Earth 
expressed  in  Earth  radii.  The  L-shell  is  a  surface  of  revolution  of  this  line  about  the  dipoie  axis. 

Unlike  most  Earth/Sun  parameters  which  remain  constant  for  hundreds  of  years,  the  magnetic 
field  strength  changes  significantly  on  the  scale  of  5  to  10  years.  At  present,  the  magnetic  field  strength 
IS  decreasmg  at  approximately  0.1  percent  per  year.  Thus,  it  is  important  to  use  a  recent  version  of  the 
model  to  correctly  determine  the  magnetic  field.  Note,  however,  that  as  mentioned  in  section  V  of  this 
document,  the  current  magnetic  field  (IGRF,  1991)  should  not  be  used  to  calculate  trapped  radiation 
with  the  APS  and  AES  radiation  models  because  they  do  not  correctly  account  for  the  effect  of  the 
atmosphere  at  orbital  altitudes. 

In  addition  to  the  general  decrease  in  the  magnetic  field  strength,  geomagnetic  storms  caused  bv 
solar  acnvity  can  change  the  magnetic  field  by  hundreds  of  gamma  ( 1  gamma  (T)  is  10-5  G). 


The  magnetic  field  can  be  accurately  represented  using  a  spherical  harmonic  expansioa  The 
expansion  m  terms  of  the  Earth’s  rotational  coordinate  system  (i.e.,  normal  latitude,  longitude)  requires 
more  terms^m  a  simple  dipole  because  the  axes  of  the  coordinate  system  are  not  coincident  with  the 
axes  of  the  field.  However,  this  coordinate  system  is  so  familiar  that  it  is  the  usual  choice.  The  specifi¬ 
cation  IS  given  in  terms  of  magnetic  potential,  U,  and  the  field  is  derived  from  the  potential  bv  taking  the 
negative  gradient  (which  is  the  vector  derivative):  ^  ® 


In  particular; 


fl  =  -VU  . 

=  sin  e)]flUldp  . 

B^^dUtdr  . 


(8-2) 

(8-3) 

(8-4) 


(8-5) 
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amer  potential:  r.  d,  and  ^  are  the  spherical  coordinates;  and  down  means  toward  the 

Given  that  repreantation,  the  potential  is  expressed  as  Schmidt-normalized  Legendre 
polynomials  with  the  following  form:  senare 


where 


^  ,li  Jo  cos  (m4,  )+h:  sin  (m  0 )1  . 


d= co-latitude  (i.e.,  measured  from  the  pole) 

0  =  east  longitude 

Re = Earth  radius  =  (637 1.2  km) 

r  =  radial  distance 

~  Schmidt-normalized  associated  Legendre  polynomials 
g”  and  h"  =  Schmidt-normalized  coefficients . 

^  Legendre  polynomials  are  related  to  the  normal  associated  Ugendre 


polynomials  as  follows 


where 


=  [«  m(n-m)y(a+m) !]  /V  , 


Rnm  =  regular  associated  Legendre  polynomial 

_  2  for  m?sO| 

"  1  for /n  =  ol  ' 

Given  to  formtom,  the  coefficients  g  and  /r  can  be  stated.  Note  that  the  h  coefficient  for  «= 0  is  not 

**  coefficients  are  given  in  units  of  gamma  (f)  which  is 

idin^rfiSf  the  value  of  ^h  coefficient  depends  on  the  radius  of  the  Ea4.  so  tte  stated 
ra*us  (6371  2  km)  must  be  used  with  these  coefficients.  The  coefficients,  according  to  the  current 
Intemanonal  Gee .aagnetic  Reference  Field  (IGRF,  1991),  are  listed  in  table  8-1. 

Since  the  polynomials  />,"  arc  of  order  unity,  the  coefficients  in  table  8-2  terms  will  represent  the 

^  ®  the  worst  case.  The  complete  model  has  been  calculated  to  degree 

order  10  which  is  more  accurate.  Computer  programs  for  determining  the  field  values  are  avtoble. 
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IX.  THERMAL  ENVIRONMENT 


9.1  General  Discussion 

solar  ral'uSn  fdescri^bSl"^^^  the  incoming 

Spacecraft  auftacS  “min'd  ^  ““  “■*  <>f  ^ace  a,  3  I' 

energy  produced  iniemally  wiuiin  the  veMctel  similar 

itself.  If  one  considers  the  Earth  anTiSlll  ns  'Sjf'.l'’"''  f®  >l«  Eanh 

the  incoming  solar  energy  and  OLR  enerpv  arp  pccpm*®  n  •  k  ^v^ges  over  long  lime  penods, 
nearly  in  radiative  equilibrium  with  the  Sun.  H^ever  ti  is  not’SSn?^  Earth/atmosphere  is  very 
and  there  are  important  variations  with  resnect  to  inrai  fSm«  *  *"  bailee  everywhere  on  the  globe. 

A  space  vehicle’s  motion  with  respect  to  ^  ^arth  rpcnil^^’  geography,  and  atmospheric  conditions, 
fuil  global  thermal  profile;  si  it  SeS  vari!,Sc  «  r  -I-  o^y  a  “swath”  across  the 

time  constants  of  the  hardware  systems  Thus  the  •*"  thermal 

dependent  upon  the  orbSmS™ 
tailored  to  allow  selection  of  a  statiSfeir^^t  de^?S 

for  most  common  LEOs.  '  esenpuon  or  this  envuonment  and  its  variations 

c  ^9lflr  Consiflni*  The  direct  solar  flux  is  the  sreatest  vinrrp  nf  hontirtn  f 

craft.  The  mean  value  of  this  solar  flux  ai  mpan  Ltk  o  ®  a-  •  for  most  space- 

Specifically,  the  solar  constant  is  defined  as  the  radiation^lLf  constant." 

the  line  from  the  Sun.  per  unit  time  outside  the  atmncnhPivf^  ^  **"**  surface  norma]  to 

Earth-Sun  distance).  In  aeSy  astronomical  unit  (the  mean 

spacecraft  is  not  quite  constant.^Tliere  are  two^farmrfinn  impinges  on  an  E^h-orbiting 
radiant  energy  that  is  emitted  bv  the  Sun  is  known  ''^®^*Ety.  First,  the  amount  of 

and  differs  frL  cycL  to  cy^^^^^^  slighUy  throughout  the  1 1-yr  solar  cycle 

mated  to  be  less  than  half  a  percent  Second  the  siiohOv  is  still  ^ing  studied,  but  it  is  esti- 

also  results  in  a  variation  of  approximately  ±3  4  nercen?  SiSS*  the  Earth  about  the  Sun 

solstice  when  the  Earth-Sun  distance  is  a^inimi^  the  maximum  occumng  at  the  winter 
137 1±10  W/m2.  The  ±10  W/m^  allows  for  the  slight  snia^  con^t  defined  in  this  document  is 

the  measurement  uncertainty.  The  variaUon  with  respect  ^rEartiTslJrH^I^*’  "'T  ‘™P®^‘*y- 

may  be  calculated  by  standid  methods  denen^niTn^n  ^stance  must  be  added.  It 

user.  Extreme  values  are  listed  in  table  9-L  ^  interest  to  the 

spa«cn.fl  is  sLift  Alb^  ,5lL  vT™,2  “iSS  ““'J* 

account  for  this  effect.  eapeciaUy  neXT^^tatTfL 

approximately  ihe  same  spectral  shaoe  thp  ^iin*c  .  P^^graph  A2.2.)  Albedo  radiation  has 
a  characteristic  temperature  of  5762  K  As  used  in  this  dornm’  a  blackbody  with 

specirun.  albedo.  Albedo  values  b..e  a  widS  Mae  Jtd  “  ««=  “lal  soli 

bandwidths  were  considered.  ^  g'^a'er  if  smaller 
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Albedo  is  highly  variable  across  the  globe  and  is  dependent  on  the  distribution  of  reflective 
properties  of  the  surface  and  the  amount  and  type  of  cloud  cover.  Albedo  is  generally  highest  over 
cloudy  regions  such  as  the  intertropical  convergence  zone,  deserts,  and  ice-  and  snow-covered  areas. 
If  the  sky  is  clear,  the  albedo  over  ocean  areas  in  generally  low.  Albedo  also  increases  as  the  solar 
zenith  angle  increases.  Because  of  the  snow  and  ice  cover,  decreasing  solar  elevation  angle,  and 
increasing  cloud  coverage,  albedo  tends  to  increase  slighUy  with  latitude  if  viewed  on  a  large  scale. 

1-3  Outgoing  Longwave  Radiation.  The  third  rnmpnn^nt  nf  thermal 

environment  is  the  OLR  emitted  by  the  harth  itself.  This  Earth-emitted  thermal  radiation  hag  a 
spectrum  of  a  blackbody  with  a  characteristic  temperature  of  288  K.  OLR  is  not  constant  over  the 
globe,  but  the  localized  variations  are  much  less  severe  than  for  albedo.  OLR  is  principaUy  influenced 
by  the  temperature  of  the  Earth's  smface  and  the  amount  of  cloud  cover.  A  warmer  region  of  the 
E^’s  surface  will  emit  more  radiation  than  a  colder  area.  On  a  large  scale,  highest  values  of  OLR 
will  occur  in  tropical  regions  (as  these  are  the  regions  of  the  globe  receiving  tlie  maximum  solar 
heating)  and  will  decrease  with  latitude.  Increasing  cloud  cover  tends  to  lower  OLR  by  absorbing 
up-welling  radiation  from  the  Earth’s  surface.  The  solar  elevatiort  angle  may  also  affect  OLR 
because  ot  its  influence  on  the  temperature  of  the  surface  and  lower  atmosphere.  Thus,  both  diurnal 
and  seasonal  variations  may  be  detected.  This  diurnal  variation  is  small  over  the  oceans,  but  for  local 
regions  over  the  continents,  e.specially  desert  areas,  it  can  amount  to  about  20  pe.  cenu 

9.1.4  The  Earth  RadiatioaBudget  Experiment.  The  data  used  to  deflne  the  thermal  environ¬ 
ment  were  collected  by  the  Earth  Radiation  Budget  Experiment  (ERBE).  ERBE  began  in  the  1980’s 
and  is  still  ongoing.  ERBE  is  a  multisatellite  experiment  which  has  as  its  primary  objective  global 
data  collection  of  such  Earth  radiation  budget  parameters  as  incident  sunlight,  reflected  sunlight 
(albedo),  and  OLR.  This  experiment  was  selected  oecause  of  its  thorough  coverage  and  its  up-to- 
date  information.  The  experiment  consisted  of  three  satellites,  the  low-inclination  Earth  Radiation 
Budget  Satellite  (ERBS)  and  two  NOAA  Sun-synchronous  satellites.  The  data  used  here  are  from 
the  active  cavity,  flat  plate  radiometers  in  the  fixed,  nonsca.ining,  wide  field  of  view  (FOV)  mode. 

This  type  of  instrument  was  chosen  to  most  closely  approximate  the  albedo  and  OLR  vaiiations  that 
a  spacecraft  surface  would  encounter.  The  data  are  available  in  three  separate  sets,  daily  averaged 
values  (S-4).  hourly  averaged  values  (S-IO),  and  the  direct  16-s  instrument  measurements  along 
the  ERBS  or  NOAA  satellite  trajectory  (S-7).  The  S-4  and  S-10  data  products  were  evaluated  but 
were  deemed  inappropriate  for  this  applicahon;  the  averaging  times  are  much  too  long  compared  to 
the  thermal  time  constant  of  most  satellite  systems.  Therefore,  the  following  design  criteria  are 
based  on  28  files,  representing  1  mo  of  16-s  data  each,  of  S-7  data.  The  measurements  were  made 
from  November  of  1984  through  July  of  1987,  with  all  seasons  represented. 

9.2  Thermal  Environment  Criteria 

9-2.1  Temporal  Variaiions-  Satellite  systems  have  thermal  response  times  ranging  from  a 
few  minutes  to  hours.  To  enable  the  analyst  to  model  the  thermal  response  of  a  particular  system  as 
simply  as  possible,  a  running  mean  analysis  was  developed  for  the  albedo  and  OLR  variations  as  a 
vehicle  would  view  them  moving  along  its  orbit  track.  Data  are  presented  for  30*.  60*.  and  90*  orbit 
inclinations.  The  dependence  on  orbit  inclination  is  not  especially  strong,  so  the  data  presented  here 
are  expected  to  be  adequate  for  most  applications.  Running  averages  were  computed  for  averaging 
times  ranging  from  a  few  seconds  to  an  hour  and  a  half.  Percentile*  statistics  of  the  running  means 
describe  the  low  and  high  extremes,  1. 3. 5, 95, 97.  and  99  percentile,  and  the  expectation  value  (50 
percentile)  for  each  parameter  as  a  function  of  averaging  time  periods  along  an  orbit  track.  The  data 


*  A  peiwntile  value  is  the  probability  that  this  value  will  not  be  exceeded  (i.e.,  is  equal  to  or  less 
than  this  value).  ^ 


arc  adjusted  to  the  top  of  the  atmosphere  (30  km)  and,  for  albedo,  zero  solar  zenith  angle. t  (The 
diurnal  (solar  zenith  angle)  variations  of  OLR  were  not  significant  enough  to  warrant  similar 
treatment.)  These  statistics  are  presented  in  tables  9-2  through  9-4  and  are  illustrated  by  figures 
9-2  through  9-4. 

9.2.2  7i»nith  Angle  Correction  for  Albedo.  To  a  first  approximation,  albedo  may  be 
assumed  to  be  independent  of  the  solar  zenith  angle;  that  is,  the  scattering  is  “Lambertian”  or  equal 
in  all  directions.  This  approximation  has  been  assumed  for  most  engineering  applications  in  the  past. 
However,  the  data  quality  and  the  capability  of  current  engineering  analysis  methods  warrant  an 
improved  approach.  Albedo  actually  exhibits  an  appreciable  solar  zenith  angle  dependence.  Treat¬ 
ment  of  this  topic  in  the  scientific  literature  (ref.  30)  is  generally  “scene  specific,”  e.g.,  it  depends  on 
the  geographic  features  in  the  FOV,  data  which  are  not  readily  available  to  the  design  engineer,  and 
the  general  algorithms  tested  did  not  fully  remove  the  zenith  angle  dependence  from  this  data  set 
Therefore,  a  zenith  angle  correction  (fig.  9-1)  was  derived  2cifically  for  this  data  set  in  a  manner 
which  is  specifically  matched  to  the  analysis  tools  most  commonly  used  for  engineering  analysis,  the 
Thermal  Radiation  Analysis  System  (TRASYS)  and  Thermal  Synthesizer  System  (TSS). 

Most  satellite  orbits  encounter  all  solar  zenith  angles.  The  exceptions  are  certain  Sun-syn¬ 
chronous  orbits  which  never  reach  the  smaller  angles.  The  standard  TRASYS  and  TSS  codes  contain 
the  primary  solar  zenith  angle  dependence,  the  cosine  dependence  of  the  view  factor,  but  they  treat 
albedo  as  a  constant  Since  maximum  albedo  energy  occurs  at  small  zenith  angles  (near  but  not  at 
local  noon)  when  the  correction  factor  is  small  and  not  rapidly  changing,  it  will  often  be  satisfjwtory 
to  calculate  an  upper  limit  albedo  and  use  it  as  a  constant  in  the  TRASYS  or  TSS  code  for  estimating 
the  extreme  hot  case.  This  is  especially  true  for  subsystems  with  short  thermal  time  constants.  For 
analysis  of  the  temporal  variations  and  time  integrated  effects  of  the  thermal  environment,  the 
thermal  analysis  code  should  be  modified  to  account  for  the  albedo  variations.  Orbital  average  solar 
zenith  angles  for  the  daylight  portion  of  several  orbits  are  illustrated  in  figure  9-5.  For  any  single 
orbit* 


1^  ^  sza  ^  a/2  , 

that  is,  the  solar  zenith  angle  is  limited  by  the  beta  angle  at  local  noon  and  a/2  at  ground  sun¬ 
rise/sunset. 

The  correction  term  was  derived  from  4  mo  of  data  restricted  to  the  -30  to  +30  latitude  band. 
It  was  verified  by  testing  another  4  mo  of  data  and  testing  to  wider  latitude  bands.  It  removes  the 
solar  zenith  angle  dependence  to  within  ±0.04.  The  correction  is  expressed  as  an  additive  term,  thus: 

Albedo(SZA)  =  Albedo(S2:A  =  0)  +  Correction 

Correction  =  [C4(SZA)<  +C3(SZA)3  +  C2(SZA)2  +  Cl(SZA)] 

where  SZA  is  the  Solar  Zenith  Angle  in  degrees  and  the  albedo  is  expressed  as  a  fraction. 

C4  = +4.9115  E-9, 

C3  =  +6.0372  E-8, 

C2  = -2.1793  E-5, 

Cl  =  +1.3798  E-3. 


t  Solar  zenith  angle  is  the  angle  between  the  Earth  center-satellite  vector  and  the  Earth  center- Sun 
vector. 
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Hgure  9-1  illustrates  the  correction.  This  term  must  be  added  to  the  albedo  from  tables  9-2 
through  9-4  and  figures  9-2  through  9-4. 

9.2.3  Correlation  Analysis.  Figures  9-6  through  9-8  show  plots  of  the  running  mean  albedo 
against  the  running  mean  OLR  from  simultaneous  measurements.  The  two  parameters  are  not 
strongly  correlated  in  the  sense  that,  for  almost  any  value  of  one  parameter,  a  substantial  range  of 
values  is  encountered  for  the  other,  even  though  there  is  clearly  some  relationship  between  them. 
For  many  applications,  the  albedo  and  emitted  thermal  radiation  may  be  considered  as  independent 
parameters,  both  of  which  vary  across  a  range  of  values.  The  range  is  quite  large  for  systems  with 
short  thermal  time  constants,  e.g..  short  averaging  times,  but  as  longer  averaging  times  are  consid¬ 
ered,  the  range  narrows  toward  the  mean  for  the  region  covered  by  the  spacecraft 

In  selecting  benchmarks  for  design  purposes,  it  should  be  remembered  that  the  joint  probabil¬ 
ity  of  occurrence  for  independent  parameters  is  given  by  the  product  of  the  individual  probabilities. 
For  example.  OLR  values  in  excess  of  the  95th  percentile  value  may  be  expected  to  occur  in  com¬ 
bination  with  albedos  in  excess  of  the  mediati  0.025  of  the  time,  e.g.  (0.05)(0.50)  =  0.025.  Thus, 
assuming  these  are  independent  parameters  implies  97.5  percent  of  the  data  would  lie  outside  this 
region  in  this  example.  Eventually  joint  probability  distributions  for  these  parameters  will  be  devel¬ 
oped  for  applications  where  the  assumption  of  indepenoence  is  not  adequate,  but  this  analysis  is  not 
yet  available. 

Table  9-1.  Thermal  parameters  for  LEO. 


Solar  constant  at  1  AU 

(1  AU  =  mean  Earth-Sun  distance) 

1371±10W/m2 

Maximum  solar  energy  flux  (winter  solstice) 

(at  Perihelion,  +10  W/m^  included) 

1428  W/m2 

Minimum  solar  energy  flux  (summer  solstice) 

(at  Aphelion.  -10  W/rn^  included) 

1316  W/m2 

Characteristic  blackbody  temperature  of  solar  spectrum 
(Albedo  radiation  has  approximately  the  same 
^)ectrum.) 

5762  K 

Space  sink  temperature 

3  K  (absolute) 
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Table  9-2.  Running  mean  albedo  and  OLR  percentile  data  for  30*  inclination  orbits. 
Albedo;  (top  of  the  atmosphere.  30  km.  and  corrected  to  zero  solar  zenith  angle.*) _ 


*Angle  between  the  Earth  center-satellite  vector  and  the  Earth  center-Sun  vector. 


Outgoing  Longwave  Radiation;  (top  of  the  atmosphere,  30  km.  units  are  W/m^) 


*Set  at  0.04  percentile,  eliminating  a  very  few,  unrealistically  low  data  points. 


Note:  Percentile  is  the  probability  that  the  indicated  albedo  or  OLR  value  will  not  be  exceeded  (i.e., 
is  equal  to  o:  less  than  this  value). 
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Table  9-3.  Running  mean  albedo  and  OLR  percentile  data  for  60"  inclination  orbits. 


Albedo:  (top  of  the  aonosphere,  30  km.  Mid  corrected  to  zero  solar  zenith  angle.*) 


Ave. 
Time  (s) 

Lowest 

Observed 

1% 

3% 

5% 

50% 

95% 

97% 

99% 

Highest 

Observed 

16 

WMSESMM 

0.09 

0.10 

■301 

0.23 

0.39 

0.41 

0.44 

0.53 

64 

0.06 

0.09 

0.10 

0.11 

0.23 

0.39 

0.40 

0.43 

0.53 

128 

0.06 

0.09 

0.10 

0.11 

0.23 

038 

0.40 

0.43 

0.53 

256 

0.06 

0.09 

0.10 

0.11 

0.23 

037 

0.39 

0.42 

0J2 

512 

0.07 

0.10 

0.12 

0.13 

0.23 

0.36 

0.37 

0.40 

0.48 

896 

0.<^8 

0.12 

0.13 

0.14 

0.23 

0.34 

0.35 

0.37 

0.44 

1344 

0.10 

0.14 

■ilM 

WtBIM 

0.23 

0.33 

0.34 

0.35 

0.41 

1800 

0.11 

0.15 

0.16 

0.16 

0.23 

032 

0.33 

0.34 

0.40 

2688 

0.12 

0.15 

0.16 

0.17 

0.23 

0.32 

0.33 

0.34 

0.39 

3600 

0.13 

0.16 

0.17 

0.17 

0.23 

0.31 

0.32 

0.33 

0.38 

5400 

0.14 

0.16 

0.17 

0.18 

0.23 

0.31 

0.32 

0.33 

0.36 

♦Angle  between  the  Earth  center-satellite  vector  and  the  Earth  center-Sun  vector. 


Ave. 
Time  (s) 

Lowest  ♦ 
Observed 

1% 

3% 

5% 

50% 

95% 

97% 

99% 

Highest 

Observed 

16 

155 

178 

187 

192 

230 

279 

286 

300 

349 

64 

156 

179 

188 

193 

230 

279 

285 

300 

348 

128 

158 

181 

189 

193 

230 

278 

285 

299 

347 

256 

164 

185 

192 

196 

230 

276 

282 

296 

342 

512 

177 

191 

197 

200 

231 

270 

275 

287 

332 

896 

186 

197 

202 

205 

232 

264 

268 

278 

316 

1344 

191 

201 

206 

209 

232 

259 

263 

271 

303 

1800 

193 

203 

208 

211 

232 

255 

259 

267 

295 

2688 

196 

204 

209 

211 

232 

255 

258 

265 

284 

3600 

198 

205 

210 

213 

232 

253 

257 

263 

280 

5400 

200 

205 

211 

213 

233 

252 

255 

260 

278 

1 

« 


♦Set  at  0.04  percentile,  eliminating  a  very  few,  unrealistically  low  data  points. 


Note:  Percentile  is  the  probability  that  the  indicated  albedo  or  OLR  value  will  not  be  exceeded  (i.e.. 
is  equal  to  or  less  than  this  value). 


Table  9-4.  Running  mean  albedo  and  OLR  percentile  data  for  90*  inclination  orbits. 


A|b(^o;  (top  of  the  atmosphere.  30  km.  and  corrected  to  zero  solar  zenith  angle  *) 

Ave.  j  Lowest  I  T - 1 - ] - f- ' 

I  Time  (s)  |  Observed  I  1%  3%  I  5%  I  I  I  oncL.  I  rvw 


Highest 


♦Angle  between  the  Earth  center-satellite  vector  and  the  Earth  center-Sun  vector. 


Outgoing  Longwave  R«Uation;  (top  of  the  atmosphere.  30  km.  units  are  W/m2t 


Ave.  Lowest  * 
Time  (s)  Observed 


♦Set  at  0.04  percenUle,  eliminating  a  very  few,  unrealistically  low  data  points. 


Note:  Percentile  is  the  probability 
is  equal  to  or  less  than  this 


Highest 

Observed 


349 


348 


347 


342 


332 


315 


304 


that  the  indicated  albedo  or  OLR  value  will  not  be  exceeded  (Le., 

V2UU6}. 


ALBE 


LONGWAVE  RADIATION  (W/m2)  ALBEDO 
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AVERAGING  TIME  (min) 


AVERAGING  TIME  (min) 


Rgure  9-3.  Statistical  summary  of  albedo  (top)  and  OLR  (bottom)  observations  for  a  60*  inclination 
orbit  Solid  lines-limiting  values;  dashed  lines-1  and  99  percent;  dotted  lines-3  and  97  percent; 
dot  dasli-5  and  95  percent;  and  ‘V’-median, 
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M 


Orbit  Number 


ALTaSOOKM:  INC=>60» 


IK 
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Orbit  Number 


figure  9-5.  Orbit-average  solar  zenith  angles  for  30*.  60*.  and  90*  inclination  orbits.  dayUght  side  only 
Annual  average  is  indicated  by  soUd  horizontal  Une.  Phase  of  the  peaks  is  a  function  of  fi  angle. 
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Figure  9*6.  Albedo-OLR  correlated  pairs  for  30*  inclination  orbits  (page  1  of  2). 
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Figure  9-7.  Albedo-OLR  correlated  pairs  for  60*  inclination  orbits  (page  1  of  2). 
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Rgure  9-7.  Albedo-OLR  correlated  pairs  for  60*  inclination  orbits  (page  2  of  2). 
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X.  GRAVITATIONAL  FIELD 


With  the  advent  of  Earth  satellites,  there  has  been  a  considerable  advance  in  the  accurate  deter¬ 
mination  of  the  Earth’s  gravitational  field.  The  current  knowledge  regarding  the  Earth’s  gravitational 
field  has  advanced  far  beyond  the  normal  operational  requirement  Adequate  accuracy  for 
most  space  vehicle  design  values  of  gravitational  inieracUons  is  obtained  with  the  central  inverse  sqiuue 


where: 


(10-1) 


? = vector  force  acting  on  a  particle  in  Newtons, 
m  =  mass  of  particle  in  kg, 

= gravitational  constant  for  the  Earth, 
r  =  distance  from  the  particle  to  the  center  of  the  Earth  in  meters, 
f  =  unit  vector  from  the  center  of  the  Earth  to  the  particle,  dimensionless. 


The  a^ve  central  force  model  accurately  represents  the  gravitational  field  to  approximaielv 
0.1  percent  If  thrs  accuracy  is  insufficient,  a  more  detaUed  model  of  the  gravitational  field  can  be  ^d 
that  amounts  for  the  non-uniform  mass  distribution  within  Earth.  This  model  gives  gravitational 
potential,  V,  to  an  accuracy  of  approximately  a  few  parts  in  a  miUion.  The  gravitational  acceleration  is 
the  negauve  gradient  (vector  derivative)  of  the  potential  and  introduces  nonradial  gravitational  forces. 


?=#«2=in(-^V)  . 


The  formula  for  V  is  shown  below; 


(10-2) 


V(r,^,A)-4^(l+JjJ^(i^)  i»^sin^)[C^cos(mA)-hS^sin(mA)j) ,  (10-3) 

where  ■' 

^ = geocentric  decliruition 

A  =  east  longitude 

He  -  gravitatimial  constant  for  the  Earth 
for  the  Marsh  et  al.^^  model  discussed  below 

/t£=3.9860064EM(/Vra2/kg)  ,* 

R£= 6378.137  km.* 

♦Use  these  values  only  in  this  model.  For  other  ^plications,  use  the  values  found  in  table  2-1. 
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and  5,^  are  the  hannonic  coefficients  of  the  potential  fiinctioa  and  P 
repi^^the  a«^ated  Legendre  functions  of  the  flrst  kind  of  degree  n  and  order  hl  BecaiM^  "" 
sm  ♦= Z» = «. «  Mmply  a  coai«  and  an«da«d  as: 


/»-  = 


2"n! 


d/u* 


(10-4) 


model  widi  the  number  of  coefficients  tailored  to  the  application,  rather  than  truncate  a  larm-  «»t  a 
lecommCTded  high  accuracy  set  of  coefficients  {36x^  is  prei^ted  rSS^l jSS.5 
mformauon  on  normalization  can  be  found  in  Anon.  (1975L  3*-  Addmonal 

Table  10-1.  Gravitational  coefficients  (4x4). 


Zonal  Harmonic  Coefficients  (Multiply  by  10-4)  j 

CIO 

0.0 

C20 

-1082.6258313 

C30 

2J326124 

C40 

1.6161966 

Tesseral  Hannonic  Coefficients  (Multiply  by  10*4) 

n 

m 

Cim 

^nm 

1 

1 

0.0 

0.0 

2 

1 

0.0 

0.0 

2 

2 

+1.5743213 

-0.90359264 

3 

1 

+2.1924062 

+0.2695930 

3 

2 

+0.3086208 

-0.2119137 

3 

3 

+0.1005368 

+0.1970571 

4 

1 

-0.5060535 

-0.4507374 

4 

2 

+0.0775920 

+0.1484817 

4 

3 

+0.0592223 

-0.0119894 
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ABBREVIATIONS  AND  ACRONYl^ 

ac  alternating  current 

AL  anomalously  large 

AO  atomic  oxygen 

\  geomagnetic  activity  index  (daily) 

geomagnetic  activity  index  (3-hourIy) 

AU  astronomical  unit 

C  Celsius 

CCD  charge-coupled  device 

CH  coronal  hole 

cm  centimeter 

CME  coronal  mass  ejection 

CMG  control  moment  gyro 

dc  direct  current 

DMSP  Defense  Meteorological  Satellite  Program 

EMI  electromagnetic  interference 

EMR  electromagnetic  radiation 

EUV  extreme  ultraviolet 

eV  electron  volt 

EVA  extravehicular  activity 

Eio.7  1 0.7  centimeter  solar  radio  noise  flux 

FORTRAN  formula  translation 

G  gauss 

GCR  galactic  cosmic  ray 
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GEO 

GeV 

GHz 

GN&C 

GRAM 

GSFC 

GV 

Hg 

Hz 

IGRF 

IMF 

IOC 

IR 

IRI 

JAS 

JGR 

JSC 

K 

keV 

H 

kHz 

km 

LDEF 


geosynchronous  Earth  orbit 
gigaelectron  volt 
gigahertz 

gravitational  constant  for  the  Earth 

guidance,  navigation,  and  control 

global  reference  atmosphere  model 

Goddard  Space  Flight  Center 

gigavolt 

mercury 

Hertz 

international  geomagnetic  reference  field 
interplanetary  magnetic  field 
initial  operational  capability 
infrared 

international  reference  ionosphere 
Journal  of  Atmospheric  Sciences 
Journal  of  Geophy»cal  Research 
Johnson  Space  Center 
Kelvin 

kiloelectron  volt 
kUogram 
kiloHertz 
kUometer 

logarithm  of  the  geomagnetic  activity  index 
Long  Duration  Exposure  Facility 


* 

i 


.EO 

^RC 

-y-a 

0 

/lET 

4eV 

4Hz 

4LT 

4SFC 

iV 

fW 

lASA 

m 

IRL 

rSSDC 

m 

)R 

CA 

£ 

F 

i\M 

EU 

ZA 


low-Earth  orbit 
Lewis  Research  Center 
large-scale  integration 
Lyman-alpha 
meter 

Marshall  engineering  thermosphere 

megaelectron  volt 

megahertz 

magnetic  local  time 

Marshall  Space  Flight  Center 

millivolt 

megawatt 

National  Aeronautics  and  Space  Administration 
nanometer 

Naval  Research  Laboratory 
National  Space  Science  Data  Center 
New  Technology,  Inc. 
ordinary 

polar  cap  absorption 
radius  of  the  Earth 
radio  frequency 
second 

statistical  analysis  mode 
single  event  upset 
s<dar  zenith  angle 


TH) 

UAH 

UT 

u 

UV 

VLF 

VLSID 

W 

wrt 


traveling  ionospheric  disturbance 

University  of  Alabama,  Huntsville 

universal  time 

atomic  mass  unit 

ultraviolet 

very  low  frequency 

very  large-scale  integrated  device 

Watt 

with  re^ct  to 
westward  traveling  surges 


WTS 
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APPENDIX  B 


1.0  GALACTIC  COSMIC  RAYS  (OCR’s)* 

The  following  paragraphs  and  tables  present  an  algodthm  for  com(^g  the  differential  energy 
spectra  of  the  most  important  charged-particle  populations  in  the  Earth’s  vicinity-  These  equations  were 
de^dsed  to  fit  the  data  and  are  intended  to  have  no  physical  interpretation.  This  analytic  formula  may  be 
easUy  programmed  for  a  digital  computer  of  almost  any  size  and  is  intended  to  become  a  subroutine  in  a 
program  which  will  be  used  to  estimate  the  soft  error  rales  in  satellite-brvne  electronics. 

OCR’s  consist  of  electrons  and  the  nuclei  of  all  die  dements  in  the  periodic  table;  the  first  28 
elements  are  the  most  important  for  cosmic  ray  effects  on  microelectronics.  These  particles  are  from 
outade  the  solar  system,  and  their  flux  at  low  energies  is  anticonelated  with  solar  activity  (i.e.,  more 
cosmic  rays  at  solar  minimum).  The  differential  energy  spectra  in  particles  per  square  meter-steradian- 
second-Megaelectron  volts  per  atomic  mass  unit  (i.e.,  particlcs/(m^'sr-s*MeV/u))  are  given  in  the 
following  paragraphs. 

The  spectra  for  protons  (hydrogen  nuclei),  a-particles  (helium  nuclei),  and  iron  nuclei  are  given 
below  for  energies  above  10  MeV/u: 

F(£,r) = A(£)  sinnir(r-q))l+B(£) 

where 

W  =  0.576  radian/year. 

tQ  =  1950.6  A.D.  date. 

t  =  current  date  in  years, 

E  =  particle  energy  in  MeV/nucleon, 

B(E)  =  0.5 

A(E)  =  03\f^{E^^m. 
and/^  differ  only  by  the  choice  of  constants  in  tbeequatioo. 

/(£)=10«(fi«b)*  . 

where 

fl  =  ao(  l-exp[-Xiaogio  £)*] ) . 
and 

#»  =  CjCxp[-X2flogjo£)^-C2-  (Bl-6) 


(Bl-1) 

(Bl-2) 

(Bl-3) 

(Bl-4) 


•Part  1,  appendix  B  is  from  Adams  et  al.^^  ^ 


The  values  of  the  constants,  oq,  Eq^  b,  X\,  X2,  Ci,  and  C2,  are  given  in  table  Bl-1  for  each  of  the 

elements  hydrogen  (H),  helium  (He),  and  iron  (Fe)  for  the  conditions  of  solar  maximum  and  solar 
minimum. 


Table  Bl-1.  Constants  used  in  equations  (Bl-4)  through  (Bl-6)  to  compute  the  differential  energy 
spectra  of  H,  He,  and  Fe  at  solar  maximum  and  solar  minimum. 


Element 

<*0 

b 

^1 

— 

mam 

C2 

H-min 

r^2.20 

1.1775E+5 

1685 

0.117 

0.80 

6.52 

4.00 

K-max 

-2.20 

1.1775E+5 

2.685 

0.079 

0.80 

6.52 

4.00 

He-min 

-2.35 

2.070 

0.241 

0.83 

4.75 

5.10 

He-max 

-2.35 

8.2700E+4 

2.070 

0.180 

0.83 

4.75 

5.10 

Fe-min 

-2.14 

1.1750E+5 

1640 

0.140 

0.65 

6.63 

7.69 

Fe-max 

-2.14 

1.1750E+5 

2.640 

0.102 

0.65 

6.63 

7.69 

The  differenual  energy  spectra  for  carbon  (C),  oxygen  (O),  fluorine  (F),  neon  (Ne),  sodium  (Na), 
aluminum  (Al),  and  phosphorus  (P)  are  obtained  by  multiplying  the  helium  spectrum  (obtained  from 
equation  (Bl-1))  by  the  appropriate  scaling  factor  in  table  Bl-2. 

Table  Bl-2.  The  ratio  of  the  abundance  of  various  nuclei  to  helium. 


Element 

Ratio 

Element 

Ratio 

C 

3.04E-2 

0 

2.84E-2 

Al 

1.07E-3 

F 

6.06E-4 

Ne 

4.63E-3 

P 

2.34E-4 

Na 

HEESSI 

The  differential  energy  spectra  for  calcium  (Ca),  cobalt  (Co),  and  nickel  (Ni)  are  obtained  by 
multiplymg  the  iron  spectrum  (obtained  from  equation  (Bl- 1))  by  the  scaling  factors  listed  in  table  Bl-3. 

Table  Bl-3.  The  ratios  of  the  abundance  of  various  elements  to  iron. 
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The  spectra  of  the  elements  lithium  (Li),  beryllium  (Be),  and  boron  (B)  are  obtained  from  the 
helium  spectrum,  modified  by  the  equation: 


'  0.021FHe 

pni9E-^**^F^ 


£<3000MeV/u 
£^3000MeV/u  ’ 


^1-7) 


to  obtain  the  combined  spectrum  of  (U+Be+B).  Equation  (Bl-7)  is  then  multipUed  by  the  ratio  in  table 
Bl-4  to  obtain  the  individual  elemental  spectra. 


Table  Bl-4.  Relative  fractions  of  Li,  Be,  and  B  in  the  combined  total  abundance  Li-fBe+B. 


Element 

Ratio 

Li 

0.330 

Be 

0.176 

B 

0.480 

The  spectrum  of  the  element  nitrogen  (N>  is  obtained  by  modifying  the  helium  spectrum,  Fhc,  as 
shown  below: 


Fyv=  {8-7£-3  expI-0.4  (log,o£-3.15)2]+7.6£-3exp[-0.9  (logio£-0.8)2]}FHe  (Bl-8) 

where  £  is  in  MeV/u. 

! 

The  Spectra  of  the  elements  magnesium  (Mg),  silicon  (Si),  and  sulfur  (S)  are  Obtained  by 

modifying  the  helium  spectrum,  as  shown: 


F*  = 


^He 

,(l+1.56£-^£-2200)F„e 


£<2200 
£^2200  ’ 


(Bl-9) 


The  individual  spectra  for  these  elements  are  obtained  by  multiplying  F*  by  the  ratio  of  table 


Table  Bl-5.  Ratios  of  Mg,  Si,  and  S  to  an  adjusted  helium  spectrum. 


Element 

Ratio 

Mg 

6.02E*3 

Si 

4.63E-3 

S 

9.30E-^ 
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The  spectra  for  the  elements  chlorine  (Cl),  argon  (Ar),  potassium  (K),  scandium  (Sc),  titanium 
(TO,  vanadium  (V),  chromium  (Cr),  and  manganese  (Mn)  are  all  obtained  by  modifying  the  iron 
spectrum,  as  shown  below: 

F*~Q(E)Ff^.  (Bl-10) 

where 

Q(E)  =  16  [l-exp(-0.075£0-»)l£^.33 ,  (Bl-ll) 

where  £  is  in  MeV/u. 

The  £*,  from  the  subiron  spectrum  (equation  (Bl>10))  is  multiplied  by  the  appropriate  ratio  in 
table  Bl-6  to  obtain  the  individual  elemental  spectra. 

Table  Bl-6.  Fractional  abundance  of  each  element  in  the  subiron  group. 


Element 

Ratio 

Element 

Ratio 

Cl 

0.070 

Ti 

0.1-:? 

At 

0.130 

V 

0.070 

K 

0.090 

Cr 

0.140 

Sc 

0.042 

Mn 

0.100 

The  differential  energy  spectra  for  elements  from  copper  to  uranium  are  obtained  by  multiplying 

(he  iron  spectrum  (from  equation  (Bl-D)  by  the  scaling  factors  listed  in  table  Bl-7. 

Table  Bl-7.  Ratio  of  the  abundances  of  various  nuclei  to  iron  (page  1  of  2). 


Element 

Ratio 

Element 

Ratio 

Cu 

6.8E^ 

Pm 

1.9E*7 

Zn 

WMMM 

Sm 

8.7E-7 

Ga 

6.5E-5 

Eu 

1.5E-7 

Ge 

1.4E^ 

Gd 

7.0E-7 

As 

9.9E-6 

Tb 

1.7E-7 

Se 

5.8E-5 

Dy 

Br 

8.3E-6 

Ho 

2.6E*7 

Kr 

Er 

4.3E-7 

Rb 

l.lE-5 

Tm 

■IUBI 

Sr 

3.6E-5 

Yb 

4.4E*7 

Y 

6.8E-« 

^  Lu 

^  6.4E-* 

B-S 


Table  B 1-7.  Ratio  of  the  abundances  of  various  nuclei  to  iron  (page  2  of  2). 


Element 

Ratio 

Element 

Ratio 

Zr 

1.7E*5 

Hf 

4.0E-7 

Nb 

2.6E-<5 

Ta 

3.6E-8 

Mo 

7.1E-6 

W 

3.8E-7 

Tc 

1.6E-6 

Re 

1.3E-7 

Ru 

5.3E'6 

Os 

5.6E-7 

Rh 

1.5E-6 

Ir 

3.7E-7 

Pd 

4.5E-6 

Pt 

7.2E-7 

Ag 

1.3E-<i 

Au 

1.3E-7 

Cd 

3.6E-« 

Hg 

2.3E-7 

In 

1.4E-6 

T1 

1.8E-7 

Sn 

7JE-6 

Pb 

1.7E-« 

Sb 

9.9E-7 

Bi 

9.0E-* 

«  w 

5.7E-6 

Po 

0 

I 

1.5E-6 

At 

0 

Xe 

3.5E-6 

Rn 

0 

Cs 

5.8E-7 

Fr 

0 

Ba 

6.0E-6 

Ra 

0 

La 

5.3E-7 

Ac 

0 

ce 

1.6E-6 

111 

9.0E-* 

Pr 

3.0E-7 

Pa 

0 

Nd 

l.lE-6 

u 

5.4E-8 

The  formula  given  above  is  correct  for  quiet  periods  in  the  interplanetary  medium  when  only  the 
GCRs  are  present  These  conditions  arc  often  disturbed,  especially  at  low  energies,  by  small  solar  flares, 
co-rotating  events,  etc.  To  allow  for  typical  disturbed  conditions,  a  worst-case  spectrum  should  be 
employed.  With  90-percent  confidence,  the  instantaneous  particle  flux  should  never  be  more  intense 
than  described  by  this  case  at  any  energy. 

To  construct  the  worst-case  spectrum  for  protons,  compute  the  “H-min”  spectru'.n  (using 
equation  (Bl-4))  and  then  compute  as  shown  below: 

.  (Bl-12) 

This  applies  for  £  <  100  MeV/u.  For  higher  energies,  use  the  GCR  spectrum  for  the  appropriate  mission 
time,  r,  in  equation  (Bl-1). 
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U  S  "pSyf  «*“!»«<  fto™  eq-Bon 


28.4e-»l3.84+i.64  ^ 


for££lOOMeV/u. 


(Bl-13) 


a  muldpfc  '''  ^  "«  b, 

^woist  =  1-64  Fg^^ .  (Bl-14) 

»?««.  i^..  i. 

In  addition  to  galactic  cosaiic'rays,  some  particles  ate  believed  to  be  accelerated  in  rhp  intor 

STed^  ““"O'  “*  "y  '«B™>  fecmmXM'be 

'■  S'Si^td”iSrcSL‘S«  flTte  ® 

2.  equation  (BM)  so  that  these  maximuin  values  apply  for  aU  energies  below 

the  energy  at  which  the  maxima  occurs.  i.e..  for  solar  minimum:  ^ 


^  0.33 

l.^ie-iiiu  (e^.  (Bl-4) 


£<200MeV/u 

£^200MeV/u 


(BUIS) 


i.  Make  the  same  type  of  modification,  fw  solar  maximum: 


f*  _/  0.33 


£<300MeV/u 

£S300MeV/u 


(Bl-16) 


4.  Combine  the  resulting  spectra  as  before  using  equations  (Bl-1)  through  (Bl-3). 


f(E)  =  6£-2  cxp[-(ln  (£>-1.79)2/0.70]  (particles/(m2-sr*s  MeV/u)  . 


(Bl-17) 
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This  spectrum  crosses  over  the  galactic  spectrum  at  approximately  30  MeV/u.  The  two  ^Kctra  should  be 
matched  at  that  point  with  equation  (Bl-17)  replacing  the  galactic  spectrum  at  lower  wietpes. 

Similarly  for  nitrogen: 

/(£)  =  1.54E-2  expKln  (£)-1.79)2/0.70J  (particles/(m2-srs-MeV/u) .  (Bl-18) 

Again,  this  crosses  with  the  GCR  spectrum  at  approximately  30  MeVAi  and  should  replace  it  below  this 
energy.  The  spectra  of  the  remaining  elements  are  unaffected  or  affected  at  too  low  an  energy  to  matter. 

There  is  a  possibility  that  the  anomalous  component  is  singly  ioniaed.  If  so,  it  will  have  an 
extraordinary  ability  to  penetrate  the  Earth's  magnetosphere.  In  this  case,  the  differential  energy  ^)ectra 
shown  below  are  assumed  for  the  helium,  carbon,  nitrogen,  oxygen,  neon,  magnesiom,  silicon,  argon, 
and  iron  spectra  of  the  anomalous  component  There  probably  are  anomalous  components  in  the  spectra 
of  the  nuclei  heavier  than  iron,  but  there  are  no  data  on  them  at  this  time. 

For  singly  ionized  helium. 


(Bl-19) 


(Bl-20) 


(BI-21) 


(Bl-22) 


(Bl-23) 


(Bl-24) 
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For  singly  ionized  silicon. 


^  _  flOOxW^  expt-ln(£>-2.20)V0.4] 

£<10 

1  0.10£-2 

£S10  ‘ 

(Bl-25) 

For  singly  ionized  argon. 

^  _  15.40x10“*  exp(-ln(£)-1.79)  VO.7] 

£<20 

\  0.28£-2 

£^20  * 

(Bl-26) 

For  singly  ionized  iron. 

^  _  /6.00xl0“*  exp[-4n(£)-2.48)^/2.0] 

£<30 

~\  035£-^ 

£^30  ’ 

(Bl-27) 

Since  these  Moraalous  component  particles  are  assumed  to  be  singly  ionized,  they  will  have  a  hieher 
magnettc  ngidity  than  galactic  cosmic  rays  of  the  same  energy.  The  magnetic  rigidity  of  galactic^mic 

rflys  IS 


R  =  (A/Z)(£2+1862.324£)i«/iooo  , 


(Bl-28) 


in  GeV/ec.  The  rigidity  of  singly  ionized  nuclei  is 


^  =  (A(£2+1862324£)i'2)/1000  . 


(Bl-29) 


To  add  the  smgly  lom^  anomalous  component,  it  is  necessary  to  modulate  both  the  galactic  cosmic  rav 
spectra  and  the  anomalous  component  spectra  given  above  using  the  geomagnetic  cutoff  transmission 
areS *ad??t^^^Ui**™^  (discussed  in  Adams  et  ai.3^  35),  the  lesuldn^  modulated  spectra 
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2.0  SOLAR  n.ARE  PARTICLE  EVENTS* 

Solar  flare  particle  events  are  sporadic  occurrences  lasting  1  to  5  days.  When  these  events  occur, 
they  can  be  the  dominant  cause  of  soft  errors.  For  statistical  treatment,  they  are  broken  into  two 
ordinary  (OR)  and  anomalously  large  (AL).  The  probability  of  having  mote  than  a  number  of  events,  n. 
in  a  time,  r.  is  given  by: 


P{n^^,T)  =1-2  (i+AOKr/Di/IiWld+r/T)'  .  (B2-i) 

where  T and  r  are  in  years,  and  N  is  the  number  of  flares  that  have  occurred  in  Tyears. 

For  ordinary  events,  equation  (B2- 1 )  becomes: 

^OR  =  P(n,r,24,7)  for  the  decreasing  portion  of  the  solar  cycle 
and  (B2-2) 

^OR=  P(n,r,6,8)  for  the  increasing  portion  of  the  solar  cycle, 

where  there  is  a  probability,  Por>  of  having  more  than  n  ordinary  events  in  t  years.  Similarly  for 
anomalously  large  events: 

P^  =  P{n.tXl)  .  (B2-3) 

The  peak  proton  flux  differential  energy  spectrum  for  ordinary  events  is,  typically: 

/oR  =  2.45F*(e-^^-5+173e^)protons/(m2  sr  s  MeV)  ,  (B2-4) 

where  £  is  in  MeV,  and  no  worse  than: 


/woR  =  2.06£5(e-^^-5+63.6  tr^^)  proions/Cm^-sr  s-MeV)  ,  (B2-5) 

with  a  confidence  of  approximately  90  percent 

Using  the  August  1972  flare  as  a  model  AL  event,  the  peak  proton  flux  differential  energy 
spectrum  is: 


fKL- 


9.3£V£/</£)exp(-P/0.10) 

\.16E\dPidE)p-^ 


£<150  MeV 
£^150  MeV  ’ 


(B2-«) 


in  protons/fm^sr-s-MeV),  where 


P=[£Vl863.2£]‘'^/1000  , 


(B2-7) 


and  £  is  in  MeV. 

♦Part  2,  appendix  B  is  from  Adams  et  al.^^  34 


\ 
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To  model  the  worst  flare  that  is  ever  likely  to  occur,  use  the  composite  of  the  August  1972  flare 
and  the  February  1956  flare.  The  composite  worst-case  flare  proton  spectrum  is  taken  to  be  the  peak  of 
the  1972  spectrum,  as  given  by  equation  (B2-6),  and  the  1956  peak  spectrum,  given  below. 

/i956  =  1.116£+«(£-»-248)(0.248+2.5£*5)((i.7)(epowKEXPOW)) 

+  4.7£+*9(£-5^)(4,3(1-EXPOW>-^.32£-i5((4){EPOW)(EXPOW)))  ,  (B2-8) 


where 

and 


EP0W  =  £‘7 

EXPOW = exp(-2.5£-5  EPOW) . 


The  composition  of  flare  particles  is  also  highly  variable  from  flare  to  flare.  Table  B2-1  gives  the 
composition  relative  to  hydrogen  for  the  elements  through  nickel.  Both  mean  and  90-percent  confldence 
level  worst  cases  are  given.  To  obtain  the  spectrum  of  any  element  in  a  flare,  just  multiply  the  abundance 
ratio  from  the  Table  B2-1  by  the  appropriate  flare  proton  spectrum. 


The  worst-case  compositions  of  the  elements  fiom  copper  to  uranium  are  obtained  by 
multiplying  the  abundance  ratios  of  table  B2-2  by: 

(C«<O)/C,(O))0.48  exp(Z0-78/6.89) .  (B2-9) 


where  CjvfO)  and  Q(0)  are  the  worst-case  and  mean  abundance  coefficients  for  oxygen  in  Table  B2-1. 


Table  B2-1.  Mean  and  worst-case  flare  compositions. 


Element 

Mean  Case 

Worst  Case 

H 

1 

1 

He 

I.OE-2 

3.3E-2 

Li 

0 

0 

Be 

0 

0 

B 

0 

0 

C 

1.6E-4 

4.0E-4 

N 

3.8E-5 

l.lE-4 

0 

3.2E^ 

1.0E*3 

F 

0 

0 

Ne 

5.1E-5 

1.9E-4 

Na 

3.2E-6 

1.3E*5 

Mg 

6.4E-5 

2.5E-4 

A1 

3.5E-6 

1.4E*5 

Si 

5.8E-5 

1.9E-4 

P 

2.3E-7 

l.lE-6 

S 

8.0E-6 

5.0E-5 

Cl 

1.7E-7 

8.0E-7 

Ar 

3.3E-6 

1.8E-5 

K 

1.3E-7 

6.0E-7 

Ca 

3.2E-6 

2.0E-5 

Sc 

0 

0 

Ti 

l.OE-7 

5.0E-7 

V 

0 

0 

Cr 

5.7E-7 

4.0E-<5 

Mn 

4.2E-7 

2.3E-6 

Fe 

4.IE-5 

4.0E'4 

Co 

l.OE-7 

5.5E-7 

Ni 

2.2E-« 

HE^EH 

The  mean  case  compositions  for  the  elements  from  copper  to  uranium  arc  taken  from  reference 
36.  The  ratios  of  these  abundances  to  hydrogen  are  given  in  Table  B2-2. 


0 
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3.0  GEOMAGNETIC  CUTOFFS* 

The  modulation  of  cosmic  spectra  by  the  Earth’s  nu^netic  field  requires  a  mwe  thorough 

treatment  than  can  be  offered  here,  but  some  guidance  will  be  provided.  The  geomagnetic  cutoff  is  a 
value  of  magnetic  rigidity  below  which  cosmic  rays  will  not  reach  a  specified  point  in  the  magneto¬ 
sphere  from  a  specified  direction.  The  magnetic  ri^ty,  P,  in  GeV/ec,  may  be  computed  fiom  the 
particles*  energy  using: 


^  =  .  (B3-1) 

where  A  and  Z  are  the  atomic  ma»  and  charge  of  the  nucleus  in  question.  The  cutoff  at  any  point  for 
particles  arriving  from  the  zenith  is  most  simply  computed  with: 

Pc=15.96/f.2-«»C.eV/ec  .  (B3-2) 

where  L  is  McDwain’s  L  parameter  fte..  the  radial  distance,  in  Earth  radii,  from  the  center  of  the  Earth  to 
the  point  in  the  geometric  equatorial  plane  where  it  is  crossed  by  the  magnetic  field  line  that  also  passes 
through  the  point  of  observation). 

Design  consideration  must  be  given  to  solar  flare  spectra,  because  (1)  the  flare  particle  intensiw 
changes  on  a  time  scale  comparable  to  or  shorter  than  an  orbital  period,  (2)  there  is  no  certain  proof  that 
solar  flare  particles  are  fully  ionized,  and  (3)  the  geomagnetic  cutoff  is  suppressed  to  some  extent  during 
a  flare.  The  geomagnetic  cutoff  during  a  flare,  Pf,  should  be  computed  from  the  “quiet  time”  cutoff,  P©, 
using: 


5P/Po  =  0.54exp(-Po/2.9)  ,  (B3-3) 

and 

Pp  =  Po-SP,  (B3-4) 

where  Pp,  Pq,  and  ^  are  in  GeV/ec. 


•Part  3,  Appendix  B  is  from  Adams.3< 
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APPENDIX  C 

Table  C-1.  Median  value  of  global  maximum  densities  for  altitude  and  fio.?  bin  with 
differences  between  global  maximum  density  and  median  value 
for  several  percentile  ranges.  Densities  in  k^m^  (page  1  of  6). 


F|07b  Bin  1  Bin  2  Bin  3  Bin  4  Bin  5  All 

Range  66-102  102-138  138-174  174-210  210-246  66-246 


Altitude  250  km 


Median 

5.941E-11 

8.41  IE- 11 

1.050E-10 

1.239E-10 

1.363E-10 

8.629E-11 

1 

-1.554E-11 

-2.191E-11 

-2.827E-11 

-3.497E-11 

-2.560E-11 

-4.066E-11 

5 

-1.209E-11 

-1.707E-11 

-2.063E-11 

-2.48  IE- 11 

-2.000E-11 

-3.559E-11 

33 

-3.760E-12 

-4.970E-12 

-5.570E-12 

-6.400E-12 

-5.100E-12 

-1.718E-H 

67 

4.460E-12 

5.750E-12 

6.000E-12 

6.100E-12 

6.100E-12 

1.921E-11 

95 

1.816E-11 

2.359E-11 

2.140E-11 

2.090E-11 

2.130E-11 

5.291E-11 

99 

2.724E-1 1 

3.529E-11 

3.060E-11 

2.940E-11 

2.680E-11 

6.621E-11 

100 

6.099E-11 

6.859E-11 

4.670E-11 

4.590E-11 

3.740E-11 

8.741E-11 

Altitude  275  km 


Median 

3.123E-11 

4.72  IE- 11 

6.184E-11 

7.595E-11 

8.574E-1 1 

4.869E-11 

1 

-9.330E-12 

-1.425E-11 

-1.956E-11 

-2.546E-11 

-1.975E-11 

-2.576E-11 

5 

-7.310E-12 

-1.119E-11 

-1.446E-11 

-1.837E-11 

-1.554E-11 

-2.277E-11 

33 

-2.320E-12 

-3.330E-12 

-4.020E-12 

-4.880E-12 

-4.040E-12 

-1.135E-11 

67 

2.780E-12 

3.920E-12 

4.400E-12 

4.780E-12 

5.030E-12 

1.350E-11 

95 

1.161E-11 

1.660E-11 

1.607E-11 

1.683E-11 

1.796E-11 

"3.946E-11 

99 

1.770E-11 

2.534E-11 

2.333E-11 

2.405E-11 

2.286E-11 

5.060E-11 

100 

4.209E-11 

5.223E-11 

3.680E-11 

3.875E-11 

3.266E-11 

6.97  IE- 11 

Altitude  300  km 


Median 

1.721E-11 

2.767E-11 

3.793E-11 

4.843E-11 

5.609E-11 

2.868E-11 

1 

-5.720E-12 

-9.360E-12 

-1.356E-11 

-1.851E-11 

-1.513E-11 

-1.658E-11 

5 

-4.510E-12 

-7.410E-12 

-1.014E-11 

-1.355E-11 

-1.201E-11 

-1.476E-11 

33 

-1.450E-12 

-2.240E-12 

-2.880E-12 

-3.700E-12 

-3.200E-12 

-7.570E-12 

67 

1.760E-12 

2.690E-12 

3.210E-12 

3.700E-12 

4.060E-12 

9.510E-12 

95 

7.530E-12 

1.169E-11 

1.201E-11 

1.336E-11 

1.491E-11 

2.934E-11 

99 

1.164E-11 

1.818E-11 

1.771E-11 

1.940E-11 

1.915E-11 

3.854E-11 

100 

2.922E-1 1 

3.968E-11 

2.874E-11 

3.232E-11 

2.801E-11 

5.542E-11 

Altitude  325  km 


Median 

9.838E-12 

1.679E-11 

2.402E-11 

3.183E-11 

:.777E-11 

1.748E-11 

1 

-3.582E-12 

-6.250E-12 

-9.480E-12 

-1.349E-11 

-1.154E-11 

-1.085E-11 

5 

-2.839E-12 

-4.970E-12 

-7.150E-12 

-l.OOOE-ll 

-9.230E-12 

-9.725E-12 

33 

-9.230E-13 

-1.530E-12 

-2.070E-12 

-2.800E-12 

-2.510E-12 

-5.100E-12 

67 

1.142E-12 

1.850E-12 

2.340E-12 

2.840E-12 

3.240E-12 

6.730E-12 

95 

4.952E-12 

8.270E-12 

8.%0E-12 

1.051E-11 

1.221E-11 

2.183E-11 

99 

7.762E-12 

1.308E-11 

1.339E-11 

1.549E-11 

1.585E-n 

2.932E-11 

100 

2.046E-11 

3.012E-11 

2.232E-11 

2.662E-11 

2.368E-11 

4.397E-11 

value  of  global  maximum  densities  for  altitude  and  F,o  7  bin  with 
differences  between  global  maximum  densityand  median  value  for  «wial 
percentile  ranges.  Densities  in  kghn^  (page  2  of  6). 


Bin  4 


-2.283E-12 


-5.980E-13 


Median 


1 


5 


33 


67 


95 


99 


100 


Altitmte  375  km 


MedianI  3.479E-12 


3.307E-12 


-4.223E.12 


-3.382E.12 


-1.057E-12 


1.290E-12 


5.900E-12 


9.470E-12 


2.291E-11 


-6.687E-12 


-5.090E.12 


-1.500E-12 


1.710E-12 


2.142E-11 


-9.880E-12 


-7.400E-12 


-2.110E-12 


_ 


II 


-1.474E-I2 


•1.180E-12 


1.034E-I1 


II 


6.670E-12 


-2.886E-12 


-2.323E-12 


1.035E-11 


■■gynTaiii 

iBSSisaBsnMFai 

imwm^assaaMMRn^ 


1.826E-11 


-6.710E-12 


-5.430E-12 


-1.530E-12 


2.030E-12 


8.010E-12 


1.059E-11 


1.645E-11 


|gl:lY.i:an 


-4.404E-12 


-2.408E-12 


3.442E-12 


1.221E-11 


immsn 


Altitude  400  km 


MedianI  2.128E-12 


-9.610E-13 


-7.730E-13 


II 


-1.994E-12 


-1.613E.12 


saoi 


I.528E-12 


2.479E  12 


7.460E-12 


iKtisaai 
iKuaani 


6.990E-12 


-3.415E-12 


-2.635E-12 


-7.980E-13 


9.390E-13 


3.780E-12 


5.850E-12 


1.045E-1I 


l.a26E-ll 


4.567E-12 


-3J06E-12 


Altitude  425  km 


Mediae 


1 


2.848E-12 


SBsaai 


-1.740E-13 


2.250E-13 


1.052E-12 


1.727E-12 


5.424E.12 


4,788E-12 


-1.129E-12 


-3.660E-13 


4.650E-13 


72S7E-12 


-4.Q22E-12 


9.405E-12 


•3.957E-12 


3.019E-12 


I  fc  y/i*)  *^3  !  a  M  1  :g  KM 


7.000E-13 


3.752E-12 


1.034E-11 


8.152E-12 


5.205E-12 


6.985E-12 


1.167E-11 


1.013E-11 


Table  C-1.  Median  value  of  global  maximum  densities  fw  alutude  and  F10.7  bin  with 
differences  between  global  maximum  density  and  median  value  for  several 
percentile  ranges.  Densities  in  k^m^  (page  ^  of  6). 


Ran&e 

Bin  1 
66-102 

Bin  2 
102-138 

Bin  3 
138-174 

Bin  4 
174-210 

Bin  5 
210-246 

All 

66-246 

Altitude  4 

i50  km 

Median 

8.298E-13 

1.898E-12 

3.320E-12 

5.196E-12 

6.873E-12 

2.021E-12 

1 

-4.173E-13 

-9.744E-13 

.1.801E.12 

-3.020E-12 

.3.058E-12 

-1.570E-12 

5 

.3.393E-13 

.7.960E-13 

-1.408E-12 

-2.339E-12 

-2.510E-12 

-1.446E-12 

33 

-1.171E-13 

-2.620E-13 

-4.380E.13 

-7.130E-13 

-7.360E-13 

67 

1.530E.13 

3.350E-13 

5.250E-13 

7.750E-13 

1.009E-12 

1.339E-12 

95 

7.292E-13 

1.652E-12 

2.183E-12 

3.123E-12 

4.197E-12 

5.321E-12 

99 

1.211E-12 

2.793E-12 

3.445E-12 

4.864E-12 

5.677E-12 

7.854E-12 

100 

3.972E-12 

8.016E-12 

^.387E-12 

9.494E-12 

9.257E-12 

1.411E-11 

Altitude  475  km 


Median 

5.271E-13 

1.277E-12 

2.326E-12 

3.758E-12 

5.073E-12 

1.366E-12 

1 

.2.767E-13 

-1.320E-12 

-2.279E.12 

.2.374E-12 

-1.091E.12 

5 

.2.262E-13 

-5.634E-13 

-1.038E-12 

-1.778E-12 

-1.956E-12 

-l.OlOE-12 

33 

-7.900E-14 

•1.870E-13 

-3.270E-13 

-5.500E-13 

-5.800E-13 

-5.955E-13 

67 

1.045E-13 

2.440E-13 

3.960E-13 

6.040E-13 

8.040E-13 

9.900E-13 

95 

5.079E-13 

1.221E-12 

1.671E-12 

2.469E.12 

3.398E.12 

4.080E-12 

99 

8.539E-13 

2.091E-12 

2.661E-12 

3.883E-12 

4.622E-12 

6.122E-12 

100 

2.927E-12 

6.243E-12 

5.025E-12 

7.732E-12 

7.647E-12 

1.135E.11 

Altitude  500  km 


Median 

3.384E-13 

8.677E-13 

1.644E-12 

2.743E-12 

3.777E-12 

9.322E-13 

1 

-1.841E-13 

-4.856E-13 

-9.711E-13 

-1,728E-12 

•1.850E-12 

-7.616E-13 

5 

-1.512E-13 

.4.009E-13 

-7.689E-13 

-1.358E-12 

-1.531E-12 

-7.087E-13 

33 

-5.350E-14 

-2.450E-13 

-4.260E-13 

-4.590E.13 

-4.257E-13 

67 

7.160E-14 

1.773E-13 

3.000E-13 

4.720E-13 

6.420E-13 

7.348E-13 

95 

3.551E-13 

9.063E.13 

1.285E-12 

1.958E-12 

2.756E-12 

3.142E-12 

99 

6.048E-13 

1.572E-12 

2.065E-12 

3.107E-12 

3.771E-12 

4.793E-12 

100 

2.168E-12 

4.883E-12 

3.970E-12 

6.314E-12 

6.323E-12 

9.168E-12 

Altitude  550  km 


1.445E-13 

4.091E-13 

8.387E-13 

L493E-12 

2.134E-12 

4.442E-13 

1 

-2.445E-13 

-5.309E-13 

-L005E-12 

-1.131E-12 

-3.751E-13 

5 

-2.037E-13 

-4.255E-13 

-9.42E-13 

-3.521E-13 

33 

.2.53E-14 

-7E-14 

-1.39E-13 

-2.57E-13 

-2.85E-13 

-2.197E-13 

67 

1.778E-13 

5.04E-13 

LOME- 12 

L785E-12 

2.554E-12 

8.521E-13 

95 

3.185E-13 

9.147E-13 

1.607E-12 

2.738E-12 

3.966E-12 

L038E-12 

99 

4.503E-13 

1.309E-12 

2.091E-12 

3.5E-12 

4.667E-12 

3.413E-12 

100 

1.348E-12 

3.435E-12 

3.343E-12 

5.744E-12 

5.415E-12 

7.105E-12 

C-4 


Table  C-1.  Median  value  of  global  maximum  densities  for  altitude  and  F10.7  bin  with 
differences  between  global  maximum  density  and  median  value  for  several 
percentile  ranges.  Densities  in  (page  4  of  6). 


i'’l0.7B 

Ranee 


Binl 

66-102 


Altitude  6(X)  km 


Median  6.534E-14 


-3.759E-14 


.3.151E-14 


-1.184E-14 


8.125E-14 


1.52E-13 


2.213E-13 


7.434E-13 


Altitude  650  km 


Median!  3.179E-14 


-1.763E-14 


.1.482E-14 


-5.64E-15 


3.962E-14 


7.566E-14 


1.128E-13 


4.199E-13 


Bin  2 


102-138 


1.995E-13 


-1.247E-13 


-1.048E.13 


-3.68E.14 


2.506E-13 


4.844E-13 


7.191E-13 


2.103E-12 


1.009E-13 


Bin  3 
138-174 


4.4E-13 


-2.935E-13 


•2.384E-13 


5.421E-13 


9.033E-13 


1.214E-12 


2.04E-12 


Bin  4 
174-210 


Bins 

210-246 


All 

66-246 


8.325E-13 

1.243E-12 

2.181E-13 

-5.907E-13 

-7.072E.13 

-1.874E-13 

-4.772E-13 

-5.945E-13 

-1.771E-13 

-1.579E-13 

-1.87E-13 

-1.139E-13 

1.015E-12 

1.512E-12 

4.479E-13 

1.635E-12 

2.468E-12 

5.578E-13 

2.148E-12 

2.955E-12 

2.091E-12 

3.719E-12 

4.267E-12 

5.292E-12 

Altitude  700  km 


Median!  1.696E-14 


1 


5 


33 


67 


95 


99 


100 


-7.289E-15 


-2.78E-15 


2.088E-14 


3.961E-14 


5.967E-14 


2.413E-13 


-5.461E-14 


-1.966E-14 


1.289E-13 


2.626E-13 


4.049E-13 


1.313E-12 


5.319E-14 


.3.39E.14 


-2.894E-14 


-1.063E-14 


6.859E-14 


1.456E-13 


2.32E-13 


8.323E.13 


2.366E-13 

4.755E.13 

7.37E-13 

l.llE-13 

-1.638E-13 

-3.515E-13 

-4.436E-13 

-7.215E-13 

-1.345E-13 

-1877E-13 

.3.755E.13 

-9.077E-14 

-4.63E-14 

-9.81E-14 

-1.205E-13 

-5.992E-14 

2.971E-13 

5.907E-13 

9.136E-13 

2.411E-13 

5.199E-13 

9.953E-13 

1.909E-12 

3.083E^13 

7.186E-13 

1.343E-12 

2.041E-12 

1.304E-12 

1.269E-12 

2.43E-12 

2.821E-12 

3.447E-12 

1302E-13 


-9.212E.14 


-7.638E-14 


-2.68E-14 


1.664E-13 


3.041E-13 


4326E-13 


8.026E-13 


2.76E-13 


-2.101E-13 


-1.741E-13 


-6.09E-14 


3.495E-13 


6.16E-13 


8.539E-13 


1.647E-12 


4.445E-13 


-2.804E-13 


-2.391E-13 


.7.8IE-I4 


5.617E-13 


1.008E-12 


1.2S3E-12 


1.941E-12 


5.872E-14 


-4.976E-14 


-4.746E-14 


-3.199E-14 


1.329E-13 


1.752E-13 


8.267E-13 


2.326E-12 


Altitude  750  km  _ 


Table  C>1.  Median  value  of  global  maximum  densities  for  altitude  and  Fioj  bin  with 
differences  between  global  maximum  density  and  median  value  for  several 
percentile  ranges.  Densities  in  kg/m^  (page  S  of  6). 


C-5 


ri0.7B 

Range 

Bin  1 
66-102 

Bin  2 
102-138 

Bin  3 
138-174 

Bin  4 
174-210 

Bin  5 
210-246 

All 

66-246 

Altitude  S 

EOO  km 

Median 

6.63  IE- 15 

1.758E-14 

4.354E-14 

9.912E-14 

1.7E-13 

1.931E-14 

1 

-2.69E-15 

-1.025E-14 

-3.061E-14 

-7.752E-14 

-1.142E-13 

-1.514E-14 

5 

-2.232E-15 

-8.788E-15 

-2.576E-14 

-6.551E-14 

•9.864E-14 

-1.441E-14 

33 

-8.29E-16 

-3.3E-15 

-9.44E-15 

-2,401E.14 

-3.38E-14 

-9.775E-15 

67 

7.803E.15 

2.248E-14 

5.663E-14 

1.292E-13 

2.221E-13 

4.449E-1.4 

95 

1.339E-14 

4.905E-14 

1.105E-13 

2.475E-13 

4.363E-13 

6.268E-14 

99 

1.961E-14 

8.168E-14 

1.647E-13 

3.606E-13 

5.612E-13 

3.477E-13 

100 

8.534E-14 

3.503E-13 

3.361E-13 

7.713E-13 

9.331E-13 

1.084E-12 

Altitude  850  km 


Median 

4.735E-15 

1.122E-14 

2.67E-14 

6.153E-14 

1.079E-13 

1.223E-14 

1 

-1.747E-15 

-6.052E-15 

-1.821E-14 

-4.795E-14 

-7.37  IE- 14 

-9.087E-15 

5 

-1.434E-15 

-5.17E-15 

I  -1.536E-14 

-4.082E-14 

-6.401E-14 

-8.594E-15 

33 

-5.19E-16 

-1.941E-15 

-5.69E-15 

-1.529E-14 

-2.228E-14 

-5.737E-15 

67 

5.454E-15 

1.41E-14 

3.47  IE- 14 

8.102E-14 

1.43E-13 

2.727E-14 

95 

8.761E-15 

3.004E-14 

6.885E-14 

1.601E-13 

2.927E-13 

3.99E-14 

99 

1.24E-14 

5.04E-14 

1.044E-13 

2.388E-13 

3.824E-13 

2.297E-13 

100 

5.272E-14 

2.316E-13 

2.216E-13 

5.322E-13 

6.51E-13 

7.466E-13 

Altitude  900  km 


3.576E-15 

7.67  IE- 15 

1.713E-14 

3.916E-14 

6.967E-14 

8.284E-15 

1 

-1.243E-15 

-3.81E-15 

-1.115E-14 

-3.006E-14 

-4.788E-14 

-5.837E-15 

5 

-1.014E-15 

-3.236E-15 

-9.388E-15 

-2.571E-14 

-4.176E-14 

-5.478E-15 

33 

-3.59E-16 

-1.1%E-15 

-3.5E-15 

-9.75E-15 

-1.47E-14 

-3.565E-15 

67 

4.048E-15 

9.423E-15 

2.211E-14 

5.194E-14 

9.359E-i4 

1.748E-14 

95 

6.152E-15 

1.92E-14 

4.395E-14 

1.052E-13 

1.986E-13 

2.65E-14 

99 

8.389E-15 

3.206E-14 

6.737E-14 

1.6E-13 

2.635E-13 

1.535E-13 

100 

3.354E-14 

1.548E-13 

1.478E-13 

3.773E-13 

4.681E-13 

5.294E-13 

AJtitude  950  km 


2.804E-15 

5.518E-15 

1.16E-14 

2.578E-14 

4.609E-14 

5.977E-15 

1 

-9.44E-16 

-2.505E-15 

-7.122E-15 

-1.927E-14 

-3.152E-14 

-4.028E-15 

5 

-7.66E-16 

-2.09E-15 

-5,973E-15 

-1.652E-14 

-2.762E-14 

-3.75E-15 

33 

-2.67E-16 

-7.17E-16 

-2.228E-15 

-6.35E-15 

-9.88E-15 

-2.352E-15 

67 

3.15E-15 

6.711E-15 

1.477E-14 

3.419E-14 

6.233E-14 

1.182E-14 

95 

4.589E-15 

1.292E-14 

2.892E-14 

7.044E-14 

1.363E-13 

1.834E-14 

99 

6.044E-15 

2.114E-14 

4.453E-14 

1.089E-13 

1.833E-13 

1.043E-13 

100 

2.21  IE- 14 

1.053E-13 

1.003E-13 

2.654E-13 

3.33E-13 

3.731E-13 

^****^5*'  of  global  maximum  densities  for  altitude  and  F.aibin  with 

differences  between  global  maximum  density  and  median  value  for  22^Jeral 
Dcrcentile  npncirtac / _ ^  venu 


™®  of  achieving  a  time  interval  without  encountering  a 

thermospheric  density  level  above  a  given  pereentile  value. 


13-month 
Smoothed  Fio,7 
Ranee 

Percentile  Level  1 

66  Percent 

90  Percent 

95  Percent 

99  Percent 

66  to  102 

0.13 

0.46 

0.63 

0.86 

102  to  138 

0.21 

0.55 

0.71 

0.90 

138  to  174 

I  0.28 

0.55 

0.68 

0.89 

174  to  210 

0.35 

0.60 

0.75 

0,93 

210  to  246 

0.38 

0.65 

0.76 

0.87 

1-2.  PtobaltiliQr  of  going  at  least  30 
(Random  Start) 


days  without  exceeding  the  indicated  percentile  density  level 
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Figure  C-1.  ^bability  of  meeting  or  exceeding  a  given  time  interval  without  exceedine  the  95th 
percen^e  dei^ty  (top)  or  99th  percentile  density  (bottom).  The  bins  represent  a  range 
of  F10.7:  (1)  66-102.  (2)  102-138.  (3)  138-174.  (4)  174-210,  and  (5)  210-246?^^ 
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